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PREFACE 

The Landsai Thematic Mapper (TM) is an advanced 
remote sensing system designed to measure the In- 
tensity of Earth radiation in selected portions of the 
electromagnetic spectrum. The purpose of this docu- 
ment is to identify prospective research themes in 
various aspects of the Earth sciences that can be ad- 
dressed in new and innovative ways through the 
analysis of multispectral TM Imagery. This prospec- 
tus specifically describes the measurement capabilities 
of tlie TM, the various types of Earth information that 
can be obtained through the analysis of TM imagery, 
and potential applications of such information to con- 
temporary research issues in various Earth science 
disciplines. 

The sections of this document dealing with the 
measurement capabilities of the TM and the orbital 
performance of the TM sensor system were prepared 
by the Landsat Project Science Office at NASA’s God- 
dard Space Flight Center. The Information presented 
in this portion of the document is based in part upon 
the results of the Landsat Image Data Quality Analysis 
(LIDQA) program that was inidated by NASA in 1982, 
The purpose of the LIDQA program is to determine 
the radiometric and geometric characteristics of TM 
image data in both raw and processed form, and to 
compare the observed characteristics of orbital TM 
imagery with technical specifications for image quality 
that were established prior to the launch of the TM, 
The investigators participating in this program are as 
follows; 

James Anderson, NASA-National Space Technology 

Laboratories 

Paul Anuta, Purdue University 

G, Begui, Centre Nationale d’Etudes Spatiale (France) 

Lee Bender, U. S. Geological Survey 

Ralph Bernstein, IBM Palo Alto Scientific Center 

Eric Beyer, General Electric Space Division 

Robert Colwell, University of California at Berkeley 

Jeffrey Dozier, University of California at Santa 

Barbara 

Michael Duggin, State University of New York at 

Syracuse 

Jack Engel, Hughes Sasua Barbara Research Center 
Jon Erickson, NASA— Johnson Space Center 
John Everett, Earth Satellite Corporation 
Gary Ford, University of California at Davis 
Luigi Fusco, European Space Agency Earthnet Pro- 
gram Office 


Charlotte Gurney, Systems and Applied Science 
Corporation 

Charles Hill, NASA— National Space Technology 
Laboratories 

Warren Hovis, NOAA— National Environmental 
Satellite Service 

Michael Jackson, National Environmental Research 
Council (England) 

Hugh Kieffer, U. S, Geological Survey 

Vic Klemas, University of Delaware 

Harold Lang, NASA— Jet Propulsion Laboratory 

Don Lauer, U, S, Geological Survey 

Robert MacDonald, NASA— Johnson Space Center 

William Mallla, Environmental Research Institute of 
Michigan 

Nelson Parada, Instituto de Pesquisas Espaciais (Brazil) 
John Price, U, S, Department of Agriculture 
S. I, Rasool, Ecole Normale Superieure (France) 
John Schott, Rochester Institute of Technology 
Philip Slater, University of Arizona 

W, Murray Strome, Canada Center for Remote Sens- 
ing (Canada) 

June Thormodsgard, U, S. Geological Survey 
Roy Welch, University of Georgia 
Robert Wrigley, NASA — Ames Research Center 
George Wubelic, Battelle Pacific Northwest 
Laboratories 

Albert Zobrist, NASA— Jet Propulsion Laboratory 

LIDQA investigators are currently analyzing data ac- 
quired by the TM sensor systems on both the Land- 
sat 4 and 5 spacecraft. This document summarizes the 
collective results of these ongoing studies as of June, 
1984. The sections of this document dealing with pro- 
spective scientific applications of TM imagery sum- 
marize the findings and recommendations of the TM 
Science Working Group (TMSWG). The TMSWG was 
organized by NASA during 1983 to define innovative 
applications of the TM’s measurement capabilities for 
basic research in the Earth sciences. The TMSWG was 
composed of the following individuals: 

Mark Settle, Co-Chairman, NASA Headquarters 
Vince Salomonson, Co-Chairman, NASA — Goddard 
Space Flight Center 

Jim Irons, Executive Secretary, NASA — Goddard 
Space Flight Center 

James Anderson, NASA— National Space Technology 
Laboratories 
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Don Peering, NASA— -Godclarcl Space riight Center 
Diane Evans, NASA— Jet Propulsion Laboratory 
Harold Lang, NASA— Jet Propulsion Laboratory 

Ken Langren, NASA— National Space Technology 
Laboratories 

M. C. Trlchel, NASA— Johnson Space Center 
Robert Wrlgley, NASA— Ames Research Center 

TMSWG participants organized a series of 
workshops during the spring and summer of 1 983 to 
identify research issues in different disciplinary fields 
that could be profitably addressed through the 
analysis and Interpretation of TM imagery. The In- 
dividuals who participated in these workshops are 
listed below, 

Michael Abrams, NASA— Jet Propulsion Laboratory 
Gautam Badhwar, NASA— Johnson Space Center 
Bryan Bailey, U. S. Geological Survey 
Robert Bizzell, NASA— Johnson Space Center 
Paul Bock, University of Connecticut 
William Bull, University of Arizona 
William Clbula, NASA— National Space Technology 
Laboratories 

Carolyn Clark, Lockheed Engineering and Manage- 
ment Services Company 
James Gone!, NASA— Jet Propulsion Laboratory 
B, Dey, Howard University 
Jack Estes, University of California at Santa Barbara 
Gary Ford, University of California at Davis 
James Foster, NASA— Goddard Space Flight Center 
Tom George, Camp Dresser &. McKee, Incorporated 
Jean Gervin, NASA — Goddard Space Flight Center 
Alan Gillespie, NASA— Jet Propulsion Laboratory 
Edward Guinness, Washington University 
Robert Gurney, NASA— Goddard Space Flight Center 
Robert Haas, U. S. Geological Survey 
Keith Henderson, NASA— Johnson Space Center 
John Hill, Louisiana State University 
Thomas Jackson, U. S. Department of Agriculture 

Chris Justice, NASA — Goddard Space Flight Center 
Edward Kanemasu, Kansas State University 

Richard Kauth, Environmental Research Institute of 
Michigan 


Siamak Khorram, North Carolina State University 
Vic Klemas, University of Delaware 
Anthony Lewis, Oregon State University 
James Lewis, South Dakota State University 
Stuart Marsh, Sun Exploration Company 
John McKeon, Gulf Science & Technology Company 
Elizabeth Middleton, NASA— Goddard Space Flight 
Center 

Bradley Muslck, NASA— National Space Technology 
Laboratories 

Eugene Peck, Hydex Corporation 
Gary Petersen, Pennsylvania State University 
Russell Pettit, Texas Tech University 
David Pitts, NASA— Johnson Space Center 
Mel Podwysocki, U. S. Geological Survey 
Charles Poulton, Poulton Associates 
Robert Ragan, University of Maryland 

Doug Rickman, NASA— National Space Technology 
Laboratories 

Ray Sadowski, Amax Exploration, Incorporated 
Gary Shelton, Environmental Protection Agency 
Ray Sugiura, Woodward-Clyde Consultants 
Don Thomas, Dames & Moore 
David Thompson, NASA— Johnson Space Center 
Paul Tueller, University of Nevada at Reno 
Joseph Ulliman, University of Idaho 
Steve Wharton, NASA — Goddard Space Flight Center 
Darrel Williams, NASA— Goddard Space Flight Center 
Chuck Wood, NASA— Johnson Space Center 


The contribution of all of the individuals listed 
above to the preparation of this document is greatly 
appreciated. Special thanks go to Jim Irons (Goddard 
Space Flight Center) and Kristine Butera (NASA Head- 
quarters) for their assistance in the finai preparation 
of the report. 

Mark Settle Vincent Salomonson 

Landsat Program Landsat Project Scientist 

Scientist 


SECTION I 
INTRODUCTION 


The surface of the Earth is a fragile and complex 
Interface between the solid Interior and gaseous ex- 
terior of our planet. It displays a far wider range of 
cover materials and conditions than any other 
planetary surface found within our solar system. Fur- 
thermore, it is in a continual state of modification and 
is affected by a far wider range of processes than any 
other planetary surface known to exlsr today. 

The Earth’s surface provides a natural home for 
mankind, and it is the source of both the renewable 
and non-renewable resources needed to support our 
modern human society. At the moi.t fundamental 
level, the study of the Earth’s surface is motivated by 
our desire as a species to learn more about the pro- 
cesses that constructed the home we live in today, 
and the processes that may alter or modify the 
habitability of our home in the future. A phenomeno- 
logical understanding of the history of our planet and 
the processes that govern its future evolution is essen- 
tial for our long term survival. 

Much of the current research in the Earth sciences 
is concerned with forecasting future conditions on 
our planet. These contemporary research efforts have 
revealed that the Earth’s oceans, atmosphere, and land 
surface are coupled in many different ways. Subtle 
changes in such observable parameters as sea surface 
temperature, atmospheric chemistry, and vegetation 
distribution can be empirically correlated, suggesting 
that they may be interrelated in a phenomenological 
sense. In many instances, however, our understand- 
ing of the underlying physical phenomena that pro- 
duce correlated changes in surface characteristics is 
quite rudimentary, In many other instances such 
understanding is lacking altogether. It is generally 
quite difficult to forecast spatial and temporal changes 
in terrestrial surface conditions on the basis of our cur- 
rent understanding of terrestrial surface processes. 

Space technology provides a unique capability to 
“step back” from the Earth’s surface in both a physical 
and intellectual sense. Remote sensing systems 


mounted on Earth-orbiting spacecraft can be used to 
observe surface conditions on our planet at a uniform 
scale and format on a globally repetitive basis, This 
capability provides a means of studying and 
characterizing the full range of surface conditions 
found on the Earth in a way that has not been possi- 
ble in the past, Furthermore, this capability can be 
used to monitor temporal variations in surface con- 
ditions or to compare the characteristics of different 
regions as a means of obtaining improved insight in- 
to the physical processes that modify the Earth’s 
surface. 

As the name implies, the Landsat series of spacecraft 
have been specifically designed to conduct orbital 
observations of the Earth’s land surface. The Landsat 
program began in earnest in July, 1972, with the 
launch of Landsat 1 . The principal Earth observation 
sensor on this spacecraft was the Multlspectral Scan- 
ner (MSS). The MSS was specifically designed to 
measure the intensity of Earth radiation in the visible 
and infrared portions of the electromagnetic spectrum 
at a spatial resolution that was appropriate for detailed 
observations of land conditions. Landsats 2 and 3 car- 
ried the same sensor complement as Landsat 1 , and 
they were launched in 1975 and 1978, respectively. 
A second generation sen.sor system capable of con- 
ducting multispectral Earth surveys was recently 
developed and initially launched on Landsat 4 in July, 
1982. This sensor is known as the Thematic Mapper 
(TM). It is superior to the MSS in terms of its spectral 
range and resolution, its spatial resolution, and its 
radiometric sensitivity and accuracy. The technical 
specifications of the TM are presented in detail in Sec- 
tion II of this document, 

Landsat research activities conducted during the 
past decade have passed through several phases. Dur- 
ing the mid-1970’s, the main objective of these ac- 
tivities was to determine the various types of infor- 
mation that could potentially be extracted from 
multispectral Earth imagery acquired at a synoptic 
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scale. A wide variety of projects were conductea In 
which spatial and temporal variations in the spectral 
characteristics of the Earth’s surface were rciated to 
ground-based measurements of surface properties. 
These research efforts lead to the development of new 
techniques for the reduction and analysis of 
multlspectral Earth Imagery. They also established a 
rigorous scientific basis for deriving specific types of 
surflcial information from such imagery. 

Landsat research activities conducted during the lat- 
ter 1970’s were primarily designed to identify prac- 
tical applications of Landsat data The main objective 
of these activities was to evaluate the utility of Land- 
sat Imagery for resource management and monitor- 
ing. A variety of large scale research projects were 
conducted by NASA in collaboration with other 
government agencies and private companies. The 
results of these research efforts indicated that Earth 
information derived from Landsat imagery could con- 
tribute significantly to improved methods of crop 
forecasting, mineral exploration, spring runoff 
forecasting, flood damage assessment, etc. 

The collective results of the first ten years of Land- 
sat research have provided extensive insight into the 
utility of multispectral surveys for specific types of 
Earth observations. These results have also indicated 
that orbital observations of the Earth’s spectral 
characteristics can be applied to the study of chang- 
ing surface conditions. This historical experience has 
enabled Earth scientists to develop considerable in- 
tuition concerning the fundamental capabilities and 
limitations of orbital multispectral surveys for studies 
of the Earth. This experience, in combination with 
the advanced measurement capabilities of the TM, 
creates a unique opportunity to apply orbital 
multispectral surveys to the study of terrestrial sur- 
face conditions and processes. 

The purpose of this prospectus is to identify topical 
problems in different aspects of the Earth sciences that 
could be addressed in new and innovative ways 


through the analysis of TM imagery. These prospec- 
tive research themes are defined here In a general 
sense in relation to the technical measurement 
capabilities of the TM and the various types of Earth 
information that can potentially be derived from 
multispectral TM imagery. 

This document is organized Into two parts. Section 
II provides an overview of the total system that has 
been developed to acquire and reduce TM data. The 
technical capabilities ohhis system are presented here 
in some detail. In addition, the orbital performance 
of the TM sensor is described, based upon the analysis 
of Landsat 4 and 5 TM data that have been collected 
to date. Section III of this document outlines prospec- 
tive research themes that could potentially be ad- 
dressed through the analysis and Interpretation of TM 
imagery. This discussion Is not intended to be a com- 
prehensive listing of all of the research topics that are 
amenable to TM Image analysis, nor does it reflect a 
prioritization of research topics in the Earth sciences 
on the part of NASA. The Intent of this discussion is 
simply to illustrate generic types of scientific Issues 
that could be addressed in new and Innovative ways 
through the analysis of multispectral TM imagery. Sec- 
tion III has been subdivided into a discussion of pro- 
spective research themes in the fields of vegetation 
and soils science, geology, and hydrology. This 
organization was adopted to reflect the major classes 
of natural cover materials encountered on the Earth’s 
surface. The subdivision of this discussion into these 
categories is nor intended to preclude the application 
of TM data to the study of topical problems in related 
disciplinary fields. 

In summary, it is anticipated that the information 
provided in this prospectus will assist Earth scientists 
in identifying potential scientific applications of TM 
data and in devising experimental methodologies that 
make substantive use of TM imagery for basic research 
in the Earth sciences. 



SECTION II 

CAPABILITIES AND PERFORMANCE 
OF THE TM SYSTEM 


INTRODUCTION 


The first Lancisat satellite, originally called the Earth 
Resources Technology Satellite (ERTS), was launched 
on July 23, 1972 by the National Aeronautics and 
Space Administration (NASA). Landsat 1 was designed 
to acquire Earth imagery in digital format for a 
minimum period of one year. Data from the satellite's 
Multispectral Scanner (MSS) proved so useful that the 
satellite remained in operation until the MSS failed in 
January 1978. Four additional satellites dedicated to 
Earth observations have been launched by NASA 
subsequent to the launch of Landsat 1 : Landsat 2 on 
January 22, 1975; Landsat 3 on March 4, 1978; Land- 
sat 4 on July 16, 1982; and Landsat 5 on March 1, 
1984, 

The Landsat satellites greatly enhanced our ability 
to study the Earth by providing multispectral Earth 
imagery at a synoptic scale on a globally repetitive 
basis. The first three satellites, which were quite 
similar in form and operation, all carried the MSS as 
their primary sensor payload. These satellites were all 
launched as experimental missions, but image analysis 
techniques developed to the point where MSS data 
arc now used in an operational manner for mineral 
and oil exploration, crop monitoring and yield predic- 
tion, forest and range management, map making and 
many other applications. Ground stations which 
directly receive and process MSS data nave been in- 
stalled in eleven nations to permit worldwide access 
to the data. 

The value of MSS data was developed and demon- 
strated by extensive research, but this research also 
revealed limitations in data utility. The potential for 
such limitations was recognized before the launch of 
Landsat 1. Planning for a second generation Earth 
resource observation system to follow the initial Land- 
sat series was initiated by Ni1.SA in 1970. The con- 
figuration of a system which could enhance Earth 
observation capabilities evolved from a subsequent 


progression of study efforts and advisory groups 
representing the many facets of the remote sensing 
community. The satellite portion of the enhanced 
system was launched in July of 1982 and Is now culled 
Landsat 4, The launch of a nearly Identical satellite, 
Landsat 5, followed In March, 1984. The Landsat 4 
and 5 sensor payload consists of an MSS similar to the 
instruments carried onboard the initial three Landsat 
satellites and a new sensor, the Thematic Mapper 
(TM). 

Landsat 4 and 5 offer several technological 
refinements over the previous systems. The platform 
portion of the satellites provides superior attitude con- 
trol when compared to the first three Landsat 
satellites. These initial systems also relied upon tape 
recorders to store sensor data for subsequent 
transmission to a ground receiving station. These tape 
recorders were often the first system component to 
fail. Landsat 4 and 5 do not carry tape recorders, but 
instead transmit data directly to ground receiving sta- 
tions as the data are acquired. 

The TM and the MSS are multispectral scanners 
which both function in a similar manner. The proven 
technology of the MSS, however, has been refined to 
create improved sensor characteristics for the TM. 
These improvements Include a finer spatial resolution, 
additional spectral bands located in new regions of 
the spectrum, and a greater number of data quantiza- 
tion levels to take advantage of an enhanced radio- 
metric sensitivity. These sensor attributes are intended 
to increase both the quantity and quality of Earth in- 
formation that can potentially be derived from Land- 
sat imagery. 

Knowledge of satellite and sensor characteristics is 
essential to the effective utilization of TM data. This 
section provides a detailed description of the Land- 
sat 4 and 5 satellites, their instrument payloads, their 
orbits, and the reception, processing, and dissemina- 
tion of TM data on the ground. Information on or- 
bital sensor performance and TM data quality is also 
presented. 
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LANDSAT 4 AND 4 SPACECRAFT 
CHARACTERISTICS 

Landsai -4 and S arc (he flrM saicllitc.s in the Land- 
sat series to conform to NASA's standard design for 
a Multimission Modular Spacecraft (MMS). Each 
satellite consists of two connected frames (Figure II- 1 ) 
One frame, the universal MMS structure, suppHirts 
modules for satellite propulsion, attitude control, 
power, and data communications. The other frame, 
the mi.ssion-specific lastrument Module Structure, car- 
ries the instrument and sensor payload. For Landsat 
-4 and S, a single-w'ing array of solar panels and a large 
antenna mxst are attached to this supporting structure. 
The antenna is used for data communications with 
the Tracking and Data Relay Satellite System (TDRSS) 


and (he mxst minimiaes interference hy the solar panel 
array. The solar array consists of four I .S m hy 2.3 m 
panels which can collectively generate up to 22(N) 
watts of power. 

The MMS offers improved attitude control 
capabilities over the previous Landsat satellites. The 
MMS permits three-axLs control of attitude. The point- 
ing accuracy is specified to be within 0.01 degree (one 
sigma) and the stability is specified to he within I O'** 
degree/second (one sigma). The first three Land.sat 
satellites only provide one-axis control of attitude 
resulting in a 0.7 degree pointing accuracy and a 
stability value oi' 0 t)l degree/second (Salomonson et 
al., I OHO). The improved attitude control facilitates 
the cartographic and multitemporal registration of im- 
age data. 
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Fifiure II- 1 Exploded view of the Landsat 4 and 5 Multimission Modular Spacecraft. 


4 


OHKill^AL PAG£ .9 
OF POOR QUALITY 




0R/G!NAL page 
OF POOR QUALITY 


SENSOR PAYLOAD 
CHARACTERISTICS 


The Landsat 4 and 3 sensor payload consists of the 
Multispectral Scanner (MSS) and the Thematic Map- 
per (TM). Both sensors are object-space line scanners 
which use moving mirror assemblies to scan across 
the spacecraft ground track and the orbital motion of 
the spacecraft to scan the perpendicular (l.e„ along- 
track) direction. Both sensors generate multispectral 
Image data in a digital format. 

The MSS sensors on Landsat 4 and 5 are nearly iden- 
tical to the MSS sensors on previous Landsat satellites. 
The MSS mirror reflects Earth radiation through a 
Ritchey-Chretein telescope onto a single focal plane. 
At the focal plane, optical filters separate the incident 


electromagnetic energy Into the familiar four MSS 
spectral bands (Table IM). Fiber optic bundles 
transmit the separated energy from the focal plane to 
detectors which convert the energy to electronic 
signals (Figure II-2). A shutter wheel In front of the 
fiber optics bundle Is Involved In the radiometric 
calibration of MSS data. Six detectors are used for each 
spectral band; photomultiplier tubes are used for spec- 
tral bands M5S1 through MSS3 and silicon photodl- 
s des are employed for band MSS4. The detector, s’ 
bignals arc quantized to ,slx bits (64 digital counts) 
before the Image data are transmitted to the ground. 
With six detectors per band, the MSS generates six 
scan lines of data per scan mirror sweep. Data are ac- 
tively acquired as the scan mirror moves the sensor’s 
field of view from west to east. No data are acquired 
during the return scan (Le., east-to-west), 
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Figure II-2 Schematic view of the Landsat 4 and 5 
Multispectral Scanner, The cutaway diagram in the 
upper left corner shows how Earth radiation enters the 
sensor and is focused onto the optical fibers situated 
in the focal plane, The arrangement of the optical fiber 


endings is displayed schematically in the upper right 
comer of the diagram. Note that six scan lines of image 
data are generated simultaneously every time the scan 
mirror sweeps from west to east. 


Tablo tl-l 

Mulltlspectral Scanner Spectral Bands (Salomonsoii c4 ai„ 1980) 


Band 

Designation 


Specified 
Band Width (/tm) 


Specified 
Radiometric 
Sensitivity (NEAo)‘ 


MSS1 0.5-0.6 0.67% 

MSS2 0.6-0.7 0.57% 

MSS3 0.7— 0.8 0.65% 

MSS4 0.8-1. 1 0.70% 

* NEAc: noise equivalent change In surface reflectance. 


The across-track swath (185 km) and the 
instantaneous-field-of-vlew (80 m) of the new MvSS 
sensors are nearly equal to the swath width and IFOV 
of the earlier MSS’s. Since the earlier Landsat satellites 
all flew at a higher altitude (920 km) than Landsat 4 
and 5 (705 km), the optics and the active scan angf»' 
of the new MSS sensors were adjusted to compensaF^’ 
for the difference in altitude. The earlier sensors col- 
lected data over a total scan angle of 1 1 .56 " wliereas 
the new MSS sensors employ a scan angle of 14.92 °. 

Several major differences exist between the basic 
design of the MSS and the design of the TM. First, the 
TM acquires data during both the forward (west-to- 
east) and reverse (east-to-west) sweeps of its scan mir- 
ror. This bidirectional approach was adopted to 
reduce the rate of oscillation of the scan mirror and 
increase the dwell time of individual detectors upon 
the Earth’s surface. The TM actively scans the surface 
through an angle of 15.4® (plus-or-minus 7.7® 
measured from nadir), It completes approximately 
seven complete scan cycles per second (note that each 
complete cycle consists of a forward and reverse 
scan). 

As another example of design differences, the TM 
does not rely on fiber optics, as does the MSS, to 
direct incoming electromagnetic energy onto 
photosensitive detectors. Fiber optics are incapable 
of transmitting 100 percent of the incident energy 
from the focal plane to the sensor’s detectors. Instead 
of using fiber optics, the TM focuses incident radia- 
tion through its own Ritchey-Chretein telescope 
directly onto detectors within a prime focal plane 
assembly, and, through a set of relay mirrors, to the 
detectors of a cooled focal plane assembly (Figure 
II-3). The prime focal plane assembly contains four 
sets of 16 monolithic silicon detectors for four spec- 
tral bands In the visible and near infrared portions of 


the spectrum. The 16 detectors for each band are 
mounted In a common module with a spectral filter 
and preamplifier electronics. The cooled focal plane 
assembly contains two arrays of 16 indium an- 
tlmonide (InSb) photodiodes for two middle Infrared 
spectral bands and four mercury cadmium tellurlde 
(HgCdTe) detectors for a thermal band along with 
spectral filters and electronics. The cooled focal plane 
assembly is mounted In a dewar on the second stage 
of a passive radiative cooler. The InSb and HgCdTe 
detectors are mounted on this cooled focal plane to 
increase their radiometric sensitivity. The use of multi- 
ple detectors for each spectral band results In the 
generation of 16 scan lines of data for the six reflec- 
tive bands and four lines of data for the thermal band 
during each sweep of the scan mirror. At any one in- 
stant, surface radiance is sensed by a total of 100 TM 
detectors. 

Due to the spatial separation of the individual detec- 
tor arrays situated within the TM focal planes, an area 
on the Earth’s surface is not simultaneously scanned 
in all seven TM spectral bands. Figure II-3 illustrates 
the projection of the detector IFOV’s onto-the Earth 
surface. Accurate band-to-band registration requires 
precise knowledge of the ground projection of the 
detectors in the two TM focal planes as a function of 
time.Thls information is derived from data concern- 
ing spacecraft position and attitude, the relative posi- 
tion of the individual detector arrays with respect to 
the sensor's optical axis, and the motion of the scan 
mirror during successive scan cycles. A scan angle 
monitor on the mirror generates signals indicating the 
scan mirror’s angular position as a function of time. 
These signals are called mirror scan correction data 
and they are transmitted to the ground for geometric 
processing of TM image data. 
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ffgure 11-3 Schematic view of optical path of Earth 
radiation whhin the Thematic Mapper. Note that the 
detectors for the TM spectral bands are mounted in 
two spatially separated focal planes. The diagram in 


the upper left hand corner displays the instantaneous 
ground projection of the one hundred detectors within 
the TM sensor. 


Signals from the scan angle monitor also Initiate the 
motions of the scan angle corrector located In front 
of the primary focal plane (Figure II-3). During each 
scan mirror sweep, the scan line corrector rotates the 
TM line*of-slght backward along the Landsat ground 
track to compensate for the forward motion of the 
satellite. The compensation prevents the overlap and 
underlap of scan lines and produces straight scan lines 
which are perpendicular to the ground track (Figure 
II-4). Also located on the TM is a three-axis angular 
displacement sensor assembly. The assembly contains 
three angular displacement sensors oriented along or- 
thogonal axes. These sensors provide data to the 
ground data processing system to compensate for any 
high-frequency displacements of the 'TM optical axis 
caused by the vibrations of the TM and MSS scan 
mirrors. 

The TM also contains an internal radiometric 
calibration source located behind the primary mirror 


and in front of the primary focal plane. The calibra- 
tion source consists of three small tungsten filament 
lamps, a blackbody for the thermal band, and a flex 
pivot mounted resonant shutter. The shutter passes 
through the field of view of the TM detectors each 
time the scan mirror changes directions. The shutter 
introduces light from the lamps directly into the field 
of view of the detectors of the reflective bands (TMl 
to TM5 and TM7), The radiant energy from the ther- 
mal calibration standard is reflected into the field of 
view of the thermal band detectors from a mirror 
mounted on the oscillating shutter. The detectors in 
all seven bands also view a dark isothermal portion 
of the shutter which is used to restore the direct cu-- 
rent voltage value of each detector to a preset 
threshold value prior to the next active sweep of the 
scan mirror. This restoration of voltage is referred to 
as the "dc restore cycle” that precedes each period 
of active data collection (Engel, 1980). 
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Fi^iure //•'/ Schematic diagram showing how the Landsat spacecraft. The scan line corrector ensures 
motion of the scan line corrector within the TM com- that individual nt scan lines are oriented perpen- 
pensatesfor the overlap and underlap of adjacent scan dtcular to the spacecraft *s ground track, 

line sweeps that result from the orbital motion of the 


The three tungsten lamps operate at approximately 
the same Intensity, but different neutral density filters 
In front of each lamp cause a difference In the inten- 
sity of light from each lamp for calibration. All possi- 
ble combinations of two lamp states (on and off) for 
three lamps provide eight different light levels for 
detector calibration. These elgm levels include a dark 
level consisting of all three lamps in the "off" state. 
The light level remains constant for approximately 20 
scan cycles (l.e,, 20 successive forward and reverse 
sweeps of the scan mirror). After each 20 cycle inter- 
val, the state of one lamp is changed (l.e., turned on 
or off) to alter the llpht level. Since 16 detectors are 
used for each band, i>. r .ial of 5120 scan lines of Image 
data are p^^nerate^ di^iuig the time Interval required 
to view aH eight fight levels (l.e., corresponding to 
eight calibration levels times 40 scan mirror sweeps 
per level times 16 lines of data per mirror per sweep). 
For the thermal band, the temperature of the 
blackbody and the temperature of the shutter’s dc 
restore surface are used for a two-level calibration of 
the four thermal detectors. 

Figure II-5 provides an overview of the Thematic 
Mapper design. The TM Is horizontally mounted in 
the satellite with the sun shade pointing toward Earth. 
The scan mirror is located directly above the sun 
shade aperture along with the mirror’s drive 
mechanism, control electronics, and scan angle 
monitor. The primary mirror is secured approxi- 
mately In the center of the telescope and is preceded 
by the secondary mirror and optical baffles. The scan 
line corrector, the internal calibrator, and the primary 
focal plane are directly behind the primary mirror. 
The angular displacement .sensors are located within 
a small assembly mounted to the exterior of the TM 
telescope (the assembly is hidden In Figure II-5). The 
back end of the instrument contains the relay optics 
an ] the cooled focal plane within the radiative cooler. 
The electronics package sits in a wedge-shaped box 


above the telescope and contains the multiplexer, 
power supplies, signal amplifiers, and electronic 
filters. 

The complex design of the TM affords advanced 
sensor attributes In comparison to the MSS. For In- 
stance, the TM acquires data In seven spectral bands. 
Four of these bands are located In portions of the 
spectrum not sensed by the MSS. The location and 
width of the seven bands were carefully chosen for 
sensitivity to certain natural phenomena and to 
minimize the attenuation of surface energy oy at- 
mospheric water. Table II-2 lists the spectral bands 
along with the rationale for selecting each band. 

Another advanced feature of the TM Is Its spatial 
resolution. The TM Instantaneous-field-of-vlew(IFOV) 
is 30 m-by-30 m for the six reflective spectral bands 
(bands TMl to TM5 and TM7) and 120 m-by-120 m 
for the thermal band (band TM6). The MSS provides 
a resolution of approximately 80 m-by-80 m In all 
bands. Table II-3 lists the optical characteristics of the 
TM which govern its surface spatial resolution. 
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Figure 11-4 Schematic diagram showing how the 
motion of the scan line corrector within the TM com- 
pensates for the overlap and underlap of adjacent scan 
line sweeps that result from the orbital motion of the 


Landsat spacecraft. The scan line corrector ensures 
that individual TM scan lines are oriented perpen- 
dicular to the spacec'aft ’s ground track. 


The three tungsten lamps operate at approximately 
the same intensity, but different neutral density filters 
in front of each lamp cause a difference in the inten- 
sity of light from each lamp for calibration. All possi- 
ble combinations of two lamp states (on and off) for 
three lamps provide eight different light levels for 
detector calibration. These eign^ levels include a dark 
level consisting of all three lamps in the "off” state. 
The light level remains constant for approximately 20 
scan cycles (i.e., 20 successive forward and reverse 
sweeps of the scan mirror). After each 20 cycle inter- 
val, the state of one lamp is changed (i.e., turned on 
or ofO to alter the light level. Since 16 detectors are 
used for each band, a total of 5 1 20 scan lines of image 
data are generated during the time interval required 
to view all eight light levels (i.e., corresponding 'o 
eight calibration levels times 40 scan mirror sweeps 
per level times 16 lines of data per mirror per sweep). 
For the thermal band, the temperature of the 
blackbody and the temperature of the shutter’s dc 
restore surface are used for a two-level calibration of 
the four thermal detectors. 

Figure II-5 provides an overview of the Thematic 
Mapper design. The TM is horizontally mounted in 
the satellite with the sun shade pointing toward Earth. 
The scan mirror is located directly above the sun 
shade aperture along with the mirror's drive 
mechanism, control electronics, and scan angle 
monitor. The primary mirror is secured approxi- 
mately in the center of the telescope and is preceded 
by the secondary mirror and optical baffles. The scan 
line corrector, the internal calibrator, and the primary 
focal plane are directly behind the primary mirror. 
The angular displacement sensors are located within 
a small assembly mounted to the exterior of the TM 
telescope (the assembly is hidden in Figure II-5). The 
back end of the instrument contains the relay optics 
and the cooled focal plane within the radiative cooler. 
The electronics package sits in a wedge-shaped box 


above the telescope and contains the multiplexer, 
power supplies, signal amplifiers, and electronic 
filters. 

The complex design of the TM affords advanced 
sensor attributes in comparison to the MSS. For in- 
stance, the TM acquires data in seven spectral bands. 
Four of these bands are located in portions of the 
spectrum not sensed by the MSS. The location and 
width of the seven bands were carefully chosen for 
sensitivity to certain natural phenomena and to 
minimize the attenuation of surface energy by at- 
mospheric water. Table 11-2 lists the spectra! bands 
along with the rationale for selecting each band. 

Another advanced feature of the TM is its spatial 
resolution. The TM instantaneous-field-of-view (IFOV) 
is 30 m-by-30 m for the six reflective spectral bands 
(bands TMl to TM5 and TM7) and 120 m-by-120 m 
for the thermal band (band TM6). The MSS provides 
a resolution of approximately 80 m-by-80 m in .u* 
bands, Table II-3 lists the optical characteristics of the 
TM which govern its surface spatial resolution. 
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Tab!® 11-2 

Thematic Mapper Spectnil Bands (Salomonson et al., 1980) 


Band 

Designation 

Specified 

Band 

Width (/rlD) 

Specified 
Radiometric 
Sensitivity (NEAg)* 

Selection Rationale 

TM1 

0,45“0.52 

0.8% 

Sensitive to chlorophyll and cartonold con- 
centrations for soil/vegetation differentia- 
tion, deciduous/ coniferous differentiation. 
Coastal water mapping. 

m 2 

0,52-0.60 

0.5% 

Sensitive to green reflectance by healthy 
vegetation. 

TM3 

0.63-0.69 

0,5% 

Sensitive to chlorophyll absorption for plant 
species differentiation. 

TM4 

0.76-0.90 

0.5% 

Sensitive to near infrared reflectance of 
healthy vegetation for biomass suiveys. 

TM5 

1,55-1.75 

1.0% 

Sensitive to vegetation moisture and 
stiow/cloud reflectance dlffe'-wcpi). 

TM6 

10.4—12.5 

0.5KNETD** 

Thermal mapping. 

TM7 

2,08-2.35 

2.4% 

Sensitive to vegetation moisture and to 
hydroxyl Ions In minerals for geological 
mapping. 


*NEAe: noise equivalent change in surface reflectance 
•*NETD: noi^e equivalent temperature difference 


Table 11-3 

Themsitic Mapper Optical Characteristics 
(Engel, 1980} 


40.6 cm. Aperture 
f/6 at prime Focus 

42:5 urad Angular IFOV, Kands TM1 through TM4 
f/3 at Relay Focus 

43:8 /zrad Angular IFOV, Bands TM5 and TM7 
170 /trad Angular IFOV, Band TM6 


The TM also offers improved radiometric sensitivity 
over the MSS, even though the spectral bandwidth 
and the IFOV of the TM detectors is significantly less 
than that of the MSS. Sensitivity specifications in terms 
of noise equivalent change in surface reflectance 
(NEAg) are presented for the MSS and TM in Tables 
Il-l and II-2, respectively, In conjunction with the im- 
proved sensitivity, TM data are quantized to eight bits 
(256 digital counts) while MSS data are quantized to 
six bits (64 digital counts). This effectively cor- 
responds to a four-fold increase in the gray scale be- 
ing used to measure the intensity of Earth radiation 
in each discrete spectral band. 


In summary, the TM’s improved sensor capabilities 
relative to the MSS are a result of the advanced design 
of the TM. Design factors which significantly affect 
the TM’s measurement capabilities include; the detec- 
tors’ angular IFOV which is a function of detector size 
and telescope focal length; the optics of the scan mir- 
ror and telescope which minimize diffraction and blur 
and thus enhance radiometric sensitivity and spatial 
resolution; the use of the bidirectional active scan ap- 
proach which increases detector dwell time and thus 
enhances radiometric sensitivity; the internal cali- 
brator which permits eight level, in-flight radiometric 
calibration; the relay optics and cooled focal plane 
which enable sensing of energy within middle in- 
frared and thermal infrared spectral bands; and the 
scan angle monitor and the scan line corrector which 
facilitate the precise geographic registration of the fine 
resolution image data, The rationale for creating this 
advanced design is to generate data which can be bet- 
ter exploited for a wide variety of Earth observations. 
For future refeience, the TM aboard Landsat 4 is called 
the protoflight (TM/PF) model while the TM aboard 
Landsat 5 is referred to as the flight (TM/F) model. The 
designs of the two TM’s are identical. Their operating 
characteristics differ in some minor respects as noted 
below. 
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tANDSAT ORBIT AND COVERAGE 

Landsat 4 and 5 travel In repetitive, circular, sun* 
synchronous, near-polar orbits at a nominal altitude 
of 705 km. The orbits have an Inclination angle of 
98,2 degrees with respect to the Earth's equator, The 
satellites cross the equator at approximately 9:45 a.m. 
local solar time during the descending (north-to-south) 
portion of each orbit. The TM can acquire thermal 
data during the satellites’ cscendlng (south-to-north) 
portion of an orbit during which the satellites cross 
the equator at 9:45 p.m. local solar time, Each orbit 
takes almost 99 minutes, and each spacecraft com- 
pletes just over 14.5 orbits a day. Due to Earth rota- 
tion, the distance between ground tracks for two con- 
secutive orbits is 2752 km at the equator. Figure II-6 
Illustrates the orbital characteristics of Landsat 4 and 
Landsat 5. 

The orbital configuration described above creates 
a 16 day repeat cycle for each satellite. In other words, 
a specific satellite will travel along a particular ground 
track once every 16 days. The orbits of the two 
satellites are eight days out of phase. In other words, 
Landsat 4 and Landsat 5 alternate orbits over a par- 
ticular ground track every eight days. For either one 
of these satellites, the time interval between adjacent 
tracks Is seven days (Figure II-7). This arrangement 
contrasts with the 18-day orbital cycles of the first 
three Landsat satellites which resulted in a one day 
Interval between orbits over adjacent tracks. 

The altitude, optics, and 15.4° scan angle of the TM 
enable the instrument to view a 1 85 km swath across 
the ground track of each satellite. At the equator, ad- 
jacent swaths overlap at their east and west margins 
by 7.6 percent. The overlap increases as the ground 



Figure II-6 The sun-synchronous orbits of Landsat 4 
and 5. 


tracks move away from i;hc equator toward either 
pole. At 40 degrees north or south latitude, for ex- 
ample, the overlap is 18.5 percent. 

The Landsat 4 and 5 pass over the same ground 
track after completing 235 orbits of the Earth. This 
orbital repeat cycle is Incompatible with the 251 orbit 
Worldwide Reference System (WRS) used to Index 
MSS data from the previous Landsat satellites, A new 
system of path and row indices has been designed for 
both MSS and TM data from Landsat 4 and 5- A path 
index number is assigned to the ground track of each 
orbit. The new system consists of 233 orbital paths 
numbered east-to-west with the center point of path 
001 crossing the equator at 64.95 degrees west 
longitude, The data for each scene consist of the data 
acquired during a 25.87-second increment of the or- 
bital motion of the satellite. This increment cor- 
responds to approximately 5700 scan lines of Image 
data. Both the MSS data and the TM data for each 
scene cover a surface area approximately 170 km 
along- track by 185 km across-track, A row index is 
assigned to each scene center, Row one of each path 
starts at 80°51 ' north latitude and the numbering in- 
creases southward through the equator (row 60) to 
81 °51 ' south latitude (row 122), The satellites then 
travel northward for the ascending (nighttime) por- 
tion of each orbit, and row numbers increase from 
south-to-north through the equator (row 184) to 
81 °51 ' north latitude (row 246), Thus, row Indices 
23 through 246 refer to scenes of nighttime TM ther- 
mal data which are not acquired on a routine basis. 
Maps delineating paths and rows can be obtained 
from the National Oceanic and Atmospheric Ad- 
ministration (NOAA) by writing to the following ad- 
dress: NOAA Landsat Customer Services, Mundt 
Federal Building, Sioux Falls, South Dakota 57198, 


TM DATA ACQUISITION 

As the TM scan mirror sweeps back-and-forth sevta 
times per second, incoming Earth energy is focused 
onto a total of 100 photosensitive detectors. The 
detectors translate the incoming energy into low-level 
electrical signals which are amplified, quantized to 
eight-bit digital words, and then multiplexed into an 
84,9 niegablts-per-second (Mbps) data stream. The 
digital image data are then transmitted to the ground 
along with additional data concerning TM sensor 
operation and spacecraft conditions. 

The first three Landsat satellites relied upon on- 
board tape recorders to store data until the satellite 
passed within range of a ground receiving station. 
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Figure II-7 The Landsat 4 and 5 orbital repeat cycles. 


These recorders were often the first component to 
fail. Landsat 4 and 5 use a new communications 
sytem, called the Tracking and Data Relay Satellite 
System (TDRSS), which eliminates the need for on* 
board recorder,*. TDRS.S will eventually consist of two 
satellites (i^lus an in-orblt spare satellite) In geosyn- 
chronous orbits and a centralized ground receiving 
station located at White Sands, New Mexico. The first 
TDRSS satellite (TDRS-Rast) was launched on April 4, 
1983 aboard NASA’s space shuttle Challenger (mission 
STS-6). After overcoming problems caused by a 
malfunction of its booster rocket, the satellite reached 
Its permanent station on October 17, 1983. TDRS-East 
is currently In geosynchronous orbit over the equator 
at 41 ° west longitude at an altitude of 35,890 km. The 
second TDRSS satellite is scheduled for launch In 
February, 1985. When launched, the second TDRSS 
satellite (TDRS-West) will be located 35,890 km over 
the equator at 171 ° west longitude. 

The performance of TDRS-East Is still being tested, 
but the satellite is currently (as of July 1, 1984) relay- 
ing Landsat 5 TM and MSS data to White Sands on 
a regular basis. The relay is limited to the data from 
two or three Landsat 5 ground tracks per day during 
the current system checkout period. The daily volume 
of TM data transmitted by the system will increase 
when TDRS-East Is considered fully operational. 

With only one TDRSS satellite in orbit, the data 
relay is also limited to TM and MSS data acquired bet- 
ween approximately 160° west longitude and 60° 
east longitude (Figure II-8). This excludes data ac- 
quired over Australia, India, and Asia from TDRS-East 
relay. When the second TDRSS satellite, TDRS-West, 
is placed in orbit in its final position, only data ac- 
quired over India and the Indian Ocean (between 67® 
east and 82° east longitude and between 50° north 
and 50° south latitude) will be excluded from TDRSS 
relay (Figure II-8). It is currently anticipated that the 
complete TDRSS system will become fully operational 
during the summer of 1985. 

To handle the high data rates from Landsat 4 and 
5, the TDRSS employs a Ku-band (approximately 15 
Gigahertz) frequency for data communications. This 
frequency is transmitted from the Landsat TDRSS 
antennas (Fig. II-l) and then relayed through the 
TDRSS satellites to the ground receiving station at 
White Sands, New Mexico. The data received by the 
White Sands facility are demodulated, separated, and 
recorded on wlde-band digital data tapes. These com- 
pacted raw data are read from the tapes and then 
transmitted, by way of a domestic communications 
satellite (DOMSAT), to a data processing facility at 
NASA’s Goddard Space Flight Center (GSFC) in 
Greenbelt, Maryland. This scheme enables the rscep- 



Figure II-8 Regions where TM data cannot be relayed 
via one or both of the TDRSS spacecraft. 


tion of data at the GSFC facility within eight hours 
of data generation by the sensor. After geometric and 
radiometric corrections are performed at GSFC 
(described below), the processed TM data, on film and 
on computer compatible tapes, are sent by air ship- 
ment to the Earth Resources Observation System 
(EROS) Data Center in Sioux Falls, South Dakota, for 
distribution to the public, The EROS Data Center is 
operated for NOAA by the United States Geological 
Survey. Figure II-9 traces the critical pathway for TM 
data products. 

The TDRS data link is augmented by a capability 
to transmit TM data directly on an X-band (8.025 to 
8,4 Gigahertz) frequency. Landsat 4 and 5 both carry 
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Figure II-9 End-to-end flow of TM image data from 

X-band transmitters for this purpose. The Landsat 4 
X-band transmitter failed on February 15, 1983 after 
transmitting over 6000 TM scenes. Most of the Land- 
sat 4 TM data were received by the Transportable 
Ground Station (TGS) located at GSFC. Landsat 5 TM 
data are currently being received by the TGS and by 
six additional stations at the following locations: 
Cuiaba, Brazil; Prince Albert, Canada; Hyderabad, In- 
dia; Fucino, Italy; Tokyo, Japan; and Kiruna, Sweden. 
Receiving stations are planned for several more loca- 
tions around the world including Australia and China. 

Along with the TM image data of the surface scene, 
Landsat 4 and 5 transmit ancillary spacecraft and sen- 
sor data required for the radiometric and geometric 
correction of the image data. These additional data 
include: radiometric calibration data consisting of the 
quantized responses of the photosensitive detectors 
to energy from the internal calibration lamps de- 
scribed above; mirror scan correction data, from the 
scan angle monitor, describing the angular position 
of the scan mirror as a function of time; and payload 
correction data describing the attitude and the 
ephemeris of the sensor. Payload correction data con- 


the Landsat 4 and 5 spacecraft to ground-based users. 

sist of data from the angular displacement sensor, data 
from the satellite’s attitude control system. Earth- 
centered ephemeris, and housekeeping information 
such as timing codes and TM internal temperatures. 
The ephemeris is computed by on-board computers 
which receive information on the satellite’s position 
from ground stations of the NASA communications 
network (NASCOM). Landsat 4 and 5 are also capable 
of receiving information from the Global Positioning 
Systems (GPS) satellites, but this system is not used 
at the present time because the GPS satellite constella- 
tion has not yet been fully deployed. Multiplexers 
aboard Landsat 4 and 5 place all of the payload cor- 
rection data described above into a standard format, 
These ancillary data are then Integrated with the TM 
image data (i.e., detector outputs). The integrated data 
stream is commonly referred to as wideband data. 
Wideband data are transmitted to the ground via the 
X-band or Ku-band antennas on the Landsat satellites. 
Payload correction data are also transmitted separately 
at a rate of 32 kilobits per second by S-band (2206 
to 2300 Megahertz) telemetry to NASCOM ground sta- 
tions. NASCOM then forwards the payload correction 
data to GSFC. 
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TM DATA PROCESSING 

All of the ground stations receiving TM data have 
or plan to have associated data processing facilities. 
The operators of these stations have been organized 
into the Landsat Ground Station Operators Working 
Group (LGSOWG), The LGSOWG has agreed to some 
digital image data format standards, but the data pro- 
cessing procedures and available TM data products 
will likely vary between facilities. The discussion 
presented here describes procedures performed at the 
Landsat data processing facility located at NASA ’s 
Goddard Space Flight Center (GSFC) in Greenbelt, 
Maryland, This facility, called the Landsat Ground 
Segment, processes all the TM data relayed to the 
ground through TDRSS. The facility also processes 
data received by the Transportable Ground Station 
(TGS) located at GSFC and by the Prince Albert, 
Canada, station on the X-band frequency. 

The GSFC facilities consist of the Mission Manage- 
ment Facility (MMF), the Control and Simulation 
Facility (CSF), and the Image Generation Facility (IGF), 


The MMF receives and processes user requests for 
data covering specific ground locations and is respon- 
sible for data archiving. Data requests are then sent 
to the CSF which commands, controls, and monitors 
the Landsat satellites and also schedules TDRSS, DOM- 
SAT, and NASCOM support. Processing of MSS and 
TM image data is performed by separate systems In 
the IGF (see Figure IMO). 

The IGF sysifim for TM data is called the Thematic 
Mapper Image Processing System (TIPS). TIPS com- 
menced operations in July, 1983. Before then, TM 
data were processed for research purposes by an im- 
provised computer system popularly known as the 
"Scrounge System.” 

The Scrounge System was assembled from available 
pieces of computer hardware prior to the launch of 
Landsat 4. The system consisted of the Applications 
Development Data System (ADDS), the Landsat 
Assessment System (LAS), and elements of the Mission 
Management Facility (MMF). ADDS is primarily a test 
system designed to evaluate the use of array pro- 
cessors for the efficient processing of digital image 
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data. The central processing unit of ADDS Is a 
minicomputer which serves as a host for two array 
processors. The original intent of the LAS was to ex> 
pedlte the scientific utilization of TM data by pro- 
viding investigators with an advanced Image analysis 
facility. The system supports interactive Image 
analysis stations which Include cathode ray tubes for 
the video display of Image data. The LAS central pro- 
cessing unit is a minicomputer supported by an array 
processor (Fischel, 1983). Another two computers 
within the MMF are used to process payload correc- 
tion data for both the Scrounge System and TIPS. 
ADDS and LAS were pressed Into service for TM data 
processing until TIPS began operations. 

TIPS consists of two Identical systems operating in 
parallel, The central processing unit of each system 
is a minicomputer interfaced to an advanced, high- 
speed array processor. The advanced array processors 
enable TIPS to process an entire TM scene In under 
15 minutes (Fischel, 1983). 

The main functions of both TIPS and the Scrounge 
System are the radiometric and geometric correction 
of TM image data. The Scrounge System correction 
algorithms emulated the procedures implemented by 
TIPS, and experience gained during Scrounge System 
processing was incorporated into the final develop- 
ment of the TIPS. The Scrounge System, however, 
operated at a much slower throughput rate than the 
specified rate for TIPS. The Scrounge System process- 
ed data at an approximate rate of one TM scene per 
day. TIPS currently processes data at a rate of 12 
scenes per day to fully processed film products, and 
two scenes per day of completel)? processed data on 
computer compatible tapes (CCT’s). This rate will In- 
crease to 50 completely processed scenes per day on 
CCT’s, 

Data processing begins with the reception of a 
scene of TM wideband data at the Data Receive 
Record Transmit System (DRRTS) within the Image 
Generation Facility and the reception of the 
associated, but separated, payload correction data 
(PCD) by the Control and Simulation Facility (Figure 
II-IO). DRRTS sends the mirror scan correction data 
to the Mission Management Facility (MMF) for 
geometric processing (described below) and passes 
the rest of the wideband data to TIPS, 

To facilitate the discussion of TM data processing, 
the radiometric correction procedures will be ex- 
plained first followed by a discussion of geometric 
correction procedures. The actual sequence of data 
processing steps employed by TIPS will then be 
presented. 


Radiometric Correction Procedures 


The two primary objectives of radiometric correc- 
tion are to remove detector-to-detector striping and 
to permit conversion of TM digital counts to absolute 
radiometric units. The conversion of digital counts 
to absolute radiometric units enables Investigators to 
measure spatial and temporal variations in scene ra- 
diance in a consistent, quantitative fashion. Recall that 
the TM employs 16 detectors for each reflective spec- 
tral band (bands TMl to TM5 and TM7) and four 
detectors for the thermal band (TM6). The amplified 
and quantized signals from an individual detector are 
referred to as a data channel. Each datum from a chan- 
nel consists of an integer value between zero and 255, 
inclusive, and is called a digital number (DN). For each 
channel, a response function mathematically relates 
digital numbers to absolute spectral radiance (express- 
ed in units of power per unit area per unit solid angle 
per unit wavelength) incident on the TM entrance 
aperture. 

The response function for a data channel may vary 
in orbit due to variations in the performance of in- 
dividual detectors and other system components. The 
internal calibration lamps and thermal references 
(discussed above) at 2 used to characterize the in-orbit 
response functions of all TM data channels. The con- 
tinued cycling of the lamp states through eight calibra- 
tion light levels permits the monitoring of channel 
response functions at regular intervals. Without 
radiometric corrections, variations between the 
response functions for the 16 channels within a reflec- 
tive spectral band (or four channels within the ther- 
mal band) would be manifested as horizontal strip- 
ing in pictorial Imagery. The correction procedure for 
TM data minimizes striping by normalizing the digital 
numbers to a common response function for each 
spectral band. The corrected digital numbers can then 
easily be converted to absolute radiometric units by 
Inverting the appropriate common response function 
for a particular band. 

Several assumptions are implicit in the radiometric 
correction procedure. First, the TM detectors and data 
system were designed to generate a linear response 
to incident spectral radiance. Assumptions of linear 
response appear valid on the basis of prelaunch 
calibration of channel respon.se functions and post- 
launch analyses of channel response to the internal 
calibration lamps. Linearity simplifies the mathemati- 
cal expression of channel response functions. A func- 
tion can be defined by the slope and intercept of the 
line relating ditigal number to spectral radiance. The 
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Figure II- 11 Radiometric response function of an 
individual TM channel. 


slope and intercept of a linear response function are 
called gain and offset, respectively (Figure II-l 1), 
Equivalently, a linear response function can be de- 
fined by the spectral radiance corresponding to a 
digital number response of zero (RMIN) and by the 
minimum radiance required to generate the maximum 
digital response of 255 (RMAX; i.e., the saturation 
radlr' ce). The range from RMIN to RMAX is the dy- 
namic range of a channel. 

A second assumption involves the stability of the 
internal calibration lamps, The spectral radiance In- 
cident on each TM detector at each of the eight 
calibration levels (or two thermal reference levels) was 
measured during pre-launch calibration. The calibra- 
tion light levels Incident on each detector are assumed 
to have remained constant during launch and in orbit. 
Changes in observed channel responses to calibration 
light levels are attributed to changes in the perior- 
mance of TM detectors and/or data system compo- 
nents, not to changes in calibration lamp output. This 
assumption is difficult to test once the sensor has been 
placed in orbit. 

The actual radiometric correction procedure for TM 
data consists of several steps. These steps are sum- 
marized here and described in greater detail below. 
In the first step, output values for malfunctioning 
detectors are replaced by the output values of adja- 
cent detectors in the same band. This correction has 
only been needed for two channels from the Landsat 
4 TM/PF sensor. In the second step, channel 
responses to the on-board calibration standards are 
used to characterize the in-orbit response functions 
of the individual channels, These in-orbit response 
functions are up-dated at Intervals of 360 scan cycles 
(i.e,, every 11520 scan lines). In the third step, a 
histogram equalization procedure Is used to refine the 


channel response functions. The purpose of histo- 
gram equalization Is to prevent residual radiometric 
striping caused by Imprecision in the In-orblt calibra- 
tion of the channel response functions. Once the in- 
orbit response functions of individual channels have 
been calibrated and refined, the raw digital numbers 
can be normalized to a common response function 
for each spectral band. In the fourth and final pro- 
cessing step, the normalization is accomplished by 
defining and applying radiometric look-up tables 
(RLUT's) which map raw digital numbers to corrected 
digital numbers, For a specific spectral band, the full 
range of corrected digital numbers represents the 
same dynamic range of spectral radiance regardless 
of the Individual channel that generated the raw data. 
This common dynamic range remains constant from 
scene-to-scene. The common dynamic range for each 
band was specified on the basis of a pre-launch 
calibration of the individual channels within each 
band. The collective result of this multiple-step 
radiometric correction procedure is the production 
of TM imagery which exhibits minimal radiometric 
striping. In addition, this procedure generates 
radlometrically calibrated TM image data which can 
be used to measure differences in scene radiance In 
a quantitative fashion both within and between in- 
dividual TM scenes. 

In the radiometric processing of TM/PF data from 
Landsat 4, data from two channels are replaced; the 
third detector of band TM5 has no radiometric 
response and the fourth detector of band TM2 has a 
poor spatial resolution. The digital numbers from each 
malfunctioning detector are replaced by the raw 
digital numbers from an adjacent detector In the band. 
The replacement data are then treated in the same 
manner as the original data from the functioning 
detectors during subsequent processing. All of the 
detectors in the Landsat 5 TM/F are operating and data 
replacement is unnecessary. 

The next step determines the in-orblt linear 
response function for each data channel. The digital 
responses of individual channels to the internal lamps 
are fit to the eight calibrated (and presumably un- 
changing) levels of Incident spectral radiance from the 
Internal calibration lamps by simple linear regression 
(i.e., a least-squares fit) for the six reflective spectral 
l5ands (Figure II-ll), The two internal thermal 
references are monitored by on-board heat sensors. 
The response function of each of the four thermal 
channels is determined by a linear regression between 
the observed digital responses of a thermal channel 
and the measured temperatures of the two Internal 
thermal standards. The slopes and Intercepts of these 
regression lines described above correspond to the 
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gains and offsets, respectively, of the measured chan- 
nel response functions. 

The possibility of Imprecision in the estimation of 
Individual channel response functions was recognized 
during the design of the radiometric correction pro- 
cedure. A histogram equalization process was In- 
cluded at this stage In the correction procedure to 
refine the estimates of channel response functions. 
The purpose of histogram equalization Is to minimize 
residual llne-to-llne striping In final TM Image prod- 
ucts which could occur given imprecise estimates of 
channel gain and offset. This procedure has been 
routinely used In the past for relative radiometric cor- 
rections of MSS data, 

The equalization procedure begins with the pro- 
duction of a digital number frequency histogram for 
each channel. The histogram is constructed by sam- 
pling every fourth digital number generated by a 
channel during an Interval of 360 scan cycles. This 
sampling protocol Is based upon the assumption that 
each channel within a band views the same distribu- 


tion of Earth radiation during a 360 scan cycle Inter- 
val, The mean and standard deviation of the digital 
numbers sampled for each histogram are subsequently 
computed. The mean and standard deviation for each 
histogram are then converted to absolute units of 
spectral radiance using the original estimates of chan- 
nel gain and offset and the inverse of the radiometric 
response function (see Figure II-l 1). 

At this point in the process, if all estimates of gains 
and biases were precise and all distributions of ob- 
served radiance were equal, then the histograms for 
all channels within a particular band would have equal 
means and equal standard deviations In terms of ab- 
solute spectral radiance units (e.g., In terms of 
milliwatts per square centimeter per steradian per 
micrometer), These equalities rarely occur due to Im- 
precision in gain and offset estimates or Inequality be- 
tween observed radiance distributions. To generate 
equal histograms In terms of absolute spectral radiance 
units, the original estimates of channel gain and off- 
set are refined in the following manner:* 


(1) G^b = Gc.b ^ 

Sb 

(2) Bc,b “ ^c,b + Gc,b fMc,b “ ~ ^ 

\ Sc.b / 


where 


Gc,b and B(-^b 


^c,b 

Mb 


Sc.b 

Sb 


Gc,b and Bc,b 


are the original estimates of gain and offset, respectively, for channel c and 
spectral band b. These estimates are derived from In-orbit observations of the 
TM internal calibration standards. Gc,b is in the following units: 
DN/mW»cm'2»sr'»^m-‘, Bc,b is in units of DN. 

is the mean associated with the histogram for channel c, band b, in units of ab- 
solute spectral radiance; mW»cm'^»sr‘»/im, 

is the average of the l6 histogram means (Mc.b’s) for all channels In a spectral 
band (four means for the thermal band) In units of absolute spectral radiance: 
mW‘cm‘2«sr->*^m''. 

is the standard deviation associated with the histogram of channel c, spectral 
band b, in units of absolute spectral radiance: mW»cm-2»sr-'»ftm-‘. 
is the average of the 16 histogram standard deviations (Scab’s) for all channels in 
a spectral band (four standard deviations for the thermal band) in units of ab- 
solute spectral radiance: mW»cm'2»sr-‘»/im-'. 

are the refined estimates of Gc,b and Bc,b» respectively, and are in the same units 
as Gc,b and Bc,b, respectively. 


'Note: Data from the TM thermal band are calibrated and corrected in terms of absolute temperature. When considering the 
correction of thermal data In the discussion above, units of absolute spectral radiance (mW«cm'2*sr**/im'') should be con- 
sidered units of absolute temperature (degrees Kelvin). 
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If the refined estimates of Individual chant.el gain 
and bias were used to convert the histogram means 
and standard deviations from digital numbers to ab- 
solute radiometric units, then the histograms for all 
channels within a spectral band would have equal 
means and equal standard deviations in terms of ab- 
solute spectral radiance. The effect of this equaliza- 
tion process on the absolute radiometric accuracy of 
corrected TM data is unclear (Barker, 1984), but very 
little channel-to-channel striping is apparent in cor- 
rected TM imagery (Salomonson, 1984). 

The preceding correction steps account for 
variability In the orbital performance of individual TM 
channels and they serve to minimize image striping. 
The last procedural step transforms raw digital 
numbers (DN*s) to corrected DN’s using the refined 
estimates of channel gains and biases. The Intention 
of the transformation is as follows: for a given chan- 
nel, the corrected DN’s can be related to absolute 
units of spectral radiance by a linear response func- 
tion which remains constant from channel-to-channel 
and between scenes. In other words, the full range 
of corrected DN’s for a specific band represents a con- 
stant dynamic range of spectral radiance no matter 
which channel generated the raw data or when the 
raw data were acquired, This transformation facilitates 
the quantitative comparisons of surface radiance both 
within scenes and between scenes. 

A constant dynamic range (i.e., radiometric 
response function) was specified for each spectral 
band on the basis of prelaunch calibration of the in- 
dividual channel response functions. The smallest 
dynamic range observed for a single channel in each 
band was chosen as the standard to which all wlthln- 
band data would be normalized. This choice prevents 


data compression. In other words, this choice allows 
a one-to-one mapping between corrected DN's and 
raw DN’s except at saturation (i.e., except at a cor- 
rected value of 255). The standard dynamic ranges 
(I.e,, RMIN and RMAX values) selected for the spec- 
tral bands are shown In Table II-4. The values selected 
for Scrounge System data processing (Table II-4a) 
were based on the prelaunch calibration of the Land- 
sat 4 TM/PF sensor. The dynamic range values chosen 
for TIPS were slightly altered on the basis of the pre- 
launch calibration of the Landsat 5 TM/PF sensor 
(Table ll-4b). TIPS uses the values given in Table II-4b 
for the radiometric correction of both TM/PF and 
TM/F data. As shown in Table 11-4, the standard 
dynamic range values define the common linear 
radiometric response functions to which the data 
from all channels within a band are normalized. 


The common gains and offsets for a spectral bands 
(see Table II-4) are used in conjunction with the re- 
fined estimates of Individual channel gains and off- 
sets (see Figure II-l 1 and equation 1 and 2) to define 
the transformation of raw DN’s to corrected DN’s, 
The transformation is based on the following 
equation: 


(3) Dcor.b “ 7n~ ^raw.c.b + Bb " "pr, 

tJc.b "-^c.b 


where, 

Draw.c.b is a raw digital number generated by 
channel c of spectral band b; 

Gc,b is the refined estimate of the gain of the 
radiometric response function for chan- 
nel c, band b (see Figure II-l 1 and equa- 
tion 1); 


Table IS-4 

Dynamic Ranges of Radlometrically Corrected Thematic Mapper Data 


Table ll-4a 

Dynamic Ranges of Thematic Mapper Data Processed by the 
Scrounge System Prior to August, 1983 (Barker, 1984). 


Spectral Bands 


TM1 TM2 TM3 TM4 TM5 TM7 TM6* 

(thermal) 


RMIN (mW*cm'2*sr‘i*/«m‘i) 
RMAX (mW*cm‘2*sr’*/im'’) 


-0.15 

15.84 


-0.28 

30.82 


- 0.12 

23.46 


-0,15 

22,43 


-0.04 

3.24 


- 0.02 

1.70 


260°K 
320 “K 
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Table IMb 

Dynamic Ranges of Thematic Mapper Data Processed by TIPS 
After January 15, 1984 (Barker, 1984). 


Spectral Bands 



TM1 

TM2 

TM3 

TM4 

TM5 

TM7 

TM6* 

(thermal) 

RMIN (mW*cm’**sr‘>*/jm’1) 

•0,15 

-0.28 

-0.12 

-0.15 

-0,04 

•0.02 

200“K 

RMAX (mW*cm‘2*sr’*^m'>) 

15.21 

29.68 

20.43 

20.62 

2.72 

1.44 

304“K 


RMIN is the spectral radiance corresponding to a detector response of zero digital counts 

RMAX Is the minimum radiance required to saturate detector response (i,e., produce a response of 255 digital counts) 


Common Radiometric Response Function for All Corrected Data In a Spectral Band 
Dcor.b ” Lb + E)b 
where 

Dcor.b Is a corrected digital number from spectral band b; 

Lb Is spectral radiance (mW*cm'2»sr’i»/*m'i) within band b; 

Gb Is the common gain for band b In the following units: 

DN/mW»cm*2»sr‘i 
Gb - 255/(RMAX-RMIN) 

Bb Is the common offset for band b in units of DN; 

Bb - "(255 RMIN)/(RMAX-RMIN) 

*The thermal band, TM6, is calibrated with respect to black body temperature, When considering the corrections of TM ther- 

mal data above, units of spectral radiance (mWcm'^sr‘'/<m"') should be considered units of absolute temperature (degrees Kelvin). 


Is the refined estimate of the offset of the 
radiometric response function for chan- 
nel c, band b (see Figure IM 1 and equa- 
tion 2); 

Gb and Bb are the gain and offset, respectively, 
of the common radiometric response 
function for band b (Table II-4); 

Dcor.b Is a corrected digital number for band b. 

To expedite computer processing, equation 3 is 
used to derive radiometric look up tables (RLUT’s) 
which map raw digital counts to ?he corrected values, 
One RLUT is generated for each of the 100 TM chan- 


nels, and the RLUT’s are used to actually conduct the 
correction. Corrected DN values (Dcor.b) appear in all 
processed TM digital image data products (l.e,, on the 
computer compatible magnetic tapes containing the 
A-data and P-data described below). 

A corrected digital number can be converted to ab- 
solute units of spectral radiance by inverting the com- 
mon radiometric response function (Table II-4) for the 
appropriate spectral band. The inverted function can 
be expressed as follows: 


(4) Lb 


(RMAX-RMIN) 

255 


Dcor.b 


+ RMIN 
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where 




Dcor.b is a corrected digital number for spectral 
band b; and 

Lb, UMIN, and RMAX are defined In Table 

The total wlthln-band radiance corresponding to 
a particular digital number, Dcor.b* be estimated 
by multiplying the appropriate value of spectral ra- 
diance, Lb, by the bandwidth of the spectral band 
under consideration. A more precise estimate' of 
wlthln-band radiance requires more precise 
knowledge of spectral response curves for the TM 
bands. Markham and Barker (1984) provide 
measurements of TM bandwldtlis and spectral 
response curves. 

Geometric Correction Procedures 

The purpose of geometric correction is to create 
a nearly conformal representation of the Earth sur- 
face consisting of TM multlspectral digital image data. 
This conformal representation Is generated by first 
computing correction data and then applying the cor- 
rection data to the resampling of TM image data. The 
correction data locate the position of each radiomet- 
rically corrected TM image data sample (i.e., each 
vligltal number) onto an output coordinate system 
which corresponds to a cartographic projection. A 
resampling procedure jissigns a digital number to each 
row/column coordinate (i.e., pixel) of the output 
system on the basis of Interpolation between 
neighboring input data samples. The resampled TM 
data are evenly spaced with respect to cartographic 
location which enables the generation of nearly con- 
formal pictorial imagery of the Earth's surface. 

At any one instant, the 100 TM detectors all view 
different areas on the Earth surface. The instantan- 
eous-field-of-view (IFOV) of an individual detector 
varies in size and location with satellite altitude, 
satellite ephemeris (i.e., location with respect to the 
center of the Earth), satellite and sensor attitude, and 
scan mirror motion. Registering these unevenly- 
spaced detector IFOV’s to an evenly-spaced, carto- 
graphlcally-based coordinate system accomplishes 
three main objectives. First, the resampled data are 
registered band-to-band. This means that an object 
depicted in the corrected data from a particular band 
can be located at the same coordinates within the cor- 
rected data from any other spectral band from the 
same TM scene. Second, corrected TM data are tem- 
porally registered. An object depicted within a par- 
ticular TM scene can be located at the same coordinate 


in any other TM scene acquired over the same geo- 
graphic area on a different date (I.e., a prcvlou.s or 
subsequent overpass). Third, the corrected TM data 
are geodetlcally registered. The cartographic location 
(e.g., latitude and longitude coordinates) of an object 
can be accurately determined from the object’s 
row/column coordinates in the corrected Image data. 
These band-to- band, temporal, and geodetic registra- 
tions generally render corrected TM data more useful 
than rawTM data for Investigators. Specifications for 
registration accuracy are given In Table II-5. 


Table II-5 

Geometric Registration Accuracy Specifications 
for Thematic Mapper Data (Beyer, 1984) 


Type of Accuracy Required for 

Registration 90% of the Data Samples 


Band-to-Band for Spec;ral Bands 

on the Sr me Focal Plane* 0,2 pixels 

Band-to-Band for Soectral Bands 

on Different Focal f lanes* 0.3 pixels 

Temporal 0.3 pixels 

Geodetic 0.5 pixels 

* The detectors for spectral bands TMl toTM4 areTocated 
on the primary focal plane. The detectors for spectral bands 
TM5, TM7, and TM6 (thermal band) are located on the cool- 
ed focal plane (see Figure 11-3), 

The output coordinate system for corrected TM 
data is an evenly-spaced grid where each grid cell, 
called a pixel, represents a 28.5 m-by-28.5 m area on 
the Earth’s surface. Each pixel is identified by Its 
row/column coordinate and each pixel is assigned 
seven digital numbers for the seven TM spectral 
bands. Most corrected TM scenes are registered to a 
coordinate system corresponding to the Space 
Oblique Mercator (SOM) cartographic projection. In- 
vestigators may make special requests for use of the 
Universal Transverse Mercator (UTM) or the Polar 
Stereographic projections, 

The geometric correction data define the inter-pixel 
location of each raw TM image data sample on the 
output coordinate system. The geometric correction 
data consist primarily of three types of data matrices. 
Two types of matrices, the Benchmark Matrices and 
the High Frequency Matrices, trace the projected TM 
optical axis across the Earth surface as a function of 
time, These matrices are derived from the payload 
correction data (PCD) and the mirror scan correction 
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data (MSCD) described cariicr. The Benchmark 
Matrices trace tlie TM optical axis as a function of 
satellite ephemeris and scan mirror position assum* 
Ing nominal operating conditions (l.c., perfect sensor 
and satellite attitude and linear scan mirror motions 
and scan line corrector motions). The High Frequency 
Matrices account for frequent attitude deviations and 
non-llnearltlcs In scan mirror motions and scan line 
corrector motions as measured by PCD and MSCD. 
A third matrix type, the Focal Plane Geometry 
Matrices, geometrically relates the individual detec- 
tor IFOV’s to the TM optical axis. This relationship 
remains constant with time. The Focal Plane 
Geometry Matrices were derived from prelaunch 
measurements of TM optics and detector locations on 
the focal planes and remain constant from scene-to* 
scene, The three types of matrices arc multiplied to 
form a master set of matrices which locate the pro- 
jected IFOV of each TM detector on a cartographic 
projection as a function of time, This final set of 
matrices Is collectively called systematic correction 
data (SCD) (Beyer, 1984). 

By use of the payload correction data, the deriva- 
tion of the systematic correction data (SCD) takes Into 
account the following factors: satellite ephemerisi sen- 
sor and satellite attitude; high frequency structural 
disturbances, commonly called Jitter, caused by 
operation of the instruments aboard th. atelllte; 
nonllnearitles In the scan mirror profile and opera- 
tion of the scan line corrector; the TM optics and the 
location of the detectors on the focal planes; and the 
desired cartographic representation of the three- 
dimensional surface viewed by the sensor, Errors, 
however, may remain in the spatial transformation 
defined by the SCD due to limits In the accuracy of 
the payload correction data (PCD) used In the com- 
putation of the SCD. To nilnlmL^e these errors. TIPS 
Is capable of using ground control points to adjust the 
SCD. The adjusted data are called geodetic control 
data (GCD), 

The adjustment of the SCD requires the building 
of a control point library, for each scene area (i.e., 
path/row location) viewed by the TM. Each “point" 
in a library actually consists of a 32 plxel-by-32 scan 
line data set extracted from a base line TM scene. The 
base line scene has been radiometrically corrected and 
geometrically corrected using SCD. 

The extracted data segments are called “chips", and 
each chip is located on a standard map (e.g., USGS 
topographic sheet). The latitude, longitude and eleva- 
tion relative to the standard geoid (l.c., the idealized 
surface of the Earth which coincides everywhere with 
mean sea level) is determined for an easily Identified 
point within each chip from the map. The number 


of chips extracted from a particular scene depends on 
the number of stable geographic features (e.g., road 
Intersections) which can be easily Identified on TM 
imagery and on the map. Control point libraries arc 
not yet available for all TM scenc.s. If a library Is 
unavailable for a particular scene, SCD are used for 
the TM data resampling. 

Given a control point library with Its cartographic 
information, SCD are used to define the expected 
locations of the ground contro' points It an uncor- 
rected TM scene. A 128 pixel*by*256 scar line search 
neighborhood Is then extracted for each control point 
from the uncorrccted TM data. Each neighborhood 
Is centered on the expected control point location and 
the neighborhood Is sufficiently la: ge to contain the 
associated 32 pixel- by-32 scan line control point chip 
even for large errors lr> SCD. Each neighborhood is 
also radiometrically corrected using RLUT’s and 
geometrically corrected using SCD. The precise loca- 
tion of the control point within a neighborhood Is 
then determined by correlation between the base line 
chip and the search neighborhood in the spatial fre- 
quency domain. The correlation process Involves 
gradient enhancement of edge features to reduce tem- 
poral differences between the base line data and the 
neighborhood data (Beyer, 1984). The difference be- 
tween tne precise location of the control point In the 
neighborhood, as determined by correlation, and the 
expected location is considered mislocatlon error. 
The mislocatlon errors for all available ground con- 
trol pol ts are used to adjust the SCD to create 
geodetl control data (GCD), The use of SCD to define 
TM data resampling produces acceptable band-to- 
band registration accuracy, but GCD are generally re- 
quired for accurate temporal and geodetic 
registration. 

The GCD (or If unavailable, the SCD) define the 
Inter-pixel lOcafis-Jii of each uncorrected TM data sam- 
ple on the output coordinate system. For resampling, 
the GCD (or SCD) project the TM data samples onto 
the standard geoid and topographic effects are not 
considered (Beyer, 1984). Resampling then assigns 
digital numbers to the discrete, evenly-spaced pixels 
of the output system to create the digital format for 
the final, fully-processed TM data products. Two 
resampling procedures are available. Investigators 
may specially request the use of nearest-neighbor 
resampling. This procedure simply assigns the digital 
number value of the closest TM image data sample 
to the output pixel. 

Cubic convolution Is the standard procedure used 
to resample TM data, This procedure assigns the 
weighted average of nearby TM samples to each out- 
put pixel. A cubic spline function defines the weighing 
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factors wliicli decrease with Increasing distance from 
the output pixel. Cubic convolution is implemented 
In a two step process (Figure The data samples 
along each uncorrected scan line arc first resampled 
to assign digital counts to temporary pixels located 
at the Intersection of the output columns and the 
Input scan lines. The teniporary pixels are then 
resampled to assign digital counts to the equally- 
spaced output pixels. 




• IXtlHSlOH P Xtl6 Am GtMdATCO WITH lyVO UMPOnAIlT Pl*U0 ADOVt 
AND TWO TtMPUnAny PIXU3 OflOW 

• I lOCATlON PfltCIOION TO ll «l PIXtk IH Pf »AyPllNC COMPUTATION 

lUguro H -15 Scan line extension for geometric 
correction. 


Sequence of TM Data Processing 
Steps 


Figure 11-12 Resampling of Thematic Mapper image 
data for geometric correction (Beyer, J984). 


The resampling procedures also compensate for 
spatial gaps which occur between the last scan line 
generated during a scan mirror sweep and the first 
scan line of the subsequent sweep. The gaps are a 
function of scene latitude and fluctuations in 
spacecraft altitude and attitude, To fill in the gaps with 
data, pixels are generated for locations on synthetic 
scan lines which extend into the spatial gap and run 
parallel to the scan lines of the upper mirror sweep 
(Figure 11-13). The extension pixels are then resam- 
pled to derive digital numbers for output pixels fall- 
ing In the gaps of the input image data. This resam- 
pling procedure is rather complex and is described 
in more detail by Beyer (1984), 

Geometric correction of the TM thermal data re- 
quires special consideration. Recall that only four 
detectors are used to acquire thermal data, and one 
image sample of raw thermal data represents an area 
equivalent to four image samples from any of the 
reflective bands. To facilitate multlspectral analyses 
of TM data, the coarse resolution thermal data is also 
resampled to form a registered grid of 28.5m- 
by-28,5m pixels. Thus, all bands of geometrically cor- 
rect TM data contain the same number of pixels per 
unit area. 


The sequence of TM data processing steps Is Il- 
lustrated In Figure II-14, and the flow of TM data 
through the GSFC data processing facilities Is shown 
In Figure 11-10. To initiate the sequence, TIPS accepts 
two or three input data sets. One set is the TM wide- 
band data containing the image data of the surface 
scene and the Internal calibration data. The Data 
Receive Record Transmit System (DRRTS) provides 
these data to TIPS on a high density tape (HDT), 
Another data set is the systematic correction data 
(SOD) derived by the Mission Management Facility 
(Ni.rlF) and provided to TIPS on a CCT. 

TIPS may also accept a set of ground control points, 
A ground control point library is currently being built 
for TIPS at the rate of one scene per day. If points 
are available, the systematic correction (SCD) are ad- 
justed to generate geodetic correction data (GCD) for 
geometric corrections. Otherwise, TIPS uses the SCD 
for corrections. Ground control points were not 
available during Scrounge System operations and the 
Scrounge System employed the SCD exclusively, 

TIPS begins processing by extracting the internal 
calibration data and by building histograms of digital 
detector output values for radiometric correctloas. 
Also early In the processing, TIPS locates control 
point search neighborhoods corresponding to the ex- 
pected locations of available ground cofurol points. 

Next, TIPS derives the channel gains and olTsets and 
generates the radiometric look-up tables (RJ.UT’s) for 
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Figure IF14 Overview of Thematic Mapper data pro- 
cessing flow at the GSFC facility (Lyon et al, 1985). 
Note that payload correction data (PCD) are received 

radiometric correction. If ground control points are 
available, TIPS correlates the control point 
neighborhoods to the known locations of the ground 
control points and adjusts the SCD to create GCD. 

Before applying the radiometric correct, 'ons, TIPS 
commutates the pixel order of the image data, aligns 
the data, and converts the data from a band- 
interlcaved-by-pixel to a band-interleaved-by-line 
(BIL) format. The commutation is made necessary b> 
the generation of image data during both the forward 
and reverse sweeps of the TM scanning mirror. Since 
16 scan lines of data are generated during each sweep 
of the scanning mirror, the directional order of pix- 
els along a line is reversed every 16 lines, The com- 
mutation alters the order of data acquired during each 
rever.se (east-to-west) mirror sweep to a west-to-east 
orientation. Once the data are oriented in the same 
direction, TIPS aligns the beginning of each scan line 
to form a contiguous image of the surface. 

The commutated and aligned image data are next 
radiometrically corrected using the RLUT’s, and the 
corrected data are stored on a high density magnetic 
tape (28 tracks, 33.3 kilobits per inch on each track). 
These data are considered partially processed having 


as a part of the wideband data and are also received 
separately via S-band telemetry. 


undergone radiometric correction but not geometric 
correction, TIPS also appends header, ancillary, an- 
notation, and trailer (HAAT) files to the tape. The 
HAAT files contain geometric correction parameters, 
and ancillary data such as scene path and row indices, 
date of data acquisition, and orbit number. The HAAT 
files plus the partially processed image data are re- 
ferred to as archive data or A-data, In summary, the 
data archive for each TM scene consists of HAAT files 
and commutated, aligned, radiometrically corrected, 
digital Image data stored in a band- interleaved-by-line 
(BIL) format on a 28-track high density tape (HDT-AT), 

TIPS creates a fully processed TM scene by resam- 
pling the image A-data for geometric correction on 
the basis of the SCD or GCD, The resampled image 
data are reformatted to a band sequential (BSQ) for- 
mat and written to a high density tape along with 
HAAT files. The HAAT files and fully processed image 
data are called P-data and a high density tape contain- 
ing these data is called an HDT-PT. 

Final TIPS products consist of radiometrically cor- 
rected archival data (A-data) and fully processed data 
(P-data) on computer compatible tapes (CCT’s). 
Radiometrically and geometrically corrected data in 
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P-data format are also available on 241 -mm film prod- 
ucts. The processed TM data are written directly from 
the HDT-AT's or HDT-PT’s onto CCT’s (CCT-AT’s 
and CCT-PT's). A laser film recorder Is used to 
generate film products from the HDT-PT’s. The CCT- 
AT’s and CCT-PT’s are shipped to the NOAA Landsat 
data archive at the EROS Data Center in Sioux Falls, 
South Dakota. These tapes are stored and duplicated 
for users. The TIPS film products are also shipped to 
the NOAA archive where they serve as masters for 
user film products. 

The Scrounge System employed a very similar se- 
quence of steps to correct TM data before July, 1983 
(Lyon et al., 1983). The only major difference is that 
image data on HDT-AT’s generated by the Scrounge 
System are in a band sequential, rather than a band- 
interleaved-by-line, format. 

TM DATA PRODUCTS 

Standard TM data products generated by the TIPS 


facility consist of digital Image data recorded on 
magnetic CCT’s and pictorial Image data recorded on 
film. The following discussion describes the format 
of the digital and film products generated by the TIPS 
facility at GSFC. Readers of this document should 
realize that the format of data products produced by 
other TM data processing facilities may be slightly difi 
ferent. The products described below are available to 
the general public and can be purchased from the 
NOAA Landsat data archive at the EROS Data Center 
in Sioux Falls, South Dakota, 

TM Digital Image Data Products 

The standard TM CCT product contains raultispec- 
tral image data in digital format for approximately one- 
quarter of a Rill 170 km along-track- by- 185 km across- 
track TM scene, A full scene of radiometrically and 
geometrically corrected TM data (P-data) consists of 
5965 scan lines with over 6000 pixels per line. A stan- 
dard scene quadrant contains approximately one- 
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Figure 11-15 Pixel and scan line dimensions of TM scene quadrants for P-data store on CCT's. 
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Figure 11-16 Pixel and scan line dimensions of TM scene quadrants for A-data stored on CCT’s. 


quarter of the pixels within a full scene, The number 
convention employed in identifying TM scene 
quadrants is presented in Figure 11-15 along with the 
pixel and scan line dimensions of scene quadrants. 

TM P-data is resampled onto a square pixel grid 
aligned with respect to true north before being 
recorded on magnetic tape. As described earlier, the 
orbit of the Landsat 4 and 5 spacecraft are Inclined 
98.2° with respect to the Earth's equator. Conse- 
quently, there are blank areas at the edges of full TM 
scenes and scene quadrants when they are projected 
onto the square pixel grid described above (see Figure 
11-15). Pixels with digital counts of zero for all spec- 
tral bands are assigned to these blank areas to pro- 
vide numerical values for all pixel locations within the 
square grid projection. These zero value pixels are 
commonly referred to as “fill pixels’’. Fill pixels are 
also employed in dividing TM A-data into scene 
quadrants for storage on digital CCT’s (Figure II- 1 6). 

An investigator is offered the following options 
when ordering a scene quadrant of digital image data; 


1) Extent of data processing: Tapes of partially pro- 
cessed image data (CCT- AT’s) or fully processed 
image data (CCT-PT’s) may be ordered. 

2) Cartographic projection: Investigators may order 
P-data registered to the Universal Transverse Mercator 
(UTM) or Polar Stereographic projections instead of 
the standard Space Oblique Mercator (SOM) projec- 
tion. 

3) Data resampling; Investigators may request that a 
nearest-neighbor resampling be performed instead of 
a cubic-convolution resampling for geometric 
correction. 

4) Digital image format: The digital image data may 
be obtained in either a band sequential (BSQ) format 
or a band-interleaved-by-line (BIL) format. 

5) Data density: The data may be packed onto a CCT 
at l600 bits-per-inch (bpi) or 6250 bpi. All of the TM 
data for one quadrant are contained on one 6250 bpi 
CCT. Three l600 bpi CCT’s are required per 
quadrant. 
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TM Film Products 

All TM film products available from the NOAA ar- 
chive are derived from the film masters generated by 
TIPS, Seven film masters, one per spectral band, are 
generated for each fully processed TM scene, Each 
master contains the black-and-white image of an en- 
tire scene at a scale of 1 : 1 ,000,000 along with annota- 
tion, tic marks showing cartographic coordinates, and 
gray scaies. 

The film masters are produced by a laser film 
recorder having an instantaneous spot size of a 0.0285 
mm square. The film masters are used to generate im- 
ages for investigators at the following three scales: 
1:750,000 (18.5 cm film); 1:375,000 (37.1 cm film); 
and 1:187,509 (74.2 cm film). All standard film prod- 
ucts represent an entire TM scene. Black-and-white 
film products may be obtained for each of the seven 
spectral bands. Investigators may also obtain color 
composites derived from three spectral bands. The 
investigator must specify which bands are to be 
depicted as red, green, and blue. Black-and-white im- 
ages are available as positive transparencies, negative 
transparencies, and positive paper prints. Color com- 
posite Images itre available as positive transparencies, 
positive paper prints, and 35 mm, positive slides. 


TM Data Availability 

Between the launch of Landsat 4 Quly 16, 1982) and 
the failure of its X-band transmitter (February 15, 
1983), over 6000 TM scenes were received by the 
Prince Albert, Canada, ground station and the 
Transportable Ground Station (TGS) in Greenbelt, 
Maryland. Scenes received by the Prince Albert sta- 
tion cover western Canada, the western United States, 
and a few areas in northern Mexico. Scenes received 
by the TGS cover eastern Canada, eastern United 
States, and the Bahamas. More than 1000 of these 
scenes are nighttime scenes consisting of thermal 
(TM6) data. Masuoka (1984) provides maps showing 
the number of Landsat 4 TM scenes acquired for each 
path-row location in North America. The utility of 
many of these images is limited by cloud cover and 
poor solar illumination conditions. 

After the launch of TDRS-East, Landsat 4 TM data 
were occasionally relayed through TDRS-East for 
engineering test purposes. Sixty TM scenes weie 
relayed on August 12 and 13, 1983, and approxi- 
mately 40 scenes were relayed on November l6 and 
18, 1983. The scenes covered areas in the United 
States, Africa, and Argentina. Although the Landsat 4 
TM/PF is still functional, the sensor has not been ex- 


tensively operated since the launch of Landsat 5 
(March 1, 1984). TM/PF data acquisition was sus- 
pended because of power problems aboard the Land- 
sat 4 spacecraft. The problems are discussed below. 

The Landsat 5 TM/F is currently acquiring data on 
a regular basis. Over 120 scenes per day are sent to 
TIPS at GSFC for processing. These data are either 
relayed by TDRS-East or directly received on the X- 
band frequency by the TGS or the Prince Albert sta- 
tion. Landsat 4 and 5 TM data received by the Prince 
Albert station are shipped to GSFC for processing as 
a result of an agreement between the United States 
and Canada. Processed data are sent from GSFC to 
the EROS Data Center for distribution, Canada also 
maintains an Independent archive of TM data received 
by the Prince Albert station. 

Only a small portion of the TM scenes sent to the 
GSFC facility are eventually processed by TIPS. The 
utility of many scenes is limited by cloud cover or 
low solar elevation during acquisition. These scenes 
are generally not processed unless specifically re- 
quested by an investigator. Of the over 6000 available 
Landsat 4 TM scenes, only 1900 scenes are currently 
scheduled for TIPS processing. 

An investigator may request the processing of any 
available TM scene through NOAA Landsat Customer 
Services. This organization may be contacted at the 
following address and telephone number: 

NOAA Landsat Customer Services 

Mundt Federal Building 

Sioux Falls, South Dakota 57198 

Telephone: (605) 594-6151 

This service organization offers assistance in the 
location of available TM scenes. An investigator can 
request a computer listing of all available TM scenes 
for a particular location by specifying geographic 
coordinates (e.g., latitude/longitude) or path/row 
coordinates. The list provides the following informa- 
tion: scene identification number, acquisition date, 
geographic coordinates, extent of cloud cover, and 
data product avaiIabi,Uty. Listings of available TM 
scenes may also be obtained on microfiche film cards. 
An investigator can initiate processing of a particular 
TM scene by ordering specific data products from the 
NOAA archive. 

The most readily accessible TM scenes are those for 
which fully processed digital data and film masters 
have already been archived by NOAA. As of June 20, 
1 984, TIPS had provided both fully processed data 
and a film master for 286 TM scenes to the EROS Data 
Center. Eighty-seven of these scenes were acquired 
by the Landsat 5 TM/F and the remaining 199 scenes 
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were acquired by the Landsat 4 TM/PF. In addition, 
the Scrounge System generated 287 scenes of fully 
processed TM/PF data on CCT’s prior to the Initiation 
of TIPS operations in August, 1983. 

Inquiries regarding the availablliiy of TM data 
received by foreign ground stations should also be 
directed to NOAA Landsat Customer Services with the 
following exceptions. The Canada Centre for Remote 
Sensing (CCRS) has an archive of digital data for over 
3000 TM scenes. The data were received by the Prince 
Albert station and were acquired primarily over 
Canada. Inquiries should be directed to: 

Mr, I. Press 

Landsat Production Coordinator 
Canada Centre for Remote Sensing 
2464 Sheffield Road 
Ottawa, Ontario 
KIA OY7, Canada 

The European Space Agency’s Earthnet Program 
Office has approximately 1200 scenes of TM data 
available. These data were received by the Kiruna, 
Sweden, and Fucino, Italy, stations and were acquired 
over Europe. Inquiries should be directed to: 

Earthnet User Services 
Earthnet Program Office 
Galileo Galilei 
00044 Frascati, Italy 

TM SENSOR PERFORMANCE AND 
DATA QUALITY 

NASA is conducting an extensive characterization 
of the orbital performance of the TM/PF and TM/F 
sensors and the quality of the data they are produc- 
ing, The characterization program, called the Land- 
sat Image Data Quality Analysis (LIDQA) program, in- 
volves investigators selected from the full breadth of 
the terrestrial remote sensing community. Results 
have shown that the Landsat 4 TM/PF performed 
beyond expectations and produced data that met or 
exceeded most prelaunch specifications. Preliminary 
results indicate that the Landsat 5 TM/F is also per- 
forming well and generating high quality data. 

Problems with the Landsat 4 spacecraft have cur- 
tailed operation of the Landsat 4 TM/PF. The first 
problem was the failure of the X-band transmitter. As 
described earlier, the X-band transmitter was origi- 
nally designed for direct readout of TM data to ground 
receiving stations within view of the spacecraft’s X- 
band antenna transmission pattern. Following the 
failure of this system in February, 1983, the only 


means of acquiring TM imagery was by the KU'band 
transmitter which was designed to relay TM data to 
the ground through TDRSS. The first TDRSS satellite 
(TDRS-East) was launched April, 1983 and it was not 
available for use prior to August, 1983. By this time, 
a second Landsat 4 technic.tl problem had emerged 
which placed additional restrictions on TM data ac- 
quisition. This problem involves the loss of electrical 
power from two of the spacecraft’s four solar panels. 
This problem has been attributed to the mechanical 
disruption of power cables onboard the spacecraft. 
Thermal stress within a cable potting compound has 
apparently caused disconnections in two of the elec- 
trical power cables connected to the satellite’s solar 
array, The remaining power is currently being used 
only to maintain the orbit and intermittently operate 
the MSS, If another panel should fail, the satellite will 
be lowered to a 500 km orbit to create the possibil- 
ity of repair and retrieval using the Shuttle. Landsat 
5 was modified before launch to prevent reoccur- 
rence of the problems experienced by Landsat 4, 

An adequate amount of Landsat 4 TM/PF data were 
received prior to the failure of the X-band transmit- 
ter to examine the orbital performance of the TM/PF 
sensor. The radiometric performance of the sensor 
is a characteristic of prime interest. Investigations have 
demonstrated that the relative radiometric responses 
of the indlviduai data channels are intercalibrated to 
a high degree of accuracy. The radiometric correc- 
tion procedure described earlier effectively accounts 
for radiometric response variations between channels, 
and only minor striping remains in radiometrically 
corrected data. The striping is usually on the order 
of the least significant bit (one or two digital counts) 
out of eight bits (256 digital counts) (Salomonson, 
1984). This residual striping is not visually evident in 
pictorial TM imagery unless extreme contrast 
enhancement is applied to data acquired over uniform 
areas. The absolute radiometric accuracy of TM data 
is still under study, but on the basis of early results, 
the absolute accuracy is believed to be within 
specifications (i.e., 10% of the full dynamic range 
within each spectral band) (Barker, 1984). 

Much of the residual striping may be attributable 
to within-scene variation in channel radiometric 
response. As described earlier, the channel radiomet- 
ric response functions are monitored at intervals of 
360 scan cycles (i.e., every 11,520 scan lines). 
Therefore, the procedure does not correct for varia- 
tion in radiometric response during acquisition of the 
approximately 5700 scan fines of data which con- 
stitute a scene. Several sources of within-scene 
radiometric variation have been identified (Barker, 
1984). 
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One important source of within-scene radiometric 
response variability is the anaiog-to-digitai (A/D) con- 
verters in the TM/PF. The A/D converters do not 
distribute the analog voltage from the detectors into 
256 equal Increments, In other words, the increment 
of radiance represented by a particular digital count 
does not necessarily equal tlie radiance increment 
represented by another count. This problem is re- 
ferred to as “bin radiance dependence" The ine- 
qualities range from radiance increments of nearly 
zero to twice the specified radiance increment and 
cause errors up to one count when relating absolute 
radiance to digital counts (Barker, 1984). Fortunately, 
the TM/PF digital response functions (digital counts 
versus spectral radiance) remained fairly linear over 
the entire 256 count range of the A/D converters. A 
least-squares fit of known spectral radiance from an 
integrating sphere to the prelaunch digital response 
of the TM/PF had a root-mean-square error of approxi- 
mately 0.5 counts (Barker, 1984). The inequal bln sizes 
caused local deviations from a perfectly linear 
response, but did not destroy the linearity over the 
entire response range, 

The effects of the A/D converter performance are 
apparent in partially processed data (A-data) by 
inspecting the frequency histograms of data from each 
detector. A particular digital value will occur much 
more or much less frequently than neighboring 
values. This effect is masked in P-data by the resam- 
pling procedure used to geometrically correct the 
data. 

An additional radiometric problem was en- 
countered with TM/PF thermal (TM6) data. The gains 
of the thermal band detectors are observed to 
decrease with time (Lansing and Barker, 1984). The 
decrease is attributable to condensation of water 
vapor on the window of the dewar containing the 
cold focal plane. Internal heaters were temporarily 
turned on and the sensor was outgassed to vaporize 
the condensed water and restore the thermal detec- 
tor gains. The radiometric correction procedure com- 
pensates for the loss of detector gain with respect to 
relative and absolute calibration. The loss, however, 
reduced the thermal sensitivity of the detectors. 

The geometric fidelity of fully processed TM/PF 
data has proved to be excellent. Early problems 
registering data from the primary focal plane (spec- 
tral bands TMl, TM2, TM3, and TM4) to data from 
the cold focal plane (TM5, TM6, and TM7) have been 
overcome. The early problems were attributable to 
a slight shift in the relative position of the two focal 
planes due to the stress of launch. The problems were 
corrected by adjusting the Focal Plane Geometry 
Matrices, described earlier, to account for the shift in 


the geometric correction procedure. Data are now 
registered to within 0. 1 pixel between any two bands 
(Salomonson, 1984), well above specifications (Table 
II-5). The temporal and geodetic registrations of TM 
data are also within specifications. An analysis of 
several corrected TM scenes by Brooks et al. (1984) 
demonstrated that the data consistently met specifica- 
tions for temporal registration (0.3 pixels) and usu- 
ally met geodetic registration accuracy specifications 
(0.5 pixels), Colvacoresses (1984) has demonstrated 
the ability to generate TM pictorial imagery which 
meets national map accuracy standards at 1:100,000 
scale. In comparison, MSS Imagery cannot generally 
meet accuracy standards at scales greater than 
1:200,000 (Salomonson, 1984). 

Before the launch of Landsat 4, much concern was 
directed towards the potential effects of high fre- 
quency structural disturbances, commonly called jit- 
ter, on the geometric quality of TM imagery. These 
disturbances are caused by the mechanical operations 
of spacecraft and sensor components such as antenna 
motors and scanning mirrors. These disturbances are 
sensed by the angular displacement sensors described 
earlier, which measure deviations in the TM sensor 
attitude occurring between frequencies of two to 125 
Hertz. Data from the angular displacement sensor are 
incorporated into the payload correction data (PCD) 
for geometric corrections. The concern regarding the 
effects of jitter have been alleviated by post-launch 
data analyses. The high band-to-band, temporal, and 
radiometric registration accuracies observed for cor- 
rected TM data (Brooks et al., 1984; Salomonson, 
1984) indicate that geometric processing effectively 
minimizes the effect of jitter on data quality, 

Markham and Barker (1984) have completed an ex- 
tensive evaluation of the spectral resolutions of both 
the TM/PF and TM/F, Their analyses were based on 
prelaunch measurements by Santa Barbara Research 
Center of the spectral transmission and response of 
individual components (l.e., optics, spectral filters, 
and detectors). They have determined spectral 
response curves, band edges (50% of peak response), 
and the spectral matching between the detectors for 
each band. Their studies showed that both sensors 
met all spectral response specifications for the reflec- 
tive bands with the following exceptions: (1) spectral 
bands TM2, TM3, and TM7 did not meet spectral 
flatness specifications in both sensors; and (2) the 
upper edge for TM5 was greater than the specified 
edge for both sensors. Table II-6 lists the specified and 
measured band edges. In addition, the spectral 
response curves for the reflective bands were observ- 
ed to be quite similar between the two models. 
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Table 11-6 

Calculated TM Spectral Band Edges for the TM/PF and TM/F (Markham and Barker, 1984). 


Spectral 

Band 


Specification 

TM/PF 

TM/F 

TM1 

Lower Edge (nm) 

450 + 10 

452 

452 


Upper Edge (nm) 

520 + 10 

518 

518 

TM2 

Lower Edge (nm) 

520 + 10 

529 

528 


Upper Edge (nm) 

600 + 10 

610 

610 

TM3 

Lower Edge (nm) 

620 + 10 

624 

626 


Upper Edge (nm) 

690 + 10 

693 

693 

TM4 

Lower Edge (nm) 

760 + 20 

776 

776 


Upper Edge (nm) 

900 + 10 

905 

904 

TM5 

Lower Edge (nm) 

1550 + 20 

1568 

1567 


Upper Edge (nm) 

1750 + 20 

1784* 

1784* 

TM7 

Lower Edge (nm) 

2080 + 30 

2097 

2097 


Upper Edge (nm) 

2350 -F 30 

2347 

2349 

TM6 

Lower Edge (nn\) 

10.4 + 0.01 

10.42 

10.45 

(Thermal) 

Upper Edge (am) 

12.5 -F 0.1 

11.6* 

12.43 


* Out of Specification 


The spectral resportses of the thermal band are not 
similar between the TM/PF and the TM/F, The upper 
band edge of the Landsat 4 TM/PF was much less than 
the specification (Table II-6). The reduced band width 
did not cause concern because the thermal sensitiv- 
ity of the band remained greater than specifications 
(Markham and Barker, 1984), The TM/F thermal band 
edges were within specifications. 

Evaluations of postlaunch Landsat 5 TM/F data have 
just begun. Early results indicate that TM/F data are 
comparable to TM/PF data with respect to radiometric 
and geometric fidelity. Evidence of problems with the 
TM/F A/D converters and with the thermal detector 
gains has again been encountered. These problems 
do not appear any more severe with TM/F data than 
with TM/PF data. Investigators can reasonably an- 
ticipate a reliable flow of high quality TM data from 
Landsat 5. 

Researchers have also begun to study the informa- 
tion content of TM data, Several investigators have 
noted the capability to use TM data, particularly the 
middle infrared bands (TM5 and TM7), to identify land 
cover categories which could not be accurately 


discriminated using MSS data (Degloria, 1984; Pitts et 
al., 1984), Other classification studies have shown that 
the prevalent per-pixel, maximum likelihood decision 
rules do not fully exploit the high TM spatial resolu- 
tion and they suggest the use of alternate approaches 
to classification (Williams et al„ 1984), The dimen- 
sionality of TM data in spectral data space has also 
been explored by eigenvector analyses. Results in- 
dicate that TM reflective band data typically contain 
four independent dimensions of useful terrestrial in- 
formation as compared to two dimensions in MSS data 
(Anuta et al., 1984; Bernstein et al., 1984; Crist and 
Cicone, 1984). The utility of TM data will become 
even more apparent as investigators extend TM data 
applications to a broad range of scientific research. 


USEFUL DOCUMENTS FOR TM DATA 
USERS 

This document presents a broad overview of infor- 
mation an investigator might need to effectively use 
TM data. Several references are recommended for 
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more detail and In-depth coverage of the relevant 
topics. The "Landsat 4 Data Users Handbook" 
describes the spacecraft, TM and MSS sensors, orbits 
and coverage, data communications, TM and MSS data 
products, and procedures for ordering data products. 
Forthcoming appendices to the handbook will 
describe TM data processing and CCT formats In 
detail. The handbook can be purchased from the 
following address: 


Distribution Branch 
Text Product Section 
U.S. Geological Survey 
604 South Pickett Street 
Alexandria, VA 22304 

Details on data processing and tape formats can also 
be found in the scries of NASA reference documents 
listed in Table II-7. A document de.scriblng the for- 


Table 11-7 

List of Landsat Reference Documents 


Accession Number 
N83-21469 

N83-21472 

N83-21479 

N83-21482 

N83-21473 and 
N83-21474 

N83-21466 

N83-15750 

N83-21467 

N83-21481 

N83-21468 

N83-21470 

N83-21480 

N83-21471 

N83-21477 

N83-21478 

N83-21475 


Document 

"Applications Notice for Participation In the Landsat-D Image Quality Analysis 
Program" 

"Landsat-D to Ground Station Interface Description" 

"Landsat-D World Reference System Users Guide" 

"Landsat-D Investigations Workshop" 

"Landsat-D Mission Operations Review" 

"LGSOWG CCT Format CCB Document: The Standard CCT Family of Tape 
Formats” 

"Spectral Characterization of the Landsat-D Multispectral Scanner Subsystems" 

"Landsat-4 MSS Radiometric Characterization Report" 

"Partially Processed Multispectral Scanner High Density Tape (HDT-AM)” 

“Landsat Multispectral Scanner Computer Compatible Tape Format" 

"Landsat-D Assessment System Library Computer Compatible Tape" 

"Thematic Mapper High Resolution 241mm Film" 

"Landsat-D Project ADDS/LAS Scrounge Interface Control Document" 

"Landsat-D Data Format Control Book, Volume VI, Appendix A, Partially Processed 
Thematic Mapper High Density Tape (HDT-A7)" 

"Landsat-D Data Format Control Book, Volume VI, Appendix D, Thematic Mapper 
Computer Compatible Tape (CCT-AT, GGT-PT), Format Specification: 

"Landsat-D Accelerated Payload Correction Subsystem, Output Computer Compatible 
Tape Format" 
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mat of TM CGT-AT’s and CCT-PT‘s Is particularly 
useful anti Is listed under accession number 
N83-21478 (Table II-7), 

These documents can be obtained from the follow- 
ing address: 

National Technical Information Service 

5285 Port Royal Road 

Springfield, VA 22161 

Results of the NASA UDQA program are summa- 
rized in a two volume document, “Landsat-4 Science 
Investigations Summary” (NASA Conference Publica- 
tion 2326). This document contains valuable informa- 
tion on data processing and quality. Selected articles 
from LIDQA investigations are also present in a special 
issue of lEBE Transactions on Geoscience and 
Remote Sensing (Volume GE-22, Number 3) published 
in May, 1984. These articles include useful reports on 
the Information content of TM data. 

Up-to-date information on the status of the Land- 
sat satellites, sensors, and data can be obtained from 
the “Landsat Data Users Notes" published quarterly 
by NOAA. Subscriptions to this publication may be 
obtained at n< > charge by contacting NOAA Landsat 
Customer Services at the address or phone number 
provided earlier. 
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SECTION III 


PROSPECTIVE THEMATIC MAPPER RESEARCH 
APPLICATIONS IN THE EARTH SCIENCES 


VEGETATION AND SOILS SCIENCE 
Introduction 

Vegetation and soils are principal components of 
the Earth’s biosphere, a cosmic description of the liv- 
ing organisms found on our planet. Vegetation and 
soils are continually changing in form and distribu- 
tion over a wide range of time scales. Yet they play 
a pervasive role in the dynamics of the biosphere. 
Through photosynthesis, green vegetation manufac- 
tures food and fiber, and it releases the oxygen essen- 
tial to man and all other living organisms. In counter- 
balance, the process of respiration causes the release 
of organically-derived energy critical to biological 
growth, and it liberates carbon dioxide to the at- 
mosphere. Soils provide a reservoir of energy and 
essential elements, and the physical substrate to sup- 
port these metabolic processes. 

In addition to these primary molecular processes, 
vegetation and soils play an integral role in 
blogeochemical cycling, global and local energy 
balance, and the hydrologic cycle. The type, condi- 
tion, and distribution of the vegetation/soil complex 
determines the habitat and distribution of the Earth’s 
animal population, as well. Vegetation and soils are 
also among the most sensitive environmental 
parameters responding to the activities of man and 
changing climatic conditions. 

As changes on the Earth’s surface increase with in- 
creasing human population and advancing 
technology, the opportunity to investigate the natural 
functioning of terrestrial ecosystems may become in- 
creasingly limited. It is important to understand the 
natural relationship between plants and soils In dif- 
ferent environmental settings to better understand the 
response of ecosystems to the activities of man. In 
developing this understanding, knowledge of the 
distribution and condition of vegetation and soils in 


time and space is critical. Because of its unique 
discrimination capabilities and synoptic view, the 
Thematic Mapper (TM) may provide a new means to 
Improve our understanding of tlie role of vegetation 
and soils In the biosphere. 

Background 

Remote sensing systems placed on Earth-orbiting 
spacecraft are uniquely able to conduct surveys of ter- 
restrial surface conditions at uniform scales and for- 
mats on a globally repetitive basis. Our current 
understanding of the utility of such surveys Is based 
largely upon the analysis of Landsat Multispectral 
.Scanner (MSS) data coile-cted by the first generation 
of Landsat .spacecraft. Analysis of MSS imagery has 
shown that multispectral surveys of reflected solar 
radiation can be used to Identify the type, condition, 
and stage of growth of various forms of surface 
vegetation. MSS surveys can also be used to 
discriminate various types of .soils on the basis of their 
compositional and textural properties. These earlier 
studies have shown that the accuracy of the informa- 
tion derived from MSS imagery can vary significantly 
depending upon local relief, solar illumination con- 
ditions, and seasonal variations in the form and nature 
of the vegetation-soil complex. Certain types of in- 
formation such as the growth stage or phenol ogic 
development of individual plant species can best be 
determined through the analysis of multitemporal 
MSS imagery. 

The measurement capabilities of the TM are 
significantly more diverse and sophisticated than 
those of the MSS. The TM holds great promise for ob- 
taining more accurate assessments of many of the 
biospheric parameters described above. In addition, 
the TM potentially provides a capability to obtain 
wholly new types of information concerning the 
nature and condition of surface vegetation and soils. 
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Our current understanding of the utility of TM data 
for studies of the vegetatlon«soil complex is based 
upon the analysis of data collected In the past by air- 
borne sensors that simulate the general measu .ement 
capabilities of the Landsat TM. In addition, a number 
of Investigators have employed TM Imagery obtained 
by landsat 4 to compare the utility of MSS and TM 
data for specific types of Earth observations. On the 
basis of these earlier studies, the following observa- 
tions can be made concerning the use of TM data for 
vegetation and soil investigations: 


a) Far more classes of vegetation can be 
discriminated than ever before during certain 
times of the year. This is true both for natural 
and cultivated vegetation and has been at- 
tributed particularly to the new shortwave in- 
frared bands (TM 5 and 7) and the Improved 
spatial resolution of the TM. It appears that 
several common southern pine species can be 
spectrally differentiated, for example, as can 
sugar beets and alfalfa (see Figure III-l). The 
enhanced spectral and, in particular, spatial 
resolution of TM data allows the discrimination 
of various forms of riparian vegetation and small 
agricultural fields, The identification of these lat- 
ter features Is achieved by invoking an 
understanding of their geometric shape and 
geographic context In the analysis of a TM 
image. 

b) Patterns correlated with soil type are clearly 
discernible in many images, even in areas con- 
taining dense vegetation. 

c) The thermal detection capability of the TM is 
useful In the discrimination of some vegetation 
types as, for example, tree identification in 
forests and crop differentiation, such as soy- 
beans and sorghum, which can not be achieved 
by the use of the visible and short wave infrared 
bands alone. 

d) The spatial resolution of the TM provides far 
more information on the variability of the land- 
scape, particularly In areas disturbed by man and 
those that have been subjected to dramatic and 
sudden climatic changes. 

e) Geometric fidelity of the TM data is excellent, 
implying that high quality Image rectification 
and multitemporal registration can be readily 
achieved. 

f) The radiometric quality of the TM data exhibits 
superior stability and sensitivity suggesting an 
unparalleled capability to derive measurements 
of absolute Earth radiance on a repetitive basis. 
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figure IIl-J 

This graph illustrates the relative spectral response 
of different evergreen tree species observed in a Land- 
sat 4 TM image of northern Virginia acquired in the 
fall of 1982, These species have a relatively similar 
spectral response in the red portion of the visible spec- 
trum (TM band 3), but markedly different spectral 
characteristics in the near infrared (TM band 4), The 
level of species discrimination achieved in this scene 
is due in part to the improved spatial resolution of the 
TM. (Illustration courtesy of Nancy Milton, U, S. 
Geological Survey). 

The contribution that TM data might make to the 
enhancement of surface observations, the refinement 
of existing models, and the development of altogether 
new ones Is yet to be tested, Certainly, Improvements 
in the identification and measurement of both natural 
and cultivated vegetation and associated soils could 
have a significant impact on the study of animal 
habitats, the condition of arid lands, soil moisture and 
water holding capacity, and surface heating and cool- 
ing. The same direct measurements of the surface may 
serve, by Inference, to support the development of 
more definitive models of land processes that will 
ultimately increase our understanding of ecological 
Interactions at both regional and global scales, 
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Prospective Research Themes 

The purpose of the following discussion is to 
stimulate creative thinking about various topical prob- 
Icms in the field of biospheric research that could 
potentially be addressed in new and innovative ways 
through the analysis of orbital TM imagery. This 
discussion is not intended to be a comprehensive 
compendium of all of the research topics that could 
potentially be conducted, nor is it intended to repre* 
sent a prioritization of research topics to be explored 
by NASA in the future. Rather, the goal of the follow* 
ing dlscu,ssion is to illustrate generic types of scien- 
tific investigations that are potentially amenable to TM 
image analysis, 

The following dlscu.sslon has been divided into 
three sections that are broadly concerned with studies 
of terrestrial ecology, the response of ecosystems to 
environmental change, and fundamental biospheric 
processes. A series of topical problems is presented 
in each section. One should note that these problems 
arc not mutually exclusive nor are they categorized 
by traditional scientific discipline. Rather, they serve 
to Illustrate generic processes that have a multidisci- 
plinary relevance. 

Prospective investigators should carefully consider 
the identification of specific problems and the 
development of experimental methodologies in rela- 
tion to the spatial and temporal measurement capa- 
bilities of the TM which are described in Section II. 
If current surface conditions in one geo/jraphical loca- 
tion can be used to simulate the predicted conditions 
in another geographical location for some past or 
future time, such as in the study of soil erosion, for 
example, then the constraints on the study due to 
temporal requirements for TM data could potentially 
be alleviated. Other kinds of data may also be enlisted 
to satisfy the need for Information that is temporally 
and/or spatially specific and significant to the conduct 
of an investigation. Other remotely sensed data, such 
as that which exists from the past twelve years of 
Landsat MSS data collection, are available. Also, other 
records and sources of data and maps should be 
considered. 

Terrestrial Ecology 

The term “ecology” is derived from the Greek 
word "oikos” meaning house. In the most literal 
sense, the study of terrestrial ecology involves the 
characterization and understanding of the system that 
houses mankind and ali the other forms of life on our 
planet. The availability of TM data on a near global 
basis provides an opportunity to observe ecological 


relationships between regional flora, fauna, and soils 
in a wide range of environmental .settings. In fact, the 
TM provides an opportunity to conduct comparative 
studies of ecological conditions and relationships in 
selected regions at a level of detail that has not been 
possible in the past. The observational capabilities of 
the TM can also be exploited to .study the relation- 
ships within and between individual ecosystems, 
leading to an Improved understanding of the global 
ecology of the Earth. A variety of topical research 
issues in this field are presented below. 

Plant Phenology, Commutlty Maturation, 
and Ecological Succesion, Plant succession is 
basically a community process in which the structure 
and function of individual plant species changes over 
time. All types of vegetation undergo an inherently 
natural kind of change, ranging from the rapid 
phenological development of annual crops, to the 
decadal growth of trees, to the still slower processes 
of community maturation, to the process of succes- 
sion, which may take thousands of years. These 
changes are complex functions, many as.soclated with 
changes in localized and/or regional climate which af- 
fect vegetation and the soil substrate. TM imagery may 
be instrumental in testing models of plant succession 
which require knowledge of the distribution of In- 
dlvlaual plant species and their associations in dif- 
ferent euvlronments. An analysis of TM imagery of 
multiple geographic sites representing various 
ecological seres of a successlonal paradigm may 
enable us to learn more about ecological succession 
without having to depend on long term observations 
of a specific area. Such an approach might be applied 
to the study of forest succession, for example. 

A forest is a complex biological community in a 
continual state of flux. The forest community is com- 
posed of numerous components, among which are 
trees, shrubs and herbaceous vegetation, as well as 
insects and higher order animals and the micro-flora 
and -fauna of the underlying soils. These components 
interact over time, and their relative health and abun- 
dance change in ways that are often predictable in ac- 
cordance with the laws of natural selection. The suc- 
cessional processes that govern the evolution of the 
different components of this community are poorly 
understood at the present time. TM data may provide 
valuable information in cliaracterizlng both autogenic 
succession (the natural replacement of one biotic 
community by another in the absence of catastrophic 
disturbances) and allogenic succession (replacement 
resulting from an external force such as man, floods, 
etc.) within the forest ecosystem. 

In the study of wetlands ecology, TM imagery may 
allow us to improve our understanding of community 
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maiuraiion and ecological aucccMion by better condition, diMribuiion, and areal extent (»ee Figure 

underManding the coupling of the biological and III-2). AsMKlation^ of given xpecieft and changes in 

phymcal environment. TM data are expected to their compoMtion are meaningful environmental in- 

enhance the Identification of vegetatum indigenous dicators. Improved observations should increase the 

to wetlands and improve the measurement of their inlerential information available to better understand 
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This is a land cover classification map of Peelfoot separabU.ty that has been achieved amonff the tree. 
Lake, Tennessee based upon the statistical analysis of grass, and aquatic species present in this area The 
a Landsat 4 TM image acquired in August, l^)H2. The portion of the lake outlined in this classification map 
classification of different forms of natural vegetation is roughly 6 x 13 kilometers in size. (Illustration 
(Kcurring in this wetland environment is based upon courtesy of Dale Quallrochi, NASA— National Space 
variations in spectra! response observed in TM bands Technology' Laboratories). 
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basic questions about water inundation patterns, soii 
chemistry, and plant growth and productivity, and 
the interdependence of these plicnomcna. For exam- 
ple, a current nuzzle to many wetlands scientists is 
the ca js^e for growth of two distinct forms—the short 
and tail forms—of Spartina altmUJJora, a dominant 
salt marsh grass species, in what are thought to be 
non-limiting nutrient environments. Unraveling this 
puzzle may help to understand the factors Influenc- 
ing plant productivity and its role in the estuarine 
trophic chain. 

The natural succession of plant communities is nor- 
mally characterized by an increase in species diver- 
sity, the structural complexity of individual species, 
and total community biomass. In wetland environ- 
ments, the local accumulation of organic matter that 
naturally accompanies community succession can 
produce significant changes in the local hydrologic 
regime. The accumulation of organic matter can serve 
to trap suspended sediments leading to a cot.* oml- 
tant decrease in water depth. As the local water level 
decreases, the frequency and depth cf water inunda- 
tion decreases so that less organic matter is removed 
during sub.sequent runoff cycles. TM imagery can be 
used to monitor spatial and temporal variations in 
species distribution and water depth within the 
wetlands environment. Observations of this nature 
could directly contribute to studies of the factors 
governing community succession within coast 
estuaries. 

Soii Genesis. Soils are an integral component of 
vegetation communities. Soils form a structured 
matrix of minerals and organic matter which anchors 
and physically supports vegetation while serving as 
a reservoir of plant nutrients and moisture. Soil 
formation results from the interaction of climatic 
conditions, living matter, and geologic materials, as 
conditioned by regional variations in topographic 
relief and the activities of man. TM data can be used 
to map regional vegetation patterns which are a 
consequence of and factor in soil development. As 
discussed In subsequent sections, TM data can also 
be used to characterize and map geologic parent 
materials and geomorphoiogic features which play a 
significant role in soil formation. The TM provides the 
spatial •'..'olution and synoptic view required to relate 
spatial variations in these soil forming factors to the 
variability of natural soils as depicted in soil maps. 

Conventional soil maps divide the natural soil 
continuum into discrete units with sharp boundaries. 
In reality, individual soil units are rarely separated by 
sharp boundaries. Soil characteristics commonly vary 
both within and between units. Soil scientists are 


beginning to use matltematlcal descriptions of spatial 
variability to better characterize the i.'ontinuum of soil 
units commonly encountered In nature. The use of 
TM data to remotely measure soil parameters would 
greatly facilitate the gathering of sufficient contiguous 
data to support the development of such 
mathematical techniques. Image texture, which is the 
characterization of image spectral data within the 
spatial context of the image, may be useful for 
defining soil boundaries over broader scales. The 
discrimination of soil types from TM imagery, 
particularly for areas where soil maps are lacking, 
when combined with other kinds of environmental 
data such as geologic, topographic, and climatic data, 
may advance our understanding of soil genesis. 

TM data may be of particular value for studies of 
remote areas or areas otlierwlse lacking detailed soil 
information. Detailed .soil maps are already available 
for much of the United States, but detailed maps are 
not available for many other regions around the 
world. Furthermore, there is a lack of standardization 
among the maps which are available. TM data offers 
the potential for discrimina^i.g different soil types 
either by direct observation of bare soil or by in- 
ference from observations of vegetation, geomor- 
phology, and geologic materials. Information un the 
different soil types may then be applied to studies of 
bioproductivity, climatology, and hydrology. 

Ecological Competition. Certain biological and 
physical pressures induce members of a plant 
community to exhibit complementary growth 
patterns. Such a response is characteristic of ecological 
competition, but the true cause and effect relationship 
is not well understooa. This applies to the biological 
interaction among members of a plant community 
and between cornmunitie as well. Plants have been 
found to emit chemical substances as a defensive 
mechanism for survival. The soil environment can act 
as a carrier of such substances. The high spatial, as 
well as spectral, resolution capabilities of the TM may 
provide valuable insight into the study of 
phytosociology and allow interferences to be made 
about the underlying blochemical/phy cial processes. 

As an example, knowledge of the distribution and 
condition of individual tree species within a forest 
canopy can contribute significantly to an 
understanding of the ecological competition that 
occurs among the forest flora. Th&spatial variability 
of individual species is an important factor in 
understanding how forest members relate to one 
another in both natural and managed situations. TM 
data may also be used to detect variations in the 
density and thermal characteristics of forest canopies. 
Information of this nature could potentially be used 
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to study thr response of individual species and plant 
communities to various tvpes ol strevs (e g. air 
pollution and insect infestation). 

1'.M data analysis may contribute to a different level 
of understanding about the forest system by 
addressing the relationship between forest 
community morphology and environmental 
influence. The distribution of forest plant 
communities is determined by genetic, as well as 
environmental influences If .such patterns re.sult in 
a characteristic surface texture as expre.s.sed in 
multi.spectral TM imagery, it may be possible to infer 
the physiological and morphological factors that 
g«>vern the development and proliferation of specific- 
plant communities. 

A cixssic phenomenon in the study of ecological 
competition is the ca.se of species invasion, for which 
the underlying caases are n..'t well • JersUMKi. Semi- 
arid regions are particula'I) ^usce)/ >!.■ to .species in- 
va.sion. Hy virtue of such p;./sical charucteri.stics as 
li)w and erratic precipitation, rougti io|>>graphy, p<K>r 
drainage, or cold temperatures, semi arid regions are 
typically unsuited to cultivation and are characterized 
by a relatively low density of indigenous plant .species. 
The improved spatial re.solution and radiometric 
sen.sitivity of the TM could prove to be particularly 
u.seful in monitoring subtle variations in the density 
and distribution of vegetation in extended .semi-arid 
regions (.see Figure lll-A) Temporal trends and .spatial 
gradients in surface vegetation can be used to detect 
ItK'al areas of incipient desertification in such regions 
In principle, the TM’s .spectral capabilities could be 
used to monitor the growth and distribution of key 
species such as juniper LJuniperas accidentalisj which 
reduce the susceptibility of surface .soils to wind ero- 
sion. The invasion of .specific species in such envi- 
ronments may al.so produce diagnostic variations in 
surface thermal conditions which could potentially 
be .scn.sed by the TM. 

Animal Habitats. Plant communities control the 
types and numbers of animals pre.sent in a region, not 
only by contiolling the finxl supply, but al.so by 
providing shelter from weather and predation. The 
animal population, in turn, influences the IcKal 
distribution of vegetation through selcctiv; herbivory 
and the transportation of seeds. The balance between 
fiHKl supply and ani^^al populations in different 
environments is a < • tvearch topic in the .study 
of animal habitats i • - ts can be identified on the 
basis of the type, diveisity, and amount of vegetation 
present, and tend to be concentrated along interfaces 
between vegetative communi ‘ ». Pattern recognition 
algorithms that can map the ''.^ometric shapes and 
boundaries of plant communities could prove to be 


extrentely useful in delinineanng the sha|X‘ and extent 
of animal habitats 


Fifture iH-.i 

Airborne multispeclral scanner image of Lisbon 
Valley, Utah in which sfwctral variations in the red 
/Hjrtion of the visihle spectrum (corresponding to TM 
band 3) have been used to infer relative variations in 
vegetation density ihe area displayed in the image 
is appro.ximately 5 x 9 kilometers in size. The scale at 
the bottom of the figure denotes the areal density of 
surface variation in fwreent Vegetation density is 
highly correlated with surface geology in this semi-arid 
environment. For example, many of the areas 
containing 20% vegetatii-e coivr or less are correlated 
with the surface expression of a regional sandstone 
unit which is outlined in white in the figure. 
(Illustration courtesy of James Conel, Jet Propulsion 
l.ahoratoiy). 

Certain wildlife species make extensive use of the 
forest periphery for feeding and protectiem. For ex- 
ample, elk require four bxsic vegetative habitat 
categories to .survive in rangeland environments: 
hiding co /er, winter thermal cover, summer thermal 
cover, as well xs forage. Maintaining forage-to-thermal 
cover ratios and hiding-to-forage cover ratios within 
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specific distances Is critical to elk survival. The use 
of TM imagery to define the edges of forest plant com- 
munities could make an important contribution to 
studies of wildlife movement and distribution. The 
improved spectral and spatial resolution, including the 
thermal detection capability of the TM, may con- 
tribute to an improved understanding of rangeland 
habitats and the quality and quantity of life that can 
be supported in specific t\. ironments. 

Detection of Extraordinary Plant Species. The 
evolutionary development of plant species has 
depended on climatic and edaphic factors, dispersion 
mechanisms for plant germplasm, and competitive 
pressures among living organisms which rely on 
plants for nutrition. Globally synoptic views of the 
Earth’s natural plant communities, using the improved 
resolution capabilities of the TM, may significantly 
enhance the current understanding of species 
phylogeny, Many treacherous and remote areas of the 
world remain essentially unexplored at the present 
time. By virtue of their isolation, these regions may 
contain new and as yet unknown plant species. The 
use of TM data in a predictive or inferential mode may 
provide a unique capability to identify areas with 
anomalous spectral and textural characteristics which 
may contain extraordinary species. The discovery of 
such species could lead to major advances in the 
sciences of genetics and ecology. Such a prediction 
mode might involve the incorporation of TM data in 
descriptive models of environments where genetic 
anomalies are likely to occur. 

Integrated Functioning of Complex 
Ecosystems. The concept of holism represents the 
view that an integrated ecosystem can behave in 
certain ways that are independent of and greater than 
the additive behavior of its component parts. For 
example, in a river-floodplain system, the river, the 
floodplain, and the vegetation and soils represent 
hierarchical components of progressively lower 
order. The characteristics of the river floodplain 
system are distinct from those of any of its 
components. In a similar manner, the processes of a 
coastal estuary that result from the interaction of 
wetlands and the nearshore marine environment 
become more than the processes associated with low- 
lying vegetative growth and tidal flooding. The 
higher-order holistic processes of an estuary may be 
viewed in terms of the exchange of nutrients between 
the land and ocean within an area of active 
sedimentation. 

In the holistic study of a system, it may not be 
necessary to study the functional parts in fine detail 
to gain insight into their mutual interdependence. 


Holistic studies are concerned with the dynamic 
relationship between ecosystem components for the 
purpo.se of obtaining improved understanding of tlie 
functioning of the overall system. Holistic studies can 
be conducted on many scales. Synoptic TM imagery 
acquired on a repetitive basis can potentially support 
such studies from the local to the regional scale. 

Ecosystem Response to Dynamic 
Change 

The response of ecosystems to changing 
environmental conditions can be exceedingly 
complex. Certain components of the integrated 
system may be altered in diagnostic ways by specific 
environmental changes. Identification of the system 
components that manifest the greatest short term 
change in response to changing environmental 
conditions is especially important.These components 
provide a basis for identifying ecosystems that may 
be experiencing major changes in form and function. 

Long term studies of ecosystems that are subjected 
to changing environmental conditions can provide 
considerable Insight into the physical processes that 
connect different system components. In particular, 
many of the more subtle processes occurring within 
a system may become more apparent or obvious 
when the system as a whole is striving to reach a new 
state of Internal equilibrium. Ecosystems responding 
to changing environmental conditions can serve as 
natural laboratories for studies of the relationships 
between system compc'^'cnts. Such studies are 
particularly important in attempting to anticipate the 
impact of man’s activities upon the condition and 
evolution of individual ecosystems. This section 
discusses several topical problems related to the study 
of ecosystem response to dynamic change, 

Botanical Stress. Plants are subjected to many 
kinds of stress including droughts, floods, nutrient 
deprivation, insect infestation, thermal stress, soil 
toxicity, disease, and animal predation. The response 
of certain plants and communities to various types of 
stress has been studied in some detail, particularly for 
agricultural species of economic interest. However, 
current understanding of how stress conditions are 
manifested and ameliorated in plants is far from 
complete, even for many major crop species. Our 
knowledge of the type and magnitude of stresses 
experienced by natural vegetation is even more 
rudimentary at the present time. In fact, our limited 
ability to detect and delineate stress in natural 
vegetation is based upon our ignorance of what 
constitutes normal conditions for individual species 
in different environmental settings. TM observations 
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can be used to obtain selected types of Information 
concerning the density, distribution, and condition 
of plant species whicii could potentially be employed 
in studies of plant stress (see Figure III-4). 
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are normally subjected to the same weather condi- 
tions, Nonetheless, substantial variations in plant 
biomass, stage of development, and yield occur 
within agricultural fields. While remote sensing does 
not offer the possibility of directly assessing yield, it 
provides the most practical approach available to 
assessing variability in biomass and stage of develop- 
ment as they relate to healthy and stressed conditions 
within fields. Since these factors can be quantified, 
it permits the use of each field as a “natural” 
laboratory to understand environmental effects on 
crop growth. In principle, such experiments would 
enable plant ecologists to determine the range of 
variability in plant morphology and biomass produc- 
ed under healthy conditions in a controlled environ- 
ment. Information of this nature could, in turn, be 
used to establish criteria for identifying abnormal con- 
ditions in the morphology or biomass of selected 
species that commonly occur In nature. The ability 
to recognize abnormal conditions in natural vegeta- 
tion canopies, and to relate the occurrence of such 
conditions to local environmental conditions could 
provide considerable insight into the causes, 
manifestation, and long term effects of stress in natural 
environments (see, for example, Figure III-5). 


This graph displays variations in oak leaf reflectance 
and soil geochemistry observed within a forested test 
site in northern Virginia during September, 1980. Leaf 
reflectance measurements were obtainea with a por- 
table field spectrometer over a spectral interval cor- 
responding to TM bands. Soil samples were obtained 
at the base of individual trees that were selected for 
detailed observation. The general trend of the plotted 
data suggests that oak leaf reflectance increases with 
increasing concentration of lead and copper in the 
local soils during this stage of the growing season. This 
trend may result from early leaf senescence in areas 
characterized by toxic soil conditions. (Illustration 
courtesy of Mark Labovitz and Robin Bell, NASA — 
Goddard Space Plight Center). 

The characterization of vegetation stress could 
potentially benefit from the use of agricultural fields 
as controlled test sites. For example, these fields could 
be used to improve our understanding of the effects 
of air and mineral pollution on plants, and the natural 
variability of healthy vegetation at the most elemen- 
tary level. An individual crop field is ordinarily sub- 
jected to the same management practices throughout, 
which can be either beneficial or deleterious to 
growth. The field usually contains a single crop, 
planted nearly simultaneously with uniform spacing, 
with similar applications of fertilizer and other 
treatments. In addition, contiguous agricultural fields 


Soil Destruction. Soil destruction by erosion, 
laterization, or other processes can be exceedingly 
rapid compared to soil genesis because the former 
phenomena are generally attributed to dramatic 
changes in the physical environment whereas the 
latter phenomenon occurs over geological time scales. 
Soil distribution historically has been a major problem 
for countries undergoing agricultural and commercial 
development. Tropical forest soils, in particular, are 
very fragile. Certain agricultural techniques currently 
in use actually induce erosion and accelerate nutrient 
depletion. Remote sensing data may be used to map 
existing soils directly or by inference from the 
discrimination of vegetation types. TM data may 
specifically be used to monitor or predict soil erosion 
by assessing changes In soils and vegetation, hence 
providing more accurate inputs to soil erosion 
models. 

Most of the commonly used models of soil erosion, 
such as the Universal Soil Loss Equation and Chepil’s 
wind erosion equation, predict the average annual 
erosion from single, uniform fields. Although 
adequate for the recommendation of field-by-field soil 
management practices, the models cannot be directly 
extended to large, heterogeneous areas or to 
individual soil erosion events. These models are often 
based on empirical observations of erosion that 
occurs on small plots. They do not account for the 
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This graph displays temporal variations in the spec- 
tral response of two different crops during the 1^82 
growing season. The “greenness" parameter displayed 
in the graph is a statistical transformation of the spec- 
tral response of TM bands 1, 2, 3, 4, 5, and 7 which 
can be correlated with the visual greenness of the 
Earth’s surface. The data points appearing in the 
graph are derived from airborne and orbital TM 
surveys conducted over Webster County, Iowa. The 
continuous curves represent modelled variations in 
surface “greenness" anticipated for corn and soybean 
crops. Changes in the “greenness " of agricultural 
fields as a function of time have been used in the past 
to discriminate different crop species. In principle, 
multitemporal variations in canopy reflectance could 
be used in the future to detect stress conditions in both 
natural and cultivated vegetation. (Illustration 
courtesy of Gautam Badhwar, NASA— Johnson Space 
Center.) 

total variability of soil conditions and tillage practices 
encountered in different geographic areas, nor do 
they account for temporal chang.s in surface cover 
and climate, Also, such models do not always reflect 
a theoretical understanding of the physics of the soil 
erosion process, Consequently, their applicability to 
different environmental settings is uncertain, 

TM data may aid in the development and valida- 
tion of theoretical models which extend the 
understanding of soil erosion to larger geographic 
areas and a wider variety of soil conditions and 
management practices. Through the use of the 
shortest TM wavelength band (band 1), it may be 
possible to estimate the sediment bedload of flowing 
rivers. Suspended sediments in running water are a 
primary indicator of soil loss from selected regions, 
Soil erosion in semi-arid environments is intimately 
related to spatial and temporal variations in the local 


abundance of surface vegetation and soil moisture. 
Some forms of arid land vegetation serve to stabilize 
the local soil surface and deter wind erosion. In other 
instances, arid vegetation serves to disturb the ground 
surface, locally exposing underlying deposits of sand 
and gravel to aeollan erosion, The distribution and 
density of different plant species Is critically depen- 
dent upon the ability of local soils to retain moisture, 
Multispectral TM imagery can be used to monitor 
variations in canopy closure and biomass in arid 
regions over a variety of time scales. Such observa- 
tions could be related to ground based observations 
of soil transport to develop improved models of soil 
erosion potential in arid regions. These models could 
subsequently be employed to identify areas of incip- 
ient desertification In semi-arid terrain. 

Short Term Environmental Change. Immediate 
environmental change is occurring in many parts of 
the world due to the activities of man and the effects 
of storms, fire, and flooding. The destruction of 
tropical forests and desertification of the Sahel have 
received popular attention, but equally profound 
Changes are occurring in many other regions, such as 
the rangelands of Texas and New Mexico. TM data 
may be able to provide a new capability for measuring 
the impact and consequences of such events. 

The enhanced spatial resolution of the TM may be 
particularly useful for detecting incipient environ- 
mental change. For example, the TM thermal detec- 
tion capability may provide a unique means of Iden- 
tifying and measuring the effects of desertification 
over large areas, Semi-arid regions which are prone 
to desertification are commonly characterized by high 
spatial and temporal variability in surface cover 
conditions, probably more than any other ecosystem. 
Erratic precipitation, animal grazing patterns, the 
appearance of ephemeral plant species, and the use 
of defoliants all contribute to this variability. Spatial 
and temporal variations in vegetation distribution, 
moisture conditions, and soil distribution produce 
commensurate variations in the spectral reflectance 
and thermal emission properties of semi-arid terrain. 
The ability to measure and monitor this variability 
with the TM will help in understanding arid 
ecosystems and their response to environmental 
change. 

In wetland environments, the enhanced spectral 
discrimination capability of the TM should improve 
our ability to detect variations in water quality, 
planktonic biomass, and macrophytic biomass. All of 
these parameters are sensitive indicators of en- 
vironmental change. In wetlands, the land and water 
are coupled through a number of complex biological, 
geological, and physical processes. For instance, an- 
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tliropogenic activity may cause changes In river flow 
or in watershed vegetative cover, which can result in 
a dramatic change in the supply of silt and nutrients 
to the downstream environment. This alteration of 
the water regime may have profound impacts on 
biological activity downstream, both in the wetlands 
themselves and adjacent estuarine waters. Therefore, 
an improved understanding of the water regime, 
specifically the relative rate and distribution of sur- 
face water movement through the wetland, may be 
an area to which TM research may contribute. 

By improving our ability to differentiate species and 
vegetation density in wetlands, TM data may provide 
insight into overland waterflow rates that can substan- 
tially Increase our present knowledge of surface water 
movement. TM data may also enhance the detection 
and definition of drainages and the land/water inter- 
face. The TM thermal band may provide Inferential 
information regarding the mixing of fresh and 
saltwater within wetlands, depth of small water 
bodies, eutrophication levels, and potential sources 
of thermal loading. Contextural information derived 
from TM observations may provide a totally new 
capability for understanding the geornorphological 
development of wetlands, as the degree of water body 
dispersion and connectivity of the wetlands signifies 
stages of ecosystem development directed by 
response to environmental change. All of these 
capabilities are potentially applicable to studies of 
dynamic changes in the wetland environment. 

Two sources of stress induced by man cause severe 
detrimental effects to the environment and have 
received a great deal of public attention. Air pollution 
and acid rain can severely stress many species of both 
natural and cultivated vegetation, although the extent 
and spatial distribution of the damage is essentially 
unknown. Remote sensing data offers the possibility 
of detecting, quantifying, and mapping various 
manifestations of stress within extended vegetation 
canopies. TM observations, in particular, can be used 
to detect variations in canopy closure, chlorophyll 
abundance, leaf moisture, and canopy temperature 
which may be diagnostically related to particular 
forms of vegetation stress. In principle, it may be 
possible to develop spectral indices for different 
ecosystems that are diagnostic of stress conditions 
caused by the deposition of airborne pollutants. The 
occurrence of such stress conditions could, in turn, 
be rel.ated to local variations in slope, slope azimuth, 
elevation, and annual precipitation to learn more 
about the processes governing the transportation and 
deposition of airborne pollutants. 


Fundamental Biospheric Processes 

The living organisms that inhabit the outer portions 
of the Earth play an integral role in many planetary 
scale processes that occur over a wide range of time 
scales. For example, the Earth’s biota are involved in 
the chemical modification and physical destruction 
of geological materials exposed at the Earth’s surface. 
They also provide a means of collecting, storing, and 
transforming solar epergy which can have a long term 
impact upon the global climate. The potential applica- 
tion of TM observations to the study of these fun- 
damental biospheric processes is illustrated in the 
following discussion of selected biospheric 
phenomena. 

Biogeochemical Cycling. The biogeochemical 
cycling of carbon, nitrogen, sulfur and phosphorous 
is critical to all life on Earth, These elements are 
essential to the metabolism of all forms of living 
biomass, and they must be available in a wide variety 
of forms and quantities to support specific terrestrial 
llfeforms. Vegetation, both above-ground and below- 
ground phytomass, plays a major role in the 
biogeochemical cycles of all the elements listed above. 
Since vegetation can be drastically altered by fire, 
flood, and storms, the remote sensing of 
environments subjected to such catastrophic events 
can provide extremely useful and otherwise difficult 
to measure information about the storage and release 
of selected elements. TM data may improve the 
measurement accuracy of surface conditions, such as 
vegetation, that are significant input variables for 
models of elemental fluxes over large geographical 
areas. 

Energy Transfer at the Earth’s Surface. The 
distribution of soils, vegetation, and water at the 
Earth’s surface controls the amount of solar radiation 
that is absorbed by the Earth, Vegetation density and 
distribution have a major impact upon the local 
energy budget of the Earth’s surface because of the 
role of plant transpiration in latent heat transfer. Plant 
transpiration is a far more effective process of 
liberating water vapor from the Earth’s surface than 
direct evaporation from a free water surface. 
Consequently, vegetation can have a major influence 
upon the net heat flux of extended areas. Changes in 
the form and distribution of surface vegetation can 
greatly alter the energy budgets of both localized areas 
and extended regions. TM observations can be used 
to monitor .spatial and temporal variations in surface 
vegetation. These types of observations could be used 
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to place boundary conditions on models of net energy 
flux at the Earth’s surface. In addition, the thermal 
band on the TM (band 6) can be used to monitor 
variations in surface thermal emission. What is truly 
unique about the TM is its abiiity to monitor variations 
in surface cover conditions and thermal characteristics 
at a synoptic scale. The use of TM data to place 
boundary conditions on models of energy flux over 
extended areas could potentially lead to major 
breakthroughs in studies of regional meteorology and 
atmospheric circulation, particularly in environments 
that experience large variations in surface cover 
conditions on a seasonal basis. 

GEOLOGY 

Introduction 

Geology is the scientific discipline devoted to the 
study of the rocks and soils that make up the Earth. 
These materials are the ultimate sources of the gases 
that are currently found in the Earth’s atmosphere, 
the water that is found in the Earth’s hydrosphere, 
and the nutrients that support the vegetation found 
on the Earth’s surface. The geological materials that 
constitute the outer crust of the Earth play an impor- 
tant role in the biogeochemical cycling of key 
elements and organic compounds that are essential 
for life on our planet. In addition, geological resources 
found in the Earth’s outer crust are essential for the 
daily operation and future growth of our modern in- 
dustrial society. An understanding of the geological 
history of the Earth and the processes that shape its 
crust is fundamentai to being able to predict the ^ture 
behavior of our planet, and to forecast the quality of 
life for future generations of mankind, 

Geologic processes occur over a wide range of 
spatial and temporal scales, A considerable degree of 
insight into the nature of these processes can be ob- 
tained by simply observing the current composition 
and distribution of geological materials that are 
situated near the Earth’s surface. These types of obser- 
vations place important boundary conditions on 
models of crustal growth and evolution. In addition, 
they enable geologists to unravel the geologic history 
of specific crustal provinces. 

During the past forty years, geologists have 
developed considerable expertise in deriving 
geological information from Imagery of the Earth’s 
surface. Following World War II, geologists were 
quick to adopt aerial photographic techniques that 


had originally been developed for purposes of military 
reconnaissance. These techniques are now used 
routinely to detect variations in the color, texture, and 
weathering characteristics of surficial geological 
materials. Information of this nature provides impor- 
tant clues about the chemical composition and 
physical properties of surficial deposits. In vegetated 
areas, geologists employ aerial photography to detect 
variations in vegetation density, distribution, or con- 
dition that may be correlated with the composition 
or properties of underlying geological materials. Fi- 
nally, geologists employ aerial photographs to detect 
variations in surface relief that may be Indicative of 
folds and faults in the Earth’s outer crust. 

Prior to the nineteen-seventies, optical Imagery of 
the Earth’s surface was collected and analyzed 
primarily in the form of photographic prints and 
transparencies (l.e., in analog format). Two important 
events occurred during the nineteen-seventies that 
served to expand the range of photointerpretive 
techniques potentially available to geologists. One of 
these events was the advent of the Landsat series of 
spacecraft. Landsat was the first spacecraft designed 
to acquire multispectral Earth Imagery at a synoptic 
scale on a repetitive, global basis. This imagery was 
collected in digital form and stored on magnetic tapes. 
The second major technological advance that oc- 
curred during the nineteen-seventies was the pro- 
liferation of computer facilities that could be used to 
reduce and analyze data in digital format. By the end 
of the seventies, these facilities were available to 
many, if not most, geologists. 

The availability of computerized facilities for digital 
data analysis and the availability of Landsat imagery 
in digital format has lead to the development of new 
techniques for extracting geological information from 
pictures of the Earth’s surface. Following the initial 
launch of Landsat 1, Landsat data products were 
analyzed primarily in analog format, employing the 
conventional methods of photointerpretation that had 
been developed in the past. During the past ten years, 
however, there has been a growing realization among 
geologists that Landsat imagery can be digitally 
manipulated to enhance the expression of specific sur- 
face features on a highly selective basis. This realiza- 
tion has lead to the development of a wide range of 
digital image analysis techniques which complement 
and extend conventional methods of photointerpreta- 
tion. With growing sophistication in digital image 
analysis, geologists have come to treat Landsat im- 
agery as a geophysical data set, comparable to seismic 
data or potential field surveys, which must be digitally 



manipulated in various fashions to address different 
types of geological questions. 

Our current understanding of the utility of Land- 
sat imagery for geological research and applications 
is based primarily upon the analysis of imagery ac- 
quired by the Multispectral Scanner (MSS) on Land- 
sats 1,2, and 3. The Thematic Mapper (TM) on Land- 
sats 4 and 5 is a second generation sensor system 
which is far superior to the MSS. As described 
elsewhere in this document, the TM possesses ap- 
proximately twice the spectral resolution, two and 
one-half times the spatial resolution, and three times 
the radiometric sensitivity of the MSS. The TM offers 
geologists a new range of Earth observation 
capabilities which can potentially be applied to basic 
studies of crustal processes and to the exploration of 
remote crustal provinces. The purpose of the follow- 
ing discussion is to briefly summarize the anticipated 
observational capabilities of the TM in a geological 
context, and to identify prospective research topics 
that could potentially be addressed through the 
analysis of TM data. The overall goal of this discus- 
sion is to stimulate creative thinking about how TM 
imagery could be used in the future to conduct basic 
studies of geological conditions and processes on the 
Earth, 

Geological Observation Capabilities 
of the Thematic Mapper 

In attempting to exploit the measurement 
capabilities of the TM for basic research in geology, 
two questions immediately come to mind, namely; 

1) What types of geological information can poten- 
tially be derived from TM measurements?, and 

2) What is unique about the TM’s Earth observa- 
tion capabilities? 

A variety of experiments have been conducted in 
the past employing airborne scanner systems which 
simulate the measurement capabilities of the TM. In 
addition, geologists have gained some Insight into the 
characteristics of orbital TM data through the analysis 
of Landsat 4 TM imagery that was collected during 
the fall and winter of 1982 over portions of the 
western U. S. These earlier studies have indicated that 
multispectral TM imagery can potentially be used to 
obtain information about the following characteristics 
of the Earth’s surface: 

Rock and Soil Lithology. Iron-bearing minerals 
possess diagnostic spectral absorption features within 
the visible and near infrared portions of the 


electromagnetic spectrum. These features extend 
acro.ss TM bands 1-4 and they can be used to detect 
limonite stained rocks and soils in arid environments. 
Hydroxyl(OH) bearing minerals such as clays and 
carbonate(C 05 )-bearing minerals such as calcite are 
characterized by diagnostic absorption features in the 
so-called shortwave infrared spectral region. TM band 
7 is situated in the shortwave infrared and it is 
sensitive to the presence of clays and carbonate rocks 
where they arc exposed directly at the Earth’s surface. 
TM bands 2-5 and 7 can be used in combination to 
detect iron-bearing clay minerals which are 
commonly considered to be diagnostic indicators of 
hydrothermal alteration. Under certain circumstances, 
it is even possible to detect variations in the local 
concentration of surficial clay minerals (see Figure 
III-6), TM band 6 is designed to measure the intensity 
of thermal radiation being emitted by the Earth's 
surface. Previous airborne experiments have shown 
that night-time thermal infrared surveys can be used 
to discriminate certain types of rocks on the basis of 
their thermal inertia or heat retention properties. 
Quartz, in particular, has a relatively high thermal 
inertia, Consequently, quartz-bearing rocks such as 
sandstone are commonly warmer than their 
surroundings in night-time thermal infrared imagery 
and they can be readily discriminated. The various 
types of compositional information that can 
potentially be obtained through the analysis of TM 
imagery are not sufficient to characterize the lithology 
of many common rocks and soils. However, analysis 
of airborne TM simulator data has demonstrated that 
the stoichiometry and abundance of the minerals 
described above are sufficiently variable in nature to 
produce marked differences in the spectral reflectance 
properties of natural terrain. In particular, analysis of 
TM simulator imagery acquired in semi-arid regions 
has shown that variations in terrain reflectivity 
detected in the TM spectral bands can be used to 
discriminate a wide variety of natural rock and soil 
units exposed at the Earth’s surface. A high degree 
of lithologic discrimination can commonly be 
achieved in TM image analysis, even in cases where 
direct identification of individual lithologic units is 
difficult or impossible. 

Surface Landforms. The 30 meter instantaneous 
field of view of the TM is ideally suited for large scale 
studies of surface geomorphology. The TM is a nadir- 
viewing sensor, and TM imagery can be used to detect 
the presence of specific landforms and characterize 
their shape in a planimetric sense. Standard 
photoclinometric techniques can be applied to TM 
imagery obtained under different solar illumination 
conditions to determine the relative relief of selected 
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The two imafies shown above have been ciern eii from 
a /.andsal 4 TM imaf(e of the Stiver Hell, Arizona area 
that was ortfitnaliy acquired in Sovemher, I9HJ F.ach 
suhscene is approximately 15x15 kilometers in size 
The suhscene on the left is a TM hand 4 view of the 
area, the suhscene on the rif^ht is a ratio irnaye con 
structed from TM hands 5 and The mountainous 
area in the center of the suhscene has been mined in 
the fHist for copper The hriffht rhomhoidal Jeatures 
seen in the TM hand 4 suhscene are tailinffs ponds 


TM BAND 5 TM BAND 7 

NASA HO EtB4 1I77I3I 
7 76 14 

assiK iated with these earlier mininff actii dies I 'ana 
lions in the TM 5/TM " ratio ima^e are produced hy 
variations in the surface concentration of clay 
minerals Proffressively larger values of the TM 5n'M 
" ratio are indicative of prof^ressively larf^er concen 
trations of surface clays The hrif^htest areas observed 
in the ratio ima^e represent the surface expression oj 
the hydrothermal alteration zone that was produced 
when copper mineralization was orif^inally emplaced 
in this area (Illustration courtesy of Pat Chevez. T 
S (leolofiical Survey) 


bndforms Finally, TM imagery can be used to 
determine the distribution of landforms over broad 
arexs and detect alignments, offsets, and variations in 
landform shape which may be indicative of local 
folding or faulting in the crust Certain types of rtK'ks 
have a natural tendency to weather in characteristic 
fashions in different types of weathering 
environments The improved spatial rest)lution of the 
TM will enable geologists to detect more subtle 
variations in the weathering habit of rock units that 
outcrop directly at the Earth's surface. Under certain 
circumstances, landform morphology can be used to 
infer the lithology of surficial rock units. It is also 
possible to infer the strike and attitude of surficial rtx k 
units on the basis of landform shape and distribution 
under certain circumstances. Studies conducted to 
date with Landsat 4 TM data indicate that TM imagery 


will support terrain analysis and structural mapping 
at a scale of 1 24, ()()() 

Ground Temperature. The TM thermal band 
(band 6) can be used to survey areal variations in 
ground temperature at a spatial resolution of 120 
meters In principle, this capability could be used to 
study variations in surface heat tlow related to 
regional volcanism of both an extrusive and/or 
intrusive nature Night time thermal surveys have 
been used in the past to detect variations in rock 
lithology xs discussed above. Night time thermal 
surveys have also been employed to detect surface 
faults and fracture patterns. Near surface moisture 
commonly collects along fractures and joint surfaces 
Under these circumstances, near surface fractures are 
significantly cooler than their local surroundings and 
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they are readily detectable In thermal infrared 
imagery. 

Surface Vegetation. MSS Imagery has been used 
in the past to detect the occurrence and map the 
relative abundance of surface vegetation over a variety 
of temporal and spatial scales. These same types of 
studies can now be performed at an improved level 
of spatial resolution with the TM. In addition, the 
Improved spectral resolution of the TM can be used 
to differentiate a wider variety of vegetative species. 
Specifically, the TM spectral bands can be used to 
discriminate both deciduous and coniferous tree 
stands that are dominated by a single species such as 
white oak, red maple, loblolly pine, and white pine. 
Past research has shown that vegetation distribution 
and density can be correlated with discontinuities in 
the lithologic characteristics of underlying geological 
materials under certain circumstances. Past research 
has also indicated that the TM spectral bands can be 
used to detect stress conditions in natural vegetation 
which are related to soil toxicity and water 
deprivation, Vegetation stress can sometimes be 
related to anomalous geological conditions such as the 
concentration of metallic elements in surface soils or 
the local seepage of underground hydrocarbons. 

The potential utility of orbital TM surveys for basic 
geological research is also governed to a large degree 
by the spatial extent and temporal frequency of TM 
observations, Multispectral TM data are unique in 
three important respects: they are collected over a 
ground swath of 185 kilometers on a near global basis 
every 16 days. Each of these considerations has 
important implications for the types of geological 
investigations that can potentially be conducted with 
TM imagery, and the experimental methodology to 
be employed in addressing specific types of geological 
problems. These implications are briefly summarized 
below. 

Synoptic Scale Observations. The synoptic scale 
of TM imagery can provide a basis for interrelating 
various types of geological Information obtained at 
a local scale by conventional field mapping and 
sampling techniques, For example, TM imagery can 
be used to integrate maps of surflcial geology, fracture 
patterns, or soil geochemistry obtained at much more 
limited spatial scales by relating the occurrence of 
specific surface materials or geologic structures to 
spectral patterns and landform patterns observed in 
a TM image, Furthermore, the synoptic scale of TM 
imagery may be used to identify geological features 
of regional significance that are not readily apparent 
to ground based observers or aerial photointerpreters. 
For example, TM imagery can be used to recognize 


major systems of crustal fractures which may have a 
discontinuous expression at the Earth’s surface. TM 
imagery can also be used to recognize key marker 
horizons in sedimentary terrancs which may outcrop 
in a discontinuous fashion over a large area. The 
synoptic scale of TM imagery is inherently useful for 
Integrating and extending point measurements of 
geological conditions that are spatially discontinuous 
in nature. 

Near Global Observations. In principle, TM data 
can be collected over more than 90% of the exposed 
landmass area of the Earth. The availability of 
multispectral Earth imagery acquired at common 
scales and formats over most of the world poses some 
interesting opportunities for basic crustal research. 
For example, many remote parts of the world remain 
unmapped at the present time due to their 
inaccessabllity. TM imagery could be used to make 
pioneering observations of geological conditions in 
many arid, tropical, and arctic regions. Analysis and 
interpretation of TM imagery acquired over such areas 
could lead to the production of regional scale 
geological maps for the first time. The near global 
availability of TM imagery also presents an 
opportunity to test geological models or theories that 
are based upon intensive studies of specific crustal 
features. For example, our current understanding of 
continental rifting is based upon detailed studies of 
a limited number of continental rifts (e.g, the East 
African Rift, the Rio Grande Rift, and the Red Sea Rift). 
TM imagery could be used, in principle, to determine 
similarities and differences between those rift features 
that have been studied intensively in the past as a 
means of developing more generic models of rift 
formation and evolution. In addition, TM imagery 
could be used to study continental rifts which have 
received little attention in the past as a means of 
testing and refining existing models of rift formation. 
This is not to suggest that crustal processes can be 
modeled solely on the basis of TM image analysis. 
However, the TM is a powerful tool for conducting 
comparative studies of crustal features on a 
worldwide basis. Studies of this nature can make a 
significant contribution to current understanding of 
crustal processes, particularly when they are 
performed in conjunction with more conventional 
methods of crustal analysis. 

Repetitive Multitemporal Observations. There 
is a popular misconception that geologists have little 
interest in multitemporal Earth Imagery. It is true that 
many geologic phenomena occur over time scales that 
are several orders of magnitude greater than the 16 
day repeat cycle of the Landsat 4 and 5 spacecraft. 
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However, certain geologic processes such as lava flow 
emplacement, aeolian transport of desert sand, and 
surface flooding and sedimentation can be studied 
vlth multitemporal TM Imagery. Multitemporal Farth 
imagery is also useful for structural mapping because 
certain geological structures may be more 
prominently displayed under different solar illumina> 
tlon conditions. Earlier studies have demonstrated 
that selected structural features may be more readily 
apparent in Landsat Imagery acquired in late autumn 
(following crop harvesting), early winter (following 
periods of snowfall at higher elevations), and/or late 
winter (following the early melting of regional 
snowpacks at lower elevations). Consequently, 
photogeologists conducting regional scale studies of 
crustal structure commonly employ Landsat imagery 
acquired during different seasons. Finally, 
multitemporal TM observations will play a major role 
in future studies of geobotanical phenomena. Changes 
In the spectral reflectance characteristics of vegetated 
terrain during the annual growing season are 
potentially an important source of Information about 
the nature of the vegetation that occurs on the surface. 
In addition, variations in the duration of the growing 
season of selected species can potentially provide 
information about the lithology of underlying 
geological materials. 


Prospective Research Themes 

The purpose of the following discussion is to Iden- 
tify prospective research themes in the geosciences 
that could be addressed In new and innovative ways 
through the analysis of TM imagery. These themes in- 
volve the use of TM data for the exploration of remote 
areas, the study of crustal processes, and research in 
the interdisciplinary field of geobotany. A series of 
topical problems is presented to illustrate the nature 
of these themes and provide the reader with examples 
of research projects that are potentially amenable to 
TM analysis. The research themes and topical prob- 
lems discussed below are not Intended to be an ex- 
haustive or comprehensive listing of all geoscience 
investigations that could potentially be conducted 
with TM imagery. Furthermore, the various pro- 
cedures for TM data collection and analysis described 
below are not intended to be definitive blueprints for 
the conduct of specific investigations. The overall goal 
of the following discussion is to stimulate creative 
thinking about different types of geological problems 
that could be addressed through the analysis of TM 
data. This discussion is not intended to be restrictive 
or prescriptive in any way. 


Earth Exploration 

Current understanding of the surface structure and 
composition of the Earth's crust Is highly variable. 
Certain crustal provinces such as the Basin and Range 
province in the southwestern U.S. have been studied 
extensively in the past. Knowledge of the composi- 
tion and distribution of the rocks exposed in tills 
province has been accumulated through years of field 
mapping and sample analysis. This knowledge is sum- 
marized in geological maps which display the local 
distribution of rock types, the strike and attitude of 
rodk formations, and major structural features such 
as folds and faults. 

Many non-geologists consider geological maps to 
be definitive statements of geological conditions 
within a specific area. In reality, geological maps are 
"living documents" which must be changed as new 
concepts are developed concerning the nature and 
history of local rock units, the structure and composi- 
tion of adjacent districts, and the nature of the 
geologic processes that acted upon the Earth’s crust 
in the past. Geologic maps generally summarize the 
results of earlier scientific investigations, and, like any 
good scientific study, they tend to raise many more 
questions about the nature and history of a particular 
province than they necessarily ansv/er. To many 
geologists, a geologic map does not represent an end 
in itself, but rather it is a point of departure for many 
other types of geologic investigations. As discussed 
in the previous section, analysis of multispectral TM 
imagery can yield information about rock lithology 
and crustal structure at synoptic scales. This type of 
information can play an important role in develop- 
ing regional scale geologic maps, In principle, this 
type of information can be used to conduct pioneer- 
ing studies of geological conditions in both well 
studied and poorly understood crustal provinces. 
Three general categories of prospective TM mapping 
projects are briefly outlined below; 

(1) Current geological understanding of many parts 
of the world is quite rudimentary or lacking 
altogether. This is particularly true for crustal pro- 
vinces located in arid, tropical, and arctic regions 
which are remote and relatively inaccessible, The 
global availability of TM imagery provides a unique 
opportunity to conduct observations of regional 
geologic conditions and establish reconnaissance 
geological maps of such areas for the first time. By 
way of example, TM Imagery could potentially be 
used to identify local areas of rock outcrop in arid 
terrain, correlate outcrop areas of similar lithology, 
and develop preliminary maps of bedrock geology, 
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Tliis type of approach is particularly applicable to 
arid environments in wliich outcrop areas are 
discontinuous and widely separated by deposits of 
windblown sand. Alternatively, the improved 
spatial resolution of TM imagery could be used to 
conduct detailed mapping studies of drainage pat* 
terns in tropical terrain. Under certain cir- 
cumstances, these drainage patterns can be used to 
detect lithologic boundaries between different 
types of bedrock formations and to identify major 
structural discontinuities within the underlying 
crust. 

(2) TM image analysis can also refine or extend cur- 
rent geologic understanding of areas that have been 
studied extensively In the past with conventional 
field mapping techniques. As discussed earlier, 
multispectral TM surveys can provide Information 
about Iron occurrence, clay mineralogy and abun- 
dance, and the carbonate mineralogy of surface 
rocks on a synoptic scale. These various types of 
lithologic information are not always readily ob- 
tainable with conventional field mapping tech- 
niques. Such information could be used to extend 
existing geologic maps of selected areas. For exam- 
ple, the .sensitivity of TM band 7 to variations in 
carbonate mineralogy could potentially be used to 
map facies variations in carbonate formations that 
are commonly exposed in the intermontane basins 
of the western U.S. Alternatively, the sensitivity of 
TM bands 2-5 and 7 to clay mineralogy and abun- 
dance could potentially be used to map areas of 
hydrothermal alteration in arid and semi-arid 
environments, The spectral coverage of the TM is 
particularly well suited for the discrimination of 
Iron-bearing clay minerals. This capability is directly 
applicable to studies of ancient weathered soils such 
as those found In the Interior of Australia. 

(3) Finally, the synoptic scale and spectral coverage 
of TM imagery could potentially be used to tie 
together conventional geological maps prepared at 
more limited spatial scales, and to develop regional 
scale maps of selected crustal provinces for the first 
time. The synoptic dimensions of individual TM 
scenes can easily support mapping studies con- 
ducted at scales of 1 .'250,000. TM image mosaics 
could potentially be used to develop maps of 
selected provinces at even larger scales. In princi- 
ple, it should be possible to develop spectral and 
spatial indices for correlating rock units of similar 
lithology that have been independently mapped by 
different field workers in the past. These Indices 
could potentially be used to interrelate the occur- 
rence of specific rock units appearing on existing 


geological map.s\ These indices could also be used 
to interpolate between existing maps and identify 
the probable occurrence of similar rock units in ad- 
joining areas whlcli have not been mapped In the 
past. The use of TM multispectral Imagery to inter- 
relate existing geologic maps and interpolate the oc- 
currence of specific rock units could lead to the 
development of regional scale maps of selected 
crustal provinces for the first time. Maps of this 
nature could, in turn, be used to develop alternative 
models for the history and evolution of such 
provinces. 

The common feature of all of the mapping In- 
vestigations described above Is the use of TM data to 
conduct exploratory observations of geological con- 
ditions in different crustal regions. They .specifically 
make use of the TM's measurement capabilities to 
establish new and/or improved understanding of 
crustal composition and structure. As such, they 
represent the use of the TM system to explore the sur- 
face of our planet in a pioneering fashion. 


Geologic Processes 

The current di.strlbution of crustal rocks and land- 
forms is the ultimate result of a myriad of proce.sses 
that have shaped the crust in the past. In many other 
Earth science disciplines, such as oceanography or 
meteorology, it is possible to conduct repetitive 
observations of in situ conditions and deduce 
something about the nature of the processes that are 
responsible for observed changes in surface condi- 
tions. Geology differs from many of these other 
disciplines because so many geological processes 
occur over time scales that are much longer than the 
average human lifetime. In fact, many geological pro- 
cesses occur over periods of time that arc longer than 
recorded human history! Geologists examine the 
nature, sequence, and configuration of rocks exposed 
in different parts of the world, and they employ in- 
ductive reasoning to develop models of many crustal 
processes. The TM can contribute to this inductive 
method of studying the Earth by providing geologists 
an opportunity to examine the current configuration 
of crustal materials and landforms almost anywhere 
in the world, In fact, TM observations provide 
geologists with a unique opportunity to study a 
variety of crustal provinces in a comparative fashion 
that has not been possible in the past. 

At the same time, however, events such as the erup- 
tion of Mt. St. Helens in 1980 or the San Francisco 
earthquake of 1906 serve as compelling reminders 
that certain geologic phenomena do occur over time 
scales that are directly amenable to human observa- 


Uon. The TM can obviously play a role in the stuUy 
of these short term phenomena, particularly when 
they occur in remote anti relatively unpopulated 
areas. 

The purpose of the following discussion is to il- 
lustrate some of the types of geological processes that 
could potentially be studied through the analysis of 
TM data- The pltenomcna described below do not 
represent the only processes that could potentially be 
studied, nor do they necessarily represent the most 
important topics for process-related research. 

Erosloaat Processes, The current morphology of 
the Earth’s crust reflects a balance between the effects 
of constructional phenomena such as tcctonism and 
•olcanism, and destructive phenomena such as fluvial 
and aeolian erosion. The nature of the balance 
between these competitive phenomena can vary 
enormously from one crustal province to the next. 

Landscape morphology is locally governed by 
regional environmental conditions and the nature of 
the rocks that arc locally exposed at the Earth’s 
surface. Conventional geological wisdom predicts that 
certain landform assemblages are always produced 
wlien selected types of rocks are exposed to specific 
types of weathering environments. For example, 
limestone formations exposed to tropical weathering 
conditions break down in a characteristic fashion to 
produce a highly dissected landscape known as karst 
terrain (see Figure II1-7). The extent to which 
variations in limestone lithology or interbedded 
formations of different lithology affect the nature and 
development of karst terrain is poorly understood at 
the present time. Furthermore, qualitative understand- 
ing of the equilibrium landform assemblages that 
should result from certain combinations of rock 
lithology and environmental conditions is lacking 
altogether. 

TM observations could potentially be used to 
extend current understanding of how erosional 
processes control landscape development in both a 
quantitative and qualitative sense. In the case of well- 
known associations betv/een landscape morphology, 
rock lithology, and weathering conditions, TM 
observations of landform shape and distribution could 
be used to establish quantitative models of the 
landrorm assemblages that should be produced in a 
specific area. Detailed measurements of landform 
morphology could potentially be related to the 
surface exposure ages of specific formations to 
establish quantitative models of landscape evolution 
as a function of time. In theory, these types of 
relationships could potentially be inverted to infer the 
length of time a particular area has been subjected to 
a specific set of weathering conditions. 


Improved understanding of the landforms that 
should appear in a particular province would also 
enable geologists to identify relict landforms with 
greater accuracy. Relict landforms are ancient features 
which have been Inherited by the regional landscape 
and they can provide clues about regional weather- 
ing conditions that may have existed in the past. The 
shape and arrangement of relict landforms can. in 
principle, be used to reconstruct the paleociimaiic 
history of selected portions of the Earth. Information 
of this nature can, in turn, place Important constraints 
on climatic models. 

TM imagery could also be employed to conduct 
pioneering observations of landform morpliology and 
distribution in areas characterised by complex 
bedrock geology or transitional weathering condi- 
tions. These ob.servations would enable geologists to 
develop qualitative models of equilibrium landform 
assemblages for certain rock lltliologies and weather- 
ing environments for the first time. In fact, research 
of this nature would provide an interesting test of the 
theory that a specific assemblage of landforms is 
necessarily produced when certain types of rocks arc 
exposed to a particular setof weathering conditions. 
Under certain geological and environmental condi- 
tions, landscape evolution may be more probabilistic 
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equilibrium assemblage of land forms may be pro- 
duced in transitional weathering environments of 
similar geology. If true, this conclusion would have 
important implications for the geologic and climatic 
histories of selected areas. 


Tcctonism and Crustal Deformation. The solid 
outer shell of the Earth is known as the lithosphere. 
Crustal rocks constitute the outermost portion of the 
lithosphere and they are underlain by denser materials 
that form part of the Earth’s mantle. The Barth’s 
lithosphere is segmented into a collection of 
fragments that are referred to as “plates". These plate, s 
can move in both a horizontal and vertical sense with 
respect to one another, Plate collllsional processes are 
responsible for the construction of mountain chains, 
the formation of oceanic trenches, and much of the 
earthquake activity that occurs today. 

A wide variety of internal and external forces act 
upon the Earth’s lithosphere plates. These forces 
include the internal upwelling of molten material 
within the Earth’s lower mantle, the centrifugal 
acceleration produced by the rotation of the Earth, 
gravitational accelerations produced by the Earth’s 
core and lower mantle, and tidal forces exerted by 
the moon and the sun. Variations in the magnitude 
and direction of these and other forces that occur 
across individual plates produce stresses within the 
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I jrih I ruM which are |Nrri«Kik ully rrhrvcJ ihrough 
c|)is«k1u. crusul motions Thr dctornuiion ol the 
I jrih's = rust in rcsfnmsc to these ton es is expressed 
in the form ot folds and faults 

In unraveling the history of individual provincex, 
geologists seek to determine the orientation and 


magnitude ol the stress Held that is currently acting 
u|>on tlKT «. rust, and i< > determine how this stress field 
lias changed t»ver die course ol geological lime The 
orientation, extent, and sense of n'.oiion ol crustal 
folds and fault'' plat e im|Kiriani boundary conditions 
ufkin the geometry and magnitude of stresses that 



Ol vr (iistam.es oj HO kilometers sufft(estinfi that tin- 
(iissei tion oJ the renional limestone formations has 
(K'curred aloiifi pre-existinff fracture systems The 
im/iroreti spatial resolution of the TM could he used 
in thejuture to study the morphometric characteristics 
of karst landforms at a level of spatial detail that has 
not been /nissihle in the past Such studies could 
potentially he conducted in a f(lohally comparative 
fashion to determine similarities and differences in the 
morplHilofiy of different karst terrains (Imafie courtesy 
of Sick Short, A'd.Vd — floddard Space Flight <’ enter) 




Figure 111-7 

Multispectral Scanner hand 7 imaffe of karst terrain 
in Guanffxi Province, China This tmane was 
onumally acquired in November, /97J by l.andsat I 
It is a standard MSS image encompassing an area of 
IHS X IH5 kilometers The darker, highly dissected 
limestone formations exposed in this region contrast 
sharply with surrounding formations in terms of their 
tone and texture Steep-sided linear depressions 
exfxtsed in the karst terrain can he traced continuously 
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have deformed the crust In the past. The improved 
spatial resolution of TM Imagery can clearly con* 
tribute to the identification of such features, and to 
the determination of displacement along surface faults 
(see Figure 111*8). Furthermore, the synoptic scale of 
TM Itnagery enables a photogeologist to determine 
the length, alignment, and relative displacement of 
such features on a regional scale. Observations of this 
nature are extremely useful in attempting to 
distinguish regional trends from localized pattern ' of 
apparent crustal deformation. They also enable 
geologists to identify key districts within a broader 
province where ground based measurements of fault 
motion and displacement are needed to resolve alter- 
native interpretations of stress field geometry. 

Field investigations of fault displacement play a key 
role in determining the relative ages of different fault 
and fracture systems, In many instances, individual 
rock formations are multiply deformed by a series of 
tectonic events. Detailed observations of the sense of 
motion and apparent displacement that has occurred 
along crosscutting fracture systems enable geologists 
to assign relative ages to specific deformation events. 
Information of this nature is essential in attempting 
to determine how the regional stress field has chang- 
ed over the course of geologic time Synoptic scale 
TM imagery could prove to be extremely useful in 
Identifying areas of fault and fracture intersection, and 
targeting such areas for detailed field analysis. The 
distribution and orientation of surface folds and faults 
can potentially provide Important insight Into the 
structure of the deep crust. Crustal provinces are com- 
monly defined on the basis of the age or lithology of 
the rocks exposed at the Earth’s surlace. Crustal prov- 
inces can also be defined on the basis of the tectonic 
fabric or style of deformation of extended regions. 
.''iS described above, TM imagery can be used to map 
patterns of crustal deformation over a wide range of 
spatial scales, Such patterns can, in turn, be used to 
identify crustal provinces characterized by different 
styles of structural deformation. Boundaries 
separating structural provinces may not necessarily 
correspond to provincial boundaries identified on the 
basis of rock lltholo^ • or stratigraphy. Sharp discon- 
tinuities in the tectonic fabric of individual regions 
may be indicative of major structural discontinuities 
in the lower crust. In principle, TM imagery could be 
analyzed in conjunction with more conventional 
geophysical data (e.g. .seismic surveys and potential 
field measurements) to detect and characterize major 
structural features in the lower crust. 

Crustal stresses are commonly relieved at a local 
scale along pre-exioting zones of weakness. These 
zones of weakness may be aligned in a highly random 
fashion with respect to the regional stress field. 


Seismic data can provide information about the loca- 
tion, depth, sense of motion, and strain energy re- 
leased during an earthquake event. Information of this 
nature can be used to estimate the direction and 
magnitude of regional crustal stresses If the orienta- 
tion and attitude of the fault zone can be determin- 
ed. Gru.stal faults are commonly exposed at the Earth's 
surface. Under certain circumstances, the orientation 
and attitude of such features can be determined 
through th. analysis of TM imagery If historical 
seismic activity can be related to specific zones of 
weakness identified in TM Imagery, it should be possi- 
ble to develop stress field models for selected por- 
tions of the Earth's crust. Combined analysis of TM 
imagery and historical seismic data could prove to be 
particularly effective for determining the nature of the 
current stress field in many remote and poorly ex- 
plored crustal provinces. 

Volcanism. Volcanism is, by its very nature, an 
episodic process that adds material to the Earth’s 
crust. Like tectonism, it has had a dominating 
influence upon the recent geologic history of certain 
crustal provinces. Many areas of recent volcanic 
activity are located in remote parts of the world, Basic 
understanding of the nature and frequency of volcanic 
activity in many of these p;ovinces is lacking 
altogether. Furthermore, our basic understanding of 
how volcanoes work and what factors govern the 
outbreak of volcanic activity is quite rudimentary at 
the present time. 

Knowledge of the composition, age, and volume 
of lava that has been erupted by an individual volcano 
in the past places important boundary conditions on 
models of future eruptive activity. The TM is 
particularly well suited for mapping the spatial extent 
of lava flows in light of the 30 meter spatial resolution 
of TM imagery. The TM spectral bands are also well 
suited for detecting variations in the iron oxide 
chemistry and clay mineralogy of weathered flow 
surfaces. Under certain circumstances, it should be 
possible to distinguish flows of different chemical 
composition on the basis of the spectral characteristics 
of their weathered surfaces. Past research has also 
suggested that the devitrification of lava flow surfaces 
and the corresponding development of clay alteration 
products is directly related to the surface exposure 
age of individual flows. As discussed earlier, variations 
in the surface abundance of clay minerals produce 
variations in natural surface reflectivity which can be 
measured by the TM (specifically by TM band 7). In 
principle, it may be possible to exploit this 
observational capability of the TM to detect variations 
in the clay abundance, and, hence, the relative ages 
of weathered flow units (see, for example, Figure 
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Fifiure Ul-ti 

Sharp diacontmutties in the shape and dtstnbufion oj 
surface i'lnd forms can be used to identify the surface 
expression of crustal faults, as illustrated in this TM 
imafte of southern California. This image was 
originally acquired in November, 1982 over path 
40/row ■? 5 of the Landsat 4 and 5 World Reference 
System The arrow denotes the surface expression of 
a right lateral shear fault that is regionally exposed 
in this area the linear feature observed in the image 


is the eastern most extension of the tiarlock Fault 
which is a major tectonic feature in southern Califor- 
nia The section of the fault displayed in this imageis 
roughly 35 kilometers long Note the marked dif- 
ferences in alluvial fan development along various 
segments of the fault The improved spatial resolution 
of the TNi should enable geologists to detect much more 
subtle expressions of neotectonic features in many dif- 
ferent parts of the u>orld, particularly in semi-arid en- 
vironments (Image courtesy of L(Kke Stuart, NASA — 
Coddard Space Flight Center). 

^i\.ulNAL PAG£ 19 
OF POOR QUALITY 


52 


ORIGINAL PAGE IS 
OF POOR QUALITY 

1119) Intormilion coiucrning ihr com|)o.si(ion, jgc, band on the I'M provides a unique op(x>rtuniiy to 

and extent i)f individual lava Hows derived from observe and monitor variations in the thermal 

synoptic scale TM imagery could contribute to the characteristics of natural terrain in areas of incipient 

development of predictive models which would or current volcanism I^rge scale surveys of surface 

enable geologists to forecast the nature, timing, and thermal emission could prove to he partKularly useful 

volume of future eruptive events in t more for forecasting the eruption probability of manv 

quantitative fashion stratovolcanoes, particularly those that are situated at 

On a global ssale, volcanis.n is the most efficient higher latitudes and elevations. Stratovokanoes are 

terrestrial preness for transferring thermal energy steep sided, conical landforms These features are 

from the Kerth’s Interior to its surfa«. e The thermal characterized by summit craters or calderxs which are 



Fifiure III-'-) 

older u euthered Jluu has a brou ftish visual ap- 
Thts IS a firouttd-based photofiruplj of two inter- fiearance and a hiffher net reflectance in the visible 
/infeed lava flows situated in the saddle area between fHtrtion uj the s/H\ truni than the wunfter flow surflue. 
Mauna Kea and Mauna l.oa on the island oj Hawaii TM measurements will enable neolontsts to detect 

Althouf(fo both flows /xissess similar chemical com/Mtsi- variations in the sfiectral characteristics oJ lava flow 

lions, they display very difjerent reflectance profK’r- surfaces tlntt can thc'orettcally be correlated with the 
ties The surf ace of the older flow on the nf^ht has been exfMtsure af^es oJ individual flows. This caftability 
chemically altered into an assemhlafie of amorphous could, in principle, be used to study the eruptiie history 
silica ffels, iron (Kxules, and clay minerals 7hc‘ younffer oJ many historically active volcanoes. (Photofiraph 
flow on the left is relatively fresh and unaltered The courtesy of John Adams, Tniversity of VCashinffton.) 
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commonly filled with water. The TM provides a 
unique capability to monitor temporal variations in 
the thermal characteristics of such summit lakes. 
Variations in the thermal state of volcanic lakes could 
be combined with seismic data and crustal 
deformation surveys to develop refined eruption 
forecasting models. 

Thermal surveys conducted by the TM could also 
prove useful in studying the physical properties and 
cooling characteristics of recent volcanic deposits. For 
example, recent studies have suggested that lava flow- 
ing in a fixed channel can remain Isothermal over large 
distances due to the thermal energy produced by 
viscous shearing within the flow’s interior. TM ther- 
mal surveys of active lava flows could provide a 
means of testing this theory in the future. TM fher- 
mal surveys could also be used to monitor the cool- 
ing characteristics of ash deposits such as those pro- 
duced during the 1980 eruption ofMt. St. Helens. In- 
formation of this nature could be used to refine and 
quantify models of welded tuff formation that have 
been developed In the past. 

Geobotany 

The occurrence, distribution, and condition of 
natural vegetation is governed by a variety of factors 
including seasonal thermal variations, local rainfall and 
relief, nutrient availability, animal predation, and the 
physical and chemical characteristics of the local 
geological substrate. Any one of these factors can 
potentially have a controlling influence on the ger- 
mination and development of individual plant species. 

The extent to which local geological conditions in- 
fluence plant succession, proliferation, and condition 
is poorly understood at the present time. Research in 
this area is currently being pursued by a variety of in- 
vestigators for very different reasons. From a botanical 
or ecological perspective, improved understanding of 
natural associations between bedrock, soils, and 
vegetation can be used to develop better models of 
the type of vegetation that is likely to develop in dif- 
ferent geological terranes. This knowledge is needed 
to isolate the effects of local geology and determine 
the influence of other factors, such as climate and 
drainage, upon vegetation potential. From a geological 
perspective, improved understanding of natural 
associations between bedrock, soils, and vegetation 
can be used to infer local geological conditions on 
the basis of observed variations in the characteristics 
of natural vegetation. Furthermore, improved 
understanding of natural plant-substrate associations 
would potentially enable geologists to recognize 
anomalous variations in substrate geology on the basis 


of stress conditions In natural vegetation canopies. Of 
key Interest to geologists is the ability to recognize 
vegetation stress that may be produced by abnormal 
concentrations of metallic ions In surficlal soils or the 
local seepage of hydrocarbon gases from a subsurface 
oil reservoir. 

The expanded spectral range of the TM can be us- 
ed to differentiate a far wider range of natural vegeta- 
tion species than was possible In the past with the 
Landsat MSS. In addition, selected TM bands can be 
used to detect variations In the abundance of 
chlorophyll and water in leafy vegetation. Informa- 
tion concerning species distribution, density, and 
condition derived from TM imagery could, in princi- 
ple, be combined with ancillary information concer- 
ning surface slope and orientation, soil geochemistry, 
and bedrock geology to develop models of natural 
bedrock-soll-vegetatlon associations in different 
ecological settings (see Figure III-IO). Models of this 
nature could initially be developed in specific test site 
areas where a wide variety of ancillary information 
concerning local substrate characteristics is readily 
available. They could subsequently be tested in other 
areas where canopy characteristics inferred from 
multispectral TM surveys would actually be used to 
predict the chemical nature and physical properties 
of the underlying substrate materials. The ability to 
infer local geological conditions on the basis of 
observed variations in the spectral properties of 
natural vegetation could potentially revolutionize 
geological mapping methods in many temperature 
and tropical portions of the world. 

Improved understanding of natural plant-substrate 
associations could also lead to the development of 
new methods for recognizing geologically-induced 
stress conditions in vegetated terrain. Some forms of 
geological stress have an obvious impact upon natural 
vegetation. For example, soils that develop on 
ultramafic rocks contain large concentrations of 
metallic ions which most forms of vegetation cannot 
tolerate. There are significant variations in plant oc- 
currence and distribution in ultramafic terranes which 
can be directly related to the metal toxicity of local 
soils. In many other instances, however, small varia- 
tions in the geological properties of the local substrate 
produce much more subtle variations in the 
characteristics of the natural vegetation canopy. Iden- 
tification of these subtle effects could be effectively 
studied with the TM, by developing improved 
understanding of what constitutes “normalcy” in 
terms of the plant-substrate associations to be ex- 
pected in different ecological settings. The definition 
of “normal” associations between vegetation, soils, 
and bedrock is an essential first step in attempting to 
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Mudy and classify anomalous variations in canop> 
characteristics that can potentially he attributed to 
ItK'al geological factors. 



Figure III- 10 

Natural relationships between vegetation, soils, and 
bedrock characteristics are exceedingly complex and 
poorly understood at the present time. The use of TM 
data to determine areal variations in vegetation den- 
sity, distribution, and condition in areas of known 
geology could make a significant contribution to the 
study of vegetation-soil-bedrock associations. 
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Introduction 

One oi the characteristics of the Harth that 
distinguishes it from all of the other planets in our 
M)lar system is the abundance of wate. in ga.seous, licj- 
uid, and solid form near the Karth's surface. 
Hydrology has lx:en defined as "the science that treats 
the waters of the Harth, their (K'currence, circulation, 
and distribution, tneir chemical and physical proper- 
ties, and their reaction with their environment, in- 
cluding their relation to living things " Hydrologists 
are fundamentally concerned with the transport and 
storage of water in its various forms within the Earth’s 
oceans, atmo.spherc, and continents. The cycling of 
water thiough these various reservoirs has been 
descril^ed by I/vovich (1979) as "the most magnifi- 
cent prtK'ess that takes place on the Earth. The interac- 
tion between water and the other components of 
nature is .so sub.stantial and its consequences so great 
that it requires the mo.st .serious attention” 

Many re.searchers have attempted to characterize 
the range of interactions that (K'cur within the Earth’s 
hydro.sphere Figures lll-l I and 111-12 illu.strate the 
various dimensions of the hydrologic cycle as defined 
by Eagle.son ( 1 970) and h»..'ton( 19.^1). It is important 
to realize that the proces.ses illustrated in these figures 
occur over a wide range of spatial scales. 





Figure lll-l ! 

An engineering view of the hydrologic cycle after 
Eagleson (1970). 
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A qualitative representation of the hydrologic cycle 
after Horton (J93lf 


Changes in the form and distribution of water have 
important consequences on the weather, climate, and 
ecology of the Earth’s surface. Water is the primary 
constituent of terrestrial vegetation. Variations in 
water availability can have a major impact upon the 
dynamics and distribution of the Earth's biomass on 
both a regional and global scale. The transformation 
of liquid water to snow and ice can also have a major 
impact upon regional meteorology and longer term 
variations in the global climate. In light of the in- 
fluence of water on these other phenomena, the study 
of hydrological processes at all scales is a very fruit- 
ful and substantial area of scientific inquiry (see, for 
example, Dooge, 1982). 

Spaceborne remote sensing systems possess unique 
Earth observation capabilities that can make a signifi- 
cant contribution to the study of hydrologic condi- 
tions and proces.ses. In particular, orbital sensors can 
be used to conduct synoptic scale surveys of surface 
conditions at a high spatial resolution. In addition, 
these sensors provide repetitive observations of sur- 
face conditions at uniform scales and formats on a 
global basis. Orbital remote sensing surveys comple- 
ment more conventional methods of ground-based 
observation that typically obtain high accuracy 
measurements of in situ conditions on a point-by- 
point basis. Expressed in another way, conventional 
field techniques can provide information about sur- 


ficial conditions with relatively high accuracy and 
high temporal frequency; whereas, orbital remote 
sensing techniques generally provide less accurate but 
highly dense observations over large areas, Orbital 
sensors possess observational capabilities that can 
contribute to improved understanding of hydrologic 
processes occurring at all scales, but they are par- 
ticularly well suited to the study of large scale 
phenomenon. Uniform and consistent observations 
of land cover conditions are generally lacking at large 
spatial scales. These types of observations must be 
performed if we are to understand the dynamic pro- 
cesses and storage capacities of different portions of 
the hydrological cycle. 

The measurement capabilities of the Landsat 
Thematic Mapper (TM) can potentially be used to 
Study selected aspects of hydrology in new and in- 
novative ways. However, to make effective use of the 
TM for such studies, it is first necessary to appreciate 
the observational capabilities and limitations of the 
TM in a hydrological context. The Landsat 4 and 5 
spacecraft which carry the TM sensor each pass over 
the same spot on the Earth’s surface every sixteen 
days. TM observations can generally be obtained with 
a temporal frequency that varies from sixteen days 
to several months depending upon local cloud cover 
conditions. This obviously restricts our ability to 
monitor a wide range of dynamic phenomena that 
produce changes in surface conditions on the 
timescale of hours, days, or weeks, However, it is 
possible to partially overcome this limitation by ob- 
taining random observations of a particular 
phenomenon occurring simultaneously in several dif- 
ferent places (for example, the dispersion of turbidi- 
ty plumes in lakes and rivtas). In certain cir- 
cumstances, such observations can be combined to 
develop generic models of selected hydrologic 
phenomena. 

The spectral range of the TM is also a major factor 
in assessing the utility of TM observations for basic 
hydrological research. The TM is designed to measure 
the intensity of reflected sunlight and thermal radia- 
tion emanating from the Earth’s surface. Variations in 
the intensity of such radiation are generally controll- 
ed by the chemical composition and physical proper- 
ties of the outermost material exposed directly at the 
surface. Earth radiation measured by the TM general- 
ly emanates from the upper few micrometers of the 
cover materials exposed at the surface. Such 
measurements provide little direct information about 
bulk surface properties. However, the spectral range 
of the TM does provide hydrologists with new obser- 
vation capabilities that are directly applicable to 
studies of snow state, water quality, and the thermal 
state of water bodies, TM bands 1-3 and 5 & 7 can 
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specifically be used to discriminate clouds from snow 
and ice. TM band 1 is situated in the blue portion of 
the visible spectrum, and it can be used to measure 
the intensity of radiation which has been transmitted 
through liquid water and reflected from suspended 
materials or submerged surfaces. TM band 1 has also 
proven useful for characterizing the surface proper- 
ties of snow and ice. Finally, TM band 6 is situated 
in the thermal infrared and it can be used to detect 
variations in the surface temperature of lakes, rivers, 
coastal zones, and snow covered areas. The use of the 
TM’s spectral capabilities for the study of specific 
phenomena is discussed in greater detail in a later 
section. 

Finally, the TM has been designed to conduct 
synoptic scale multispectral surveys at a spatial resolu- 
tion of 30 meters (note that the spatial resolution of 
the thermai band TM 6 is 120 meters). The improved 
spatial resolution of the TM can be used to 
characterize the hydrogeomorphology of selected 
areas at a level of detail that has not been possible in 
the past. In addition, this improved resolution can be 
used to characterize the distribution of surface cover 
materials in a more detailed fashion. Information of 
this nature is directly applicable to studies of water 
storage and runoff in different hydrologic 
environments. 

In summary, the measurement capabilities of the 
TM are significantly more sophisticated than those of 
the MSS sensor carried on earlier Landsat spacecraft. 
These advanced measurement capabilities can poten- 
tially be used to study a far wider range of hydrologic 
phenomena. They can specifically be used to monitor 
the surface manifestation of such phenomena on a 
globally repetitive basis. The temporal, spectral, and 
spatial measurement capabilities place important con- 
straints on the nature and frequency of hydrologic 
observations that can actually be performed by the 
TM. Investigators seeking to employ the TM for basic 
research in hydrology must be mindful of these limita- 
tions in designing specific hydrological experiments. 

Prospective Research Themes 

The purpose of the following discussion is to iden- 
tify a series of topical problems in different aspects 
of hydrology that could be profitably addressed 
through the analysis and interpretation of TM data. 
This discussion is not intended to be all-inclusive, 
rather it is intended to illustrate generic types of 
hydrologic problems that are potentially amenable to 
TM analysis. The reader is encouraged to employ the 
following discussion as a point of departure in defin- 
ing specific experiments that make innovative use of 


the TM's measurement capabilities for basic research 
in hydrology. 

Hydrogeomorphic Analysis 

Water catchment areas (l.c. watersheds) consist of 
a network of stream channels that transport running 
water from higher to lower elevations. Watersheds 
represent the spatial transfer function between sur- 
face precipitation and runoff. They evolve in shape 
and character as a function of local geological and 
botanical conditions, regional climatic conditions, and 
the orientation of the watershed with respect to the 
sources and sinks of surface precipitation. The 
improved spatial resolution of the TM provides a 
unique opportunity to characterize the morphology 
of catchment areas at a level of detail that has not been 
possible in the past. TM imagery can specifically be 
used to determine stream order, density, and sinuosity 
at mapping scales of 1:25,000 (see Figure III-13). 
Detailed information of this nature would potentially 
enable hydrologists to define the size and distribution 
of sub-basins within broader catchment areas. Stream 
morphometric characteristics derived from TM im- 
agery could also be analyzed in conjunction with ex- 
isting maps of soil distribution, vegetation, and 
regional rainfall to determine the extent to which 
geological, botanical, and meteorological factors con- 
trol stream geometry and watershed morphology. 

The ability to acquire TM imagery on a near global 
basis also presents an interesting opportunity to study 
watershed development in different climatic settings. 
For example, the occurrence of permafrost terrain and 
seasonally frozen bodies of water at higher latitudes 
may have a significant influence upon the organiza- 
tion of regional drainage systems in polar regions. The 
seasonal release of large quantities of liquid water dur- 
ing periods of spring runoff could also have an im- 
portant influence upon the size and spatial distribu- 
tion of stream networks in high latitude regions. 

Information concerning the regional geomor- 
phology of stream networks situated in many parts 
of the world is quite limited at the present time. This 
is particularly true for many polar, tropical, and arid 
regions that are remote and relatively inaccessible. 
The TM provides a unique capability to explore the 
hydrogeomorphology of such regions at a level of 
detail that has not been achievable in the past. TM 
imagery acquired under different illumination condi- 
tions can specifically be used to determine the occur- 
rence, len^, size, and orientation of stream channels 
at both local and regional scales. Information concern- 
ing stream morphometry derived from TM imagery 
could be used to study the geomorphlc characteristics 
of different regions in a comparative fashion that has 
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Fiffure 

This IS a portion of a Landsat 4 TM imuffe acquired 
on November 16, 1982 over south central Colorado 
near Trinidad, Colorado The major drainaf^e feature 
runninfi west-to east through the center of the image 
is the Purgatory River The mountain range situated 
along the left hand margin of the image is the Sangre 


de Cristos mountains Note the sharp difference 
between the radial drainage pattern associated with 
the two snow-capped peaks in the upper left comer and 
the dendritic drainage pattern associated with the 
Purgatory River The improved spatial resolution of 
the TM could potentially b<! used to study the regional 
geomorphology of selected watersheds at a level of 
detail that has not been possible in the past 


ORIGINAL P' i_ 

OF POOR QUALITY 


ORIGiNAL PAu^ 

OF POOR i . 


whether strear.i network geometry is governed to u 
significunt degree by the si/e ot the catchment area 
(sec, (or example, (iregory and VX'alling, 1975) 


Snow and Ice Studies 

>X'ater in the form o( snow and ice is a major com 
poneni o( the hydrologic cycle Kleven percent o( the 
Karth’s land surface is |xrrmanently covered by 
Kli ciers, whereas more than fifty percent of the land 
surface is blanketed on a seasonal basis by snow 
Snow and ice arc sensitive indicators of climatic 
variability at all scales ranging from regional areas to 
the entire globe In fact, areal variations in snow cover 
that (K'cur on an annual basis are greater than the areal 
variability of surface vegetation or any other type of 
natural cover material Because snow possesses an 
exiremely high albedo, this variability produces large 
spatial and temporal variations in the surface radia 
tion budget Surface deposits of snow and ice modu- 
late regional variations in air temperature, buffer the 
moisture content of the atmosphere, and they can 
insulate underlying materials Some of the critical 


scientific questions regarding the role o( snow and 
ice in the hydrologic cycle have recenily bc’en iden 
tihed by the Polar Kesearch Btiard ol the National 
Research Council (19H5) and Meier (1982) 

In attempting to study the boundary layer 
meteorology ot snow covered areas a variety ol ques 
turns arise concerning how time varying weather con 
ditions affect the projx-rties of the snowpa*. k and vice 
versa The spectral bands of the I M offer some 
specific observational capabilities that can be applied 
to studies of snowpack extent and condition h'or 
example, frozen water is a more effective absorber 
ot solar radiation at wavelengths greater than I t) 
micrometer than liquid water In addition, large ice 
particles are more absorptive than smaller particles at 
these longer wavelengths (l)ozier, 1985) As a conse 
quence of these factors the liquid and frozen droplets 
that make up clouds tend to be much more highly 
reflective than surface snow and ice at wavelengths 
beyond 1 () micrometer TM band S is centered at a 
wavelength of 1 micrometer It can be- used in con 
junction with other TM bands situated at shorter visi 
ble wavelengths to discriminate surface snow from 
clouds in a highly accurate fashion (see Figure III N) 



Figure 111-14 

The two images shown above are segments of a full 
TM imafie of < )wens Valley in southern California This 
imagery was originally acquired during the winter of 
1982 The image on the left was constructed from TM 
hands 2, J, and 4 Snow and cloud covered areas both 
appear bright within the spectral region covered by 
these hands The image on the right was constructed 


from TM hands 2, 5, and Snow reflectance is less 
than that of clouds in TM hands 5 and ’’ Conse- 
quently, snow covered areas apfmir dark in this image 
whereas cloud covered areas remain bright Note that 
the area of extensive cloud cover situated in the upper 
right comer of each image can he clearly discriminated 
from snow covered areas situated on the lejt side of 
both images (Imagery courtesy of Jack Fngvl. 
Hughes — Santa Barbara Research Corporation) 







Several of the other TM spectral bands can be used 
to characterize snow state or condition. Snow reflec- 
tance is a sensitive function of surface grain size in 
TM bands 4 and 7, whereas snow reflectance prop- 
erties are sensitive to surface contamination at shorter 
wavelengths corresponding to TM bands 1,2, and 3. 
Finally, TM band 6 can be used to monitor spatial 
variations in snowpack thermal conditions. In prin- 
ciple, It should be possible to combine Information 
derived from these different spectral bands to form 
new spectral indices of snow condition. Such indices 
could be monitored in repetitive TM surveys to study 
snowpack response to constantly changing 
meteorological conditions. 

The inverse problem is equally Intriguing, namely 
how the snowpack influences or modulates regional 
weather and climate. TM imagery could potentially 
be used to monitor annual variations in snow accumu- 
lation, the altitude, slope, and aspect of intercon- 
nected catchment basins, and the concavity or con- 
vexity of the snowpack itself. These various types of 
observations could be used to test and refine radia- 
tion budget models for high latitude watersheds 


Glaciology 

Glaciers are unique features of high mountain and 
arctic environments that are important in a 
climatological, ecological, and hydrological sense, 
Smaller glaciers react passively to the climate, whereas 
large ice sheets interact with or even control it, On 
intermediate scales, the linkages between glacier oc- 
currence and climatology are not well understood at 
the present time, A large mountain glacier system may 
dominate the local climate; whereas, a large tabular 
Iceberg may carry climatic anomalies with it during 
its drift. At the other extreme, some modeling results 
suggest that the deterioration of global climate to full 
glacial conditions may have trailed the development 
of ice sheets in North America and Scandinavia, 
Because glaciers may be extremely sensitive indicators 
of climatic change, variations in the extent of glaciers 
have been measured sporadically for many centuries, 
Except for some reconnaissance scale estimates, the 
total area of the Earth covered by glaciers, and the 
hydrological impact of glaciers has been largely unde- 
fined. TM observations may help in defining the 
extent and impact of glaciers in selected regions. 

Glacier variations may be monitored by sateiiites 
if small scale glacial features are larger than the instan- 
taneous field-of-view of the orbital sensor system, The 
motion of calving glaciers, surging glaciers, and small 
ice caps has been successfully measured using Land- 
sat MSS imagery in the past (see Figure III-15). How- 


ever, most glacier variations are so slow that annual 
changes in glacier position cannot be seen on MSS 
imagery. Similarly, glacier motion can, in some in- 
stances, be measured but the motion must be rapid 
in order to exceed the Inherent positioning error of 
about 100 meters for the MSS .sensor on Landsats 1-3. 
However, with the higher resolution of the TM and 
smaller positioning error of Landsats 4 and 5, It may 
be possible to measure small motions of glaciers for 
the first time. 

Periodic glacier surges are an Important problem 
and in some areas constitute an economic danger. 
Surging glaciers can be identified and monitored using 
MSS Imagery, It has been suggested that there are 
three classes of glaciers: those that never surge, those 
that surge periodically, and those that surge con- 
tinuously at variable rates of movement depending 
on the available ice flux and the degree of supercavita- 
tion. A variety of theories have been advanced con- 
cerning the causes and dynamics of glacial surges. The 
TM with its seven bands in the visible, near infrared 
and thermal Infrared and with its spatial resolution of 
30 m should provide more Information regarding sur- 
face configuration, ice flow velocities, and ice 
temperatures than was heretofore possible from a 
remote sensing platform and thereby, contribute to 
improved understanding of the dynamics of surging 
glaciers. 



1973 


Figure III-15 


This graph displays estimates of downslope surface 
velocity that were inferred from Landsat Multispectral 
Scanner imagery of the Tweedsmuir Glacier acquired 
during 1973. The surface velocity of the glacier was 
estimated along three lateral transects labelled A, B, 
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muiC in this figure (soo Miller, 1973 for details). Note 
the sharp surge in the surface velocity of this glacier 
that occurred during the spring of 1973- The improved 
spatial resolution of the TM should enable hydrologists 
to estimate the surface motions of glaciers with 
improved precision in the future. 

Snowlines on glaciers have been found to be rela- 
tively easy to identify in MSS imagery. It may be possi- 
ble to monitor changes in the mass of individual 
glaciers in the future by exploiting the improved 
spatial resolution of the TM. The mass balance of a 
glacier is determined by rates of net snow accumula- 
tion or ablation, and it is controlled by precipitation, 
wind, and radiation, TM band 5 would play a critical 
role in such mass balance studies because of the 
prevalence of clouds at higher latitudes during both 
the summer and winter seasons. The synoptic scale 
of TM imagery and Its availability on a near global 
basis also provides an Interesting opportunity to deter- 
mine the degree of correlation between the mass fluc- 
tuations of individual glaciers. Earlier studies of mass 
fluctuations have generally been limited to individual 
glaciers within specific regions. TM Imagery could, 
in principle, be used to determine whether the expan- 
sion of a selected set of glaciers is compensated in a 
broader context by the attenuation of glaciers situated 
in other polar regions. Observations of this nature are 
essential in attempting to forecast the climatic 
significance of temporal fluctuations in the mass of 
individual glaciers or regional glacier fields. 

The properties of surface albedo and temperature 
of a glacier are of primary Interest to energy balance 
Studies, As previously mentioned in the discussion on 
snow, the albedo of frozen water is sensitive to grain 
size in TM bands 4 and 7, whereas it is sensitive to 
surface contamination in TM bands 1, 2, and 3. In 
principle, it should be possible to estimate the extent 
to which ice albedo is degraded by the deposition of 
atmosphere aerosols. The thermal channel on the TM 
(band 6) permits the surface temperature to be mea- 
sured at a spatial resolution of 120 meters that is 
superior to all other orbital sensors. Thus, the TM may 
prove to be useful in developing improved models 
of energy balance in polar areas on both a regional 
and global scale. 

Runoff Estimation and Modeling 

Water plays a critical role in most of the biophysical 
processes that occur on the Earth’s land surface. In 
many environments, the principal source of water for 
such processes is surface precipitation. Hydrologists 
are keenly interested in developing quantitative 


models that will enable them to predict how surface 
precipitation Is distributed and stored in surface 
vegetation and soils (sec Figure III-16). Such models 
can be used to assess the viability of selected biophysi- 
cal processes. They can also be used to study the bio- 
geochemical cycling of water-soluble elements and 
compounds which can be transported over large 
distances by running water. 

Runoff estimation and modeling is closely linked 
to studies of catchment morphology and snow con- 
dition discussed in earlier sections. TM data could 
make a significant contribution to modeling the runoff 
potential of selected regions. TM data can specifically 
be used to map the occurrence and distribution of dif- 
ferent land cover conditions (e.g. bare fields, agricul- 
tural fields, grasslands, forested areas, etc,). They pro- 
vide a means of determining the extent and condition 
of snow cover in mountainous or high latitude areas. 
They can also be used to Identify small depressions, 
ponds, and lakes that may serve as temporary storage 
sites during periods of runoff. All of these various 
types of Information are Important considerations in 
attempting to forecast the response of a particular 
catchment area to surface precipitation. Knowledge 
of the extent and nature of cover materials derived 
from TM observations could presumably be com- 
bined with ancillary Information regarding surface 
relief, soil permeability, etc. to develop models that 
predict the volume, rate, and location of runoff in 
response to specific types of precipitation events. 
Expressed in another way, land cover information 
derived from TM imagery could be used to develop 
distributed parameter models of the surface boundary 
conditions that govern water infiltration and storage. 
Models of this nature could be tested and evaluated 
with respect to more conventional methods of runoff 
estimation (see Klrkby, 1973). Quantitative runoff 
models based upon TM image analysis could. In turn. 



Figure III-16 

A simplified representation of the principal factors 
involved in predicting the runoff of surf ace precipita- 
tion (after Gray, 1970). 
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be coupled with models of aerosol deposition, plant 
respiration, and soil formation to predict the net 
transport of selected elements and compounds 
through a particular watershed that occurs In response 
to different types of precipitation events. 

Discussions of the measurement capabilities of the 
TM have suggested that TM data might lead to signifi- 
cant improvements in modeling the runoff potential 
of urban areas. The Improved spatial resolution of the 
TM is particularly effective in identifying many small 
scale impervious features such as roads and buildings 
that commonly occur in urban landscapes. The TM 
spatial resolution may also prove effective In 
characterizing the spatial variability of surface vegeta- 
tion in urban environments. It may even provide a 
basis for distinguishing different types of urban 
neighborhoods (e.g. Industrial, residential, business, 
recreational parks, etc.) which, In turn, have very dif- 
ferent runoff characteristics. Improved characteriza- 
tion of urban landcover conditions should not only 
lead to better quantification of runoff volume from 
urban watersheds, but it would also Improve our abil- 
ity to forecast water quality variations in neighbor- 
ing rivers and lakes resulting from non-point sources 
of pollution. 

As discussed in other portions of this document, 
the measurement capabilities of the TM can be used 
to monitor temporal variations in the condition and 
extent of land cover materials, For example, temporal 
variations in TM bands 3 and 4 can be used to monitor 
the budding, flushing, and subsequent senescence of 
leafy vegetation in temperate regions. TM bands 4-7 
can be used to monitor water and thermal stresses In 
vegetation canopies. Day and night rime observations 
of natural thermal emission in TM band 6 can be 
analyzed in conjunction with meteorological data to 
monitor variations in the moisture condldon.s and 
water holding capacity of surface soils. Finally, TM 
observations can be used to monitor major changes 
in cover conditions associated with logging activities, 
forest fires, crop harvesting, and man’s activities. All 
of these various types of information could potentially 
be used to re-parameterize existing runoff models in 
a fundamental fashion, and develop a new generation 
of models which predict the time-varying runoff 
potential of selected watersheds as a function of the 
current condition of the surface materials found 
within the catchment area. The solubility of selected 
elements and compounds is particularly sensitive to 
the conditions of cover materials. Consequently, time- 
varying models of runoff potential could lead to sig- 
nificant improvements in our ability to forecast the 
net flux of selected elements and compounds out of 
individual watersheds during different times of the 
year, 


In areas where soils arc exposed directly at the 
Earth's surface, TM data may be used to discriminate 
soils on the basis of their textural characteristics 
which, in turn, may be related to their infiltration 
characteristics. It may be possible to employ the tex- 
tural characteristics of soils, as expre.ssed in multispcc- 
tral TM Imagery, to map variations in soil roughnes-, 
coefficient over extended watersheds. Quantitative 
descriptions of the spatial variability of ,H>il type, tex* 
turc, and roughness properties in relation to the 
regional stream network would presumably contri- 
bute to runoff modeling. TM data may also be used 
to study the nature of runoff phenomena in different 
portions of an extended watershed (e.g. overland 
flow, sheet erosion, detentional flow, depresslonal 
storage, etc.), and to determine the factors that govern 
different runoff regimes. 

Water Quality 

Water quality refers to the nature and concentra- 
tion of materials su.spended or dissolved in water. The 
suspension or dissolution of the materials results from 
the processes of the hydrologic cycle, and the nature 
and concentration of the entrained materials deter- 
mine the suitability of the water for different applica- 
tions (e.g, irrigation and human consumption). At least 
two spaceborne Instruments, the Coastal Zone Color 
Scanner (CZCS) on the Nimbus-7 satellite and the 
Ocean Color Experiment (OCE) on board the second 
flight of the Space Shuttle (STS-2), have been used to 
study the nature and distribution of suspended 
materials within large bodies of water, The TM, 
however, offers new capabilities for the remote sens- 
ing of water quality. In particular, the TM offers a 
much finer spatial resolution than either CZCS or the 
OCE while sensing energy within spectral bands 
similar to bands of the CZCS and the OCE. TM data 
can potentially provide new water quality informa- 
tion applicable to a wide range of hydrologic research 
areas which Include the following topics. A review 
of past efforts using MSS data for water quality studies 
is provided by Salomonson et al., 1983. 

The water quality parameters which seem most 
amenable to remote measurement are the concentra- 
tions of chlorophyll and suspended solids. These 
parameters can be measured by physically sampling 
the water, but concentrations can change quickly and 
acquiring sufficient samples to characterize the 
distribution of chlorophyll and sediments within a 
large water body is often impossible. By virtue of its 
fine spatial resolution, the TM can simultaneously pro- 
vide a large number of spectral reflectance observa- 
tions and an instantaneous synoptic view of a water 
body. Characterization of the spatial variability of 
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water quality or turbidity. a» derived from ifie jpec* water temperature. Thene varKm» parameters can be 
tral data, could siKnIficantly contribute to the directly observed with TM data. Static imanes of cir- 

undersiandinK of «Hhcr hydrologic piienomena suc h culatory (Miierns can certainly aid in the development 

as sediment loading, biological production in water. and testing of water circulation models. TM data ran 

eutrophication, and circulation and mixing also he used to observe the mixing of water from dif* 

lerent sources (see Figure III- 1 7) .Mixing can 
(.In uiMtion Mtid MiMing- Models of water circula- dramatically impact water quality. For instarue. fresh 

tion are necessary to forecast the distribution and water flowing into saline bays and estuaries 

de|>osiiion of suspended materials, and, on a global significantly alters the salinity, temperature, quality, 

scale, for studies of meteorology and climatology . (^ir- and biota of the local aqratk' environment. As another 

culatory patterns are expressed by the distributitms example, the use of large volumes of water to C(m>I 

of suspended sediments, chlorophyll, and surface power plant generators can impact water quality 



Figure 111-17 tides, whereas the darker waters to the north contain 

lower concentrations of suspended particles. TM band 
This is a portion of a TM hand I image of San Fran- I is particularly uell suited for studies of water quality 

cisco Bay that was originally acquired in December, because of the optical transmission properties of water 

1982. The albedo contrast observed in the upper in the blue portion of the visible spectrum Fine scale 
reaches of the hay marks a boundary between tu<o variations that cKCur along the boundary between the 
water masses. The brighter waters to the south con- two water masses are caused by local mixing and they 
tain a significant concentration of suspended silt par- are readily obs,rvable in TM band / imagery. 
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when warm water Is mixed back Into the source of 
the cooling water. Remote obsservations of sediment, 
chlorophyll, and temperature distributions could 
potentially aid In the development of mixing models 
describing the effects of mixing on water quality. 

Biologic Productivity Eutrophication. In 
addition to the remote seihvtg of terrestrial plants, 
the TM spectral bands are appropriate for remote 
observations o" "hlorophylfcontaining algae In rivers, 
lakes, and oceans. Many of the major research Issues 
related to the study of aquatic ecosystems that are 
potentially amenable to TM analysis are analogous to 
research Issues discussed earlier for land-based 
ecosystems containing natural and cultivated vegeta- 
tion. These Issues Include: blogeochemlcal cycles of 
carbon dioxide, nitrogen, phosphorous, and 
nutrients: the energy budget In aquatic ecosystems; 
the primary biologic productivity of aquatic 
ecosystems; and algal growth in response to water 
quality and ecosystem dynamics. 

Algal growth Is controlled by a range of water 
quality factors which include: water chemistry, 
temperature, nutrient cycles, and the energy budget. 
Variations In algal growth indicate changes in water 
quality and observations of changes in algae are thus 
useful in the study of the hydrological processes 
which effect water quality. TM data can potentially 
be used to monitor spatial variations in the nature, 
density, distribution, and temporal succession of algal 
species, The temporal succession of algae in fresh 
water lakes provides a direct Indication of lake 
eutrophication. 

Sediment Loading. Ground water runoff is a por- 
tion of the hydrological cycle which directly impacts 
water quality by carrying suspended and dissolved 
materials into water bodies. TM data can provide 
Information on both the terrestrial features which 
contribute materials to the runoff and the resultant 
sediment load in the wafer. This information can aid 
in the development and testing of runoff models and 
is necessary for characterizing non-point sources of 
water pollution. Characterizing the sources of pollu- 
tion is the first step in the development of strategies 
to prevent water quality degradation. 


Water circulation plays a major role in cycling 
nutrients through coastal wetland areas. Nutrients arc 
assimilated within wetlands during the spring and 
early summer, and they are released In late summer 
and fall. The cycling of nutrients Is closely linked to 
seasonal variations in water circulation. Studies of 
nutrient fluxes in wetlands are difficult to conduct due 
to a lack of sufficient hydrological detail with which 
to help calculate nutrient budgets. The various 
capabilities of TM to penetrate shallow water, sense 
subtle soil moisture differences, discriminate vegetation 
species and density, and assess water quality var'a^les 
can all contribute to a better understanding u . 
hydrologic regime In wetland environments. TM image 
analysis can speclflcally contribute to Improved delinea- 
tion of hydrologic features, saltwater Intrusions, depth 
of small water bodies, eutrophication levels, potential 
sources of thermal loading, and water quality. Studies 
of seasonal water circulation and the Impact of major 
storms on nutrient fluxes represent Important areas for 
future research. 

Water Balance Studies 

The concept of water balance over regions of all 
sizes links together most of the ideas, concepts, or 
research areas discussed above (see again Figure III-1 1). 
To do an effective water balance study the magnitudes 
of storages in lakes, rivers, streams, soil moisture, 
atmospheric moisture, ground water, etc., must be 
evaluated as well as the principal water fluxes such 
as evapotranspiratlon, precipitation, and runoff. 

TM image analysis could contribute to studies of 
water balance through the delineation of large scale 
variations in land cover characteristics, particularly 
those variations associated with climatic change or the 
activities of man. For example, studies quantifying the 
extent of rice planting in Southeast Asia, or the reduc- 
tion of forests in the Amazon region may lead to bet- 
ter estimates of the water balance over those large 
areas through better estimation of evapotranspiratlon 
and runoff, Studies depicting the spatial patterns of 
the water budget over relatively large areas (10“* - 10^ 
km^) would make a very significant contribution to 
our understanding of the global hydrologic cycle (see 
Miller, 1977). 
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